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PREFACE 


This  is  the  8th  Texas  International  Conference  on  Approximation  Theory.  It  is  being 
held  in  College  Station  between  the  8th  and  12th  of  January,  1995,  and  is  dedicated  to 
our  distinguished  colleague  E.  Ward  Cheney  in  honor  of  his  65th  birthday.  We  extend  a 
very  warm  welcome  to  all  participants  and  wish  them  a  pleasant  and  productive  sojourn 
in  ‘Aggieland’. 

A  special  feature  of  this  conference  is  the  awarding  of  the  first  Vasil  A.  Popov  Prize 
in  Approximation  Theory  to  an  outstanding  young  mathematician.  This  prize  will  be 
awarded  every  three  years  in  this  conference  series.  We  thank  Ron  DeVore  for  establishing 
this  prize,  and  for  chairing  the  selection  committee  which  also  comprises  Charles  Chui, 
Paul  Nevai,  Allan  Pinkus,  Pencho  Petrushev,  and  Edward  Saff. 

One  distinctive  aspect  of  this  conference  is  that  a  special  effort  has  been  made  to 
include  several  interdisciplinary  topics  of  current  interest.  In  this  regard,  we  are  most 
grateful  to  the  following  special  session  chairs:  Akram  Aldroubi  (Sessions  9B  and  14B, 
Wavelets),  Chandrajit  Bajaj  (Session  19A,  Scientific  Visualization),  and  Joseph  Warren 
(Session  24A,  Subdivision  Schemes  in  CAGD). 

In  addition  to  the  invited  plenary  lectures,  177  research  presentations  have  been  or¬ 
ganized  in  two  parallel  sessions.  With  a  few  exceptions  (due  to  scheduling  constraints  and 
late  requests  for  time  changes),  these  talks  have  been  grouped  according  to  subject  areas, 
with  little  conflict  between  the  opposing  sessions.  However,  despite  our  best  efforts,  these 
arrangements  are  bound  to  fall  short  of  some  expectations,  and  for  this  we  apologize. 

As  in  the  case  of  the  previous  seven  symposia,  the  National  Science  Foundation  is  the 
major  sponsor  of  this  conference.  Additional  funding  from  the  U.S.  Army  Research  Office^ 
is  also  acknowledged  and  appreciated. 

The  abstracts  printed  in  this  program  booklet  were  either  reproduced  from  e-mail 
messages  or  typed  by  the  conference  secretary,  Ms.  Jan  Want;  no  stylistic  changes  were 
attempted.  The  booklet  itself  was  prepared  by  Margaret  Chui.  It  is  both  om  duty  and 
great  pleasure  to  thank  Jan  and  Margaret  for  their  imtiring  efforts. 


The  Organizers 
Januaiy  1995 


t  The  views,  opinions,  and/or  findings  contained  in  this  program  booklet  are  those  of 
the  author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army  position, 
pohcy,  or  decision,  unless  so  designated  by  other  documentation. 
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sia,  Unconditional  convergence  of 

many.  Algebraic  moduli  of  smooth¬ 

multivariate  D^-splines 

ness  and  best  approximations 

5:50  p.m. 

Vladimir  Vasilenko,  Novosibirsk, 

Antonio  da  Silva,  Rio  de  Janeiro, 

Russia,  Algorithms  for  multivariate 

Brazil,  Jackson-type  theorems  on  the 

spline  approximation 

real  line 

6:10  p.m. 

Thomas  Kunkle,  Charleston,  S  Car¬ 

Burkard  Bolster,  Christchurch,  New 

olina,  Generalized  exponential  box 

Zealand,  Interpolation  and  incidence 

splines 

geometry 

6:30  p.m. 

Yuesheng  Xu,  North  Dakota  State 

Theodore  Kilgore,  Auburn  Univ,, 

Univ,,  Construction  of  wavelets  on  fi¬ 
nite  domains 

Some  inequalities  for  derivatives 

6:50  p.m. 

Guanglu  Zhang,  Dongying,  China, 

Frank  Pinter,  Leuven,  Belgium,  Per¬ 

Integral  representation  of  bivariate 

turbation  of  orthogonal  polynomials 

splines 

on  an  arc 
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Wednesday,  January  11 

Session  16A 

Session  16B 

Chair:  Yuan  Xu 

Chair:  Jurgen  Prestin 

8:20  a.m. 

Knut  Petras,  Braunschweig,  Ger- 

Tian-Xiao  He,  Illinois  Wesleyan 

many,  Polynomial  interpolation  and 

Univ,,  Spline  interpolation  and  its 

Peano  kernels 

wavelet  analysis 

8:40  a.m. 

Ricardo  Estrada,  Texas  A&M  Univ,, 

Tom  Hogan,  U  of  Wisconsin,  Madi- 

Asymptotic  approximation  by  cardi- 

son.  Shifts  of  multivariate  refinable 

nal  series 

functions 

9:00  a.m. 

Armin  Iske,  Gottingen,  Germany, 

Suresh  Lodha,  U  of  California, 

Title  to  be  announced 

Santa  Cruz,  B6zier  wavelets 

9:20  a.m. 

Thomas  Sauer,  Erlangen,  Germany, 

Gerlind  Plonka,  Rostock,  Germany, 

Newton  method  for  polynomial  in¬ 
terpolation 

Multiwavelets  with  short  support 

9:40  a.m. 

Jeremy  Levesley,  Leicester,  England, 

Gina  Puccio,  Messina,  Italy,  Recog- 

Generalised  sh-spline  interpolants  on 

nition  of  biological  shapes  by  wavelet 

compact  abelian  groups 

analysis 

10:00  a.m. 

Johnnie  Baker,  Kent  State  Univ,, 

Mei  KobayashJ,  Tokyo,  Japan, 

Ohio,  Sorting  algorithm  for  reconfig- 
urable  meshes  and  networks 

Wavelet  analysis  of  speech  signals 

Session  17 

Chair:  George  Loren  tz 

11:00  a.m. 

Gilles  Pisier,  Texas  A&M  Univ,  and  Univ,  of  Paris  VI,  Real  and  complex 

interpolation  of  operator  spaces 

Session  18 

Chair:  Bjorn  Jawerth 

1:30  p.m. 

Ingrid  Daubechies,  Princeton  Univ.,  Determination  of  regularity  ofrefinable 
functions 

Scientific  Visualization  Session 

Session  19B 

Chair:  Chandrajit  Bajaj 

Chair:  Christopher  Heil 

2:30  p.m. 

B.K.  Natarajan,  Hewlett  Packard, 

Kurt  Jetter,  Duisburg,  Germany, 

Sparse  approximate  multiquadric  in- 

de  Boor’s  stability  result  and  sym- 

terpolation 

metric  preconditioning 

2:50  p.m. 

Tom  Foley,  Arizona  State  Univ., 

Joachim  Stockier,  Duisburg,  Ger- 

Tempo,  Plane  preserving  volume  de- 

many,  Wavelet  analysis  of  scattered 

formations  using  scattered  landmark 
points 

data 

3:10  p.m. 

Marian  Neamtu,  Vanderbilt  Univ., 

Fran  Narcowich,  Texas  A&M  Univ., 

Spherical  splines 

Nonstationary  spherical  wavelets  for 
scattered  data 

3:30  p.m. 

Chandrajit  Bajaj,  Purdue  Univ., 

Joe  Ward,  Texas  A&M  Univ.,  Wave- 

Modeling  and  visualizing  scat- 

lets  associated  with  periodic  basis 

tered  function  data 

functions 

Wednesday,  January  11  (Cont’d) 


Session  20A 

Session  20B 

Chair:  Xiang  Ming  Yu 

Chair:  Ewald  Quak 

4:30  p.m. 

Dany  Leviatan,  Tel  Aviv,  Israel, 

Philip  Smith,  WindWard  Technolo- 

Convex  polynomial  approximation  in 

gies.  Approximation  on  compact  sets 

the  Lp  norm 

by  radi^  basis  functions 

4:50  p.m. 

George  Anastassiou,  U  of  Memphis, 

Hein  Hundal,  Penn  State  Univ., 

Differentiated  shift-invariant  integral 

Generalizations  of  Dykstra’s  algo- 

operators 

rithm 

5:10  p.m. 

Mi  Zhou,  U  of  Memphis,  Asymptotic 

Jiinjiang  Lei,  Southern  Illinois 

expansions  for  certain  operator  semi- 

Univ,,  Integral  operators  in  shift- 

groups 

invariant  spaces 

5:30  p.m. 

Elisabetta  Santi,  Rome,  Italy,  Nu- 

Gisela  Poplau,  Rostock,  Germany, 

merical  solution  of  singular  integral 

Fast  Poisson  solvers  in  shift-invariant 

equations  of  Cauchy  type 

spaces 

5:50  p.m. 

Kurt  Riedel,  Courant  Institute,  New 

Shayne  Waldron,  U  of  Wisconsin- 

York,  Piecewise  convex  function  es¬ 

Madison,  Error  represention  in  mul¬ 

timation  and  model  selection 

tivariate  polynomial  interpolation 

6:10  p.m. 

Alex  Sidorenko,  Courant  Institute, 

Tianping  Chen,  Fudan  Univ.,  China, 

New  York,  Adaptive  kernel  estima¬ 

Radial  basis  functions  and  applica¬ 

tion  of  spectral  density 

tions  to  neural  networks 

6:30  p.m. 

Lefan  Zhong,  Beijing,  China,  Inter¬ 

Vladimir  Andrievskii,  Eichstdtt, 

polation  bases,  the  Marcinkiewicz- 

Germany,  Uniform  polynomial 

Zygmund  inequality  and  Ap  weights 

approximation  on  a  quasi  disk 

6:50  p.m. 

M.A.  Bokbari,  Dhahran,  Saudi  Ara¬ 

D,V.  Pai,  Bombay,  India,  Continu¬ 

bia,  L2  and  Loo-approximation  with 

ity  of  the  prox  and  restricted  center 

interpolatory  constraints 

maps 

Thursday,  January  12 


8:20  a.m. 

8:40  a.in. 

9:00  a.m. 

9:20  a.m. 

9:40  a.m. 

10:00  a.m. 

Session  21A 

Chair:  N.  Sivakumar 

Guanrong  Chen,  U  of  Houston,  Al¬ 
gorithms  for  matrix-valued  Nevan- 
linnarPick  interpolation 

X.-L.  Liu,  Vanderbilt  Univ,,  B4zier 
triangle  patches  on  sphere-like  sur¬ 
faces 

Sonya  Stanley,  Vanderbilt  Univ,,  Cir¬ 
cular  Bemstein-B4zier  curves  and 
degree  raising 

LJ.  Anderson,  Huddersfield,  Eng¬ 
land,  Matrix  methods  for  approxima¬ 
tion  of  multivariate  data 

Ranjan  Muttiah,  Texas  A&M  Univ., 
Climate  surfaces  using  thin  plate 
splines 

G.A.  Watson,  Dundee,  Scotland,  An 
approximation  problem  in  measure¬ 
ment  analysis 

Session  2  IB 

Chair:  Laura  Montefusco 

Jian-ao  Lian,  Prairie  View  A&M 
Univ.,  Texas,  Characterization  of 
polynomial  reproduction  for  multi¬ 
scaling  functions 

Alexander  Petukhov,  St.  Petersburg, 
Russia,  Nonorthogonal  compactly 
supported  wavelets 

BL  Sendov,  Sofia,  Bulgaria,  Com¬ 
pression  through  rotations 

C.  Guerrini,  Bologna,  Italy,  Paral¬ 
lel  wavelet-Galerkin  methods  using 
wavelet  packets 

Gerhard  Opfer,  Hamburg,  Germany, 
Cardinality  of  orthogonal  wavelet 
multiplier  matrices 

Karsten  Urban,  Aachen,  Germany, 
Divergence- free  wavelet  bases 

Session  22 

Chair:  Tom  Lyche 

11:00  a.m. 

Rong-Qing  Jia,  Edmonton,  Canada,  Refinable  shift-invariant  spaces:  from 
splines  to  wavelets 

Session  23 

Chair:  Will  Light 

1:30  p.m. 

Robert  Schaback,  Gottingen,  Germany,  Multivariate  interpolation  and 
approximation 

Subdivision  Scheme  Session 

Session  24B 

Chair:  Joe  Warren 

Chair:  Joe  Lakey 

2:30  p.m. 

Nta  Dyn,  Tel  Aviv,  Israel,  Non- 

Peter  Oswald,  Texas  A&M  Univ., 

stationary  subdivision  schemes  and 

Piecewise  linear  prewavelets  with 

multiresolution  analysis 

small  support 

2:50  p.m. 

Ulrich  Reif,  Stuttgart,  Germany,  Es- 

Johan  de  Villiers,  Stellenbosch,  S 

timation  of  curvature  continuous 

Africa,  Construction  of  fundamental 

subdivision  surfaces 

interpolants  and  wavelets 

3:10  p.m. 

Leif  Kobbelt,  Karlsruhe,  Germany, 

Laura  Montefusco,  Bologna,  Italy, 

A  variational  approach  to  subdivi- 

Adapted  wavelet  packets  for  paral- 

sion 

lelism  in  image  restoration 

3:30  p.m. 

Joe  Warren,  Rice  Univ.,  Texas, 

Sonia  Gomes,  Campinas,  Brasil, 

Smooth  functional  subdivision  meth- 

Numerical  algorithms  for  periodic 

ods  for  irregular  triangulations 

biorthogonal  wavelets 

4:30  p.m. 

4:50  p.in. 

5:10  p.in. 

5:30  p.m. 

5:50  p.m. 

6:10  p.m. 

6:30  p.m. 

6:50  p.m. 


Thursday,  January  12  (ContM) 


Session  25A 
Chair:  Alistair  Watson 

FYancois  Dubeau,  Sherbrooke,  Que¬ 
bec,  Interlacing  properties  of  roots  of 
Euler-Frobenius  polynomials 

J,  Savoie,  Richelain,  Quebec,  Opti¬ 
mal  error  bounds  for  spline  interpo¬ 
lations 

J.  Koch,  Stuttgart,  Germany,  Geo¬ 
metric  Hermite  interpolation 

Barbara  Bertram,  Michigan  Techno¬ 
logical  Univ.,  Divergent  Neumann  se¬ 
quences  for  Fredholm  integral  equa¬ 
tions 

V.V.  Zhuk,  SL  Petersburg,  Russia, 
Approximation  of  continuous  func¬ 
tions  of  bounded  variations 

Shusheng  Xu,  Jiangsu,  China,  Char¬ 
acterization  of  best  Lp  approxima¬ 
tion  with  restricted  derivatives 

Avraham  Melkman,  Beersheba,  Is¬ 
rael,  An  efficient  variation  on  the 
Oslo  algorithm 

Xiaolei  Dong,  Harbin,  China,  Best 
approximation  on  function  spaces 


Session  25B 
Chair:  Joe  Ward 

M.  Skopina,  St.  Petersburg,  Russia, 
Pointwise  convergence  of  periodized 
wavelet  series 

Shuzhan  Xu,  Edmonton,  Cana¬ 
da,  Computing  continuous  wavelet 
transforms  by  subdivision 

Valery  A.  Zheludev,  St.  Petersburg, 
Russia,  On  signal  processing  using 
periodic  spline  wavelets 

Chun  Li,  Texas  A&M  Univ.,  Mul¬ 
tivariate  interpolating  wavelets  and 
their  duals 

Jianrong  Wang,  HARC,  Texas, 
n-multiple  wavelets  with  n-multiple 
vanishing  moments 

S.  Ya.  Novikov,  Samara,  Russia, 
Canonical  injection  operators  for  re¬ 
arrangement  invariant  spaces 

Vadim  Iv.  Filippov,  Saratov,  Russia, 
Classification  of  function  systems  in 
Lp  and  wavelets 
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ABSTRACTS 


Oblique  Projections  in  Unitary  Invariemt  Spaces 
Akram  Aldroubi 

National  Institutes  of  Health,  Biomedical  Engineering 
and  Instrumentation  Program 
Bethesda,  Maryland 

Let  be  a  Hilbert  space,  and  let  O  be  a  unitary  operator  on  H.  Using  the  operator 
O,  and  a  set  of  vectors  in  H  we  construct  an  0-invariant  subspace  U  CH: 

=  c‘€/2,  i=i,2 . n\.  (1) 

li=ifcez  ) 

We  give  the  necessary  and  sufficient  conditions  for  U  to  be  a  well-defined  subspace  of  Ti. 
with  i-i2  Riesz  basis.  We  then  consider  the  obUque  projection  Pu±v 

on  the  space  U(P,{tp^}i=i,...N)  in  a  direction  orthogonal  to  V{R,  {T*}i=i,...7v)-  We  give 
necessary  and  sufficient  conditions  on  O,  R,  for  Pu±v  to  be  well 

defined.  We  then  use  the  results  to  construct  the  theory  of  biorthogonal-multi-wavelets 
which  is  the  generalization  of  the  Goodman,  Lee,  and  Tang  theory  of  wavelets  in  wandering 
subspaces. 


Differentiated  Shift-Invariant 
Integral  Operators,  Univariate  Case 

George  A.  Anastassiou 

The  University  of  Memphis 
Memphis,  Tennessee 

In  a  recent  paper  the  author,  along  with  H.  Gonska,  introduced  some  wavelet  type  inte¬ 
gral  operators  over  the  whole  real  line  and  studied  their  properties  such  as  shift-invariance, 
global  smoothness  preservation,  convergence  to  the  unit,  and  preservation  of  probability 
distribution  functions.  These  operators  are  very  general  and  they  are  introduced  through 
a  convolution-like  iteration  of  another  general  operator  with  a  scaling  type  function.  In 
this  paper  the  author  provides  sufficient  conditions,  so  that  the  derivatives  of  the  above 
operators  enjoy  the  same  nice  properties  as  their  originals.  A  sufficient  condition  is  also 
given  so  that  the  “global  smoothness  preservation”  related  inequality  becomes  sharp.  At 
the  end  several  applications  are  given,  where  the  derivatives  of  the  very  general  specialized 
operators  are  shown  to  fulfill  all  the  above  properties.  In  particular  it  is  shown  that  they 
preserve  continuous  probability  density  functions. 
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Efficient  Matrix  Methods  for  the  True  Least-Squares 
Approximation  of  Structured  Multivariate  Data 

I.  J.  Anderson*  and  J.  C.  Mason 

University  of  Huddersfield 
Huddersfield,  U.K. 

It  is  well-known  that  very  eflScient  algorithms  may  be  derived  for  approximating  data 
whose  abscissae  lie  on  a  mesh,  based  on  the  idea  of  separating  the  variables.  The  present 
authors  and  co-workers  have  extended  this  separation-of- variables  approach  in  a  number  of 
papers  and  developed  algorithms  for  other  data  sets,  which  have  structure  but  to  a  lesser 
extent,  such  as  data  on  or  near  to  a  set  of  lines,  data  on  a  family  of  curves,  and  data  which 
are  imiformly  scattered.  Generally  these  algorithms  do  not  produce  “true”  least-squares 
approximations  but  rather  “near-best”  approximations.  However,  an  alternative  approach 
to  such  problems  is  to  set  up  the  complete  overdetermined  system  of  observation  equations, 
without  separating  variables,  and  then  to  exploit  the  special  structures  that  exist  in  the 
matrix,  in  order  to  perform  a  rapid  QR  decomposition.  In  this  paper  we  illustrate  this 
approach  for  data  on  a  family  of  lines,  showing  that  an  efficient  method  can  be  devised 
which  yields  the  true  least-squares  solution.  The  method  is  valid  for  a  wide  variety  of  basis 
functions,  but  can  be  further  refined  in  the  case  of  spline  functions,  to  take  account  of  local 
support  featmres. 


Uniform  Polynomial  Approximation  on  a  Quasidisk 
Vladimir  Andrievskii 
Katholische  Universitat  Eichstatt 
Eichstatt,  Germany 

This  report  is  connected  with  the  study  of  the  values  En{f,G),  n  =  0,1,2...,  of 
best  uniform  polynomial  approximations  of  a  function  /  analytic  in  a  bounded  Jordan 
domain  G  C  C  and  continuous  on  its  closure  G.  The  rate  of  decrease  of  En{f,G)  as 
72  — >  oo,  the  geometric  structure  of  the  boundary  dG  of  G,  and  the  smoothness  of  the 
function  /  near  the  boimdary  interact  in  a  complicated  way.  The  main  subject  of  our  talk 
is  the  consideration  of  the  following  two  problems.  Let  >  0,  be  a  so-called  normal 

majorant  (for  example,  n(6)  =  6^,  c=  const  >  0). 

PROBLEM  A.  Describe  all  functions  /  satisfying 

En{f,G)  =  0{n{l/n))  as  72  — >  oo. 

PROBLEM  B.  Describe  all  fimctions  /  for  which 

Enif,  G)  ~  /2(l/n),  72  =  1, 2, ... , 

where  the  symbol  g  (f  means  that 

1/c  <g>f<c 
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holds  for  some  constant  c>  0. 

We  give  the  solution  of  Problems  A  and  B  in  the  case  when  G  is  an  arbitrary  quasidisk 
and  apply  then  these  results  to  the  study  of  a  problem  of  Turan  concerning  the  correlation 
between  polynomial  and  rational  approximations  on  the  unit  disk. 


Multicriteria  Optimization 

Alexis  Bacopoulos 

University  of  Montreal 
Montreal,  Canada 

We  define  a  convenient  framework  of  vector  norms  and  give  some  results  in  multicri¬ 
teria  optimization.  Theorems  of  existence,  characterizations,  imiqueness  and  computation 
of  best  approximations  are  given  here  in  three  mutually  related  contexts  of  optimality. 
Parallel  computation  and  numerical  schemata  of  optimization  are  discussed. 


Error  Of  An  Arbitrary  Order  For  The  Approximate  Solution 
Of  System  Of  nth  Order  Differential  Equations  With  Spline  Functions 

Badr  Saad  badr*,  Thaxwat  Fawzy,  and  Z.  Ramadan 

University  of  Tanta,  Tanta,  Eg5q)t 
University  of  Suez  Canal,  Issmailya,  Egypt 
University  of  Ein  Shams,  Cairo,  Egypt 

In  this  p)aper  we  present  a  method  for  approximating  the  solution  of  a  system  of 
nonlinear  nth  order  differential  equations  =  fi{x,y,z),z^'^\x)  =  f2{x,y,z)  with 

y^^\xo)  =  yo\z^^xo)  =  Zq\  j  =  0, 1,  ...,n—  1.  We  use  spline  functions  for  finding 
the  approximate  solution.  It  is  a  one-step  method  where  the  error  in  is  of  order 

O  provided  that  /i,  /2  €  C""  ([0, 1]  x  B?)  and  the  modulus  of  continuity  of 

2^”^  is  O  (/i“),where  r  G  iV,  i  =  0, 1, ...,  n-f  r  ,  0  <  a  <  1  and  m  >  1  is  an  arbitrary  positive 
integer  which  equals  the  number  of  iterations  used  in  computing  the  spline  functions  defined 
in  the  method.  It  is  also  shown  that  the  method  is  stable. 


Modeling  and  Visualizing  Scattered  function  Data  on  Curved  Surfaces 

C.  Bajaj*  and  G.  Xu 

Purdue  University,  West  Lafayette,  Indiana,  and 
Lab  of  Scientific  Sz  Engineering  Computing,  Beijing,  China 

We  shall  present  algorithms  for  constructing  a  cubic  or  a  quintic  piecewise 
trivariate  Bemstein-Bezier  polynomial  approximation  F  of  scattered  function  data  defined 
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on  an  arbitrary  smooth  domain  surface  D  in  three  dimensional  real  space  (IR.^)-  The 
smooth  polynomial  pieces  or  finite  elements  of  F  are  defined  on  a  three  dimensional  trian¬ 
gulation  called  the  simplicial  hull  surrounding  the  domain  surface  D.  We  shall  also  present 
techniques  for  visualizing  the  isocontours  of  the  function  F.  The  smooth  polynomial  ap¬ 
proximation  F  allows  one  to  additionally  manipulate  the  local  geometry  of  the  modeled 
function  and  supports  applications  in  data  visualization,  navigation  and  querying. 


A  Constant  Time  Sorting  Algorithm  for  the  3-D 
Reconfigurable  Mesh  and  the  Reconfigurable  Network 

Johnnie  W.  Baker*  and  Mark  S.  Merry 

Kent  State  University 
Kent,  Ohio 

Sorting  techniques  have  numerous  applications.  Current  real  number  and  integer  sort¬ 
ing  techniques  for  the  reconfigurable  mesh  operate  in  constant  time  using  a  reconfigurable 
mesh  of  size  n  x  n  to  sort  N  numbers.  This  paper  presents  a  constant  time  algorithm  to 
sort  n  items  on  a  reconfigmable  network  with  0{y/n  x  s/n  x  y/n)  switches  and  0{y/n  x  y/n) 
processors.  Also,  new  constant  time  selection  and  compression  algorithms  are  obtained. 
All  results  may  also  be  implemented  on  the  3-D  reconfigurable  mesh. 


On  Some  Constants  in  Simultaneous  Approximation 

Katherine  Balazs 
Auburn,  Alabama 

Pointwise  estimates  for  the  error  which  is  feasible  in  simultaneous  approximation  of 
a  function  and  its  derivatives  by  an  algebraic  polynomial  were  originally  pursued  from 
theoretical  motivations,  which  did  not  immediately  require  the  estimation  of  the  constants 
in  such  results.  However,  recent  numerical  experimentation  with  traditional  techniques  of 
approximation  such  as  Lagrange  interpolation,  slightly  modified  by  additional  interpolation 
of  derivatives  at  ±1,  shows  that  rapid  convergence  of  an  approximating  polynomial  to 
a  function  and  of  some  derivatives  to  the  derivatives  of  the  function  is  often  easy  to 
achieve.  The  new  techniques  are  theoretically  based  upon  older  results  about  feasibility, 
contained  in  work  of  Trigub,  Gopengauz,  Telyakovskii,  and  others,  giving  new  relevance 
to  the  investigation  of  constants  in  these  older  results.  We  begin  this  investigation  here. 
Helpful  in  obtaining  estimates  for  some  of  the  constants  is  a  new  identity  for  the  derivative 
of  a  trigonometric  polynomial,  based  on  a  well  known  identity  of  M.  Riesz.  One  of  our 
results  is  a  new  proof  of  a  theorem  of  Gopengauz  which  reduces  the  problem  of  estimating 
the  constant  there  to  the  question  of  estimating  the  constant  in  a  simpler  theorem  of  Trigub 
used  in  the  proof. 
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Haar  Theory  in  Co(T,R'^) 

Martin  Bartelt*  and  Wu  Li 

Christopher  Newport  University 
Newport  News,  Virginia 

In  1956  S.  B.  Stechkin  and  S.  I.  Zukhovickii  characterized  those  finite  dimensional 
subspaces  of  Co{T,R'‘)  from  which  best  uniform  approximations  are  unique.  Using  these 
subspaces  we  extend  the  Haar  Theory  involving  strong  imicity  and  Lipschitz  continuity 
fromCW  toCo(T,R^). 


Compactly  Supported  Radial  Basis  Functions 
B.  J.  C.  Baxter 

Imperial  College  of  Science,  Technology  and  Medicine 
London,  England 

A  radial  basis  function  approximation  has  the  form 

n 

*:=! 

where  (piTZ^  — ►  72.  is  some  given  (usually  radially  symmetric)  function,  (a*;)”  are  real 
coefficients,  and  the  centres  {bk)i  are  points  in  We  investigate  some  properties  of 
when  it  has  the  distributional  Fourier  transform  =  c||^||“^'”,  where  c  is  a  real  number 
and  m  is  a  positive  integer.  The  class  of  such  functions  includes  the  Euclidean  norm 
(p{x)  =  ||x||  when  the  dimension  d  is  odd  and  the  thin  plate  spline  (p{x)  =  |lx|p  log  ||x|| 
when  d=  2.  Our  key  observation  is  that  it  is  easy  to  construct  a  signed  measure  fj,  on 
for  which  the  function  is  an  entire  function  of  exponential  type.  Hence  the  Paley- 

Wiener  theorem  allows  us  to  conclude  that  the  convolution  fi*(p  is  a,  compactly  supported 
function.  One  consequence  of  this  result  is  a  local  stability  estimate  in  the  uniform  norm 
for  the  function  s  defined  by  (1);  If  \\bj  —  6fc||  >  1  for  j  ^  k,  then 

|aj|  <  max{|s(x)|  :  ||x  -  bj\\  <  l}/A2m-i,  j  =  1,  •  •  •  i 

where  A2m-i  denotes  the  distance  in  C7[-l,  1]  of  the  univariate  function  {(p{tu)  :  t  e  72} 
from  the  space  of  polynomials  of  degree  2m—  1,  u  being  any  unit  vector  in  72*^.  Another  is  a 
nonn  estimation  technique  generalizing  earlier  research.  Further,  we  find  that  it  is  possible 
to  choose  the  coefficients  and  the  centres  (5*)*=!  such  that  the  spherical  average 

As  of  s  is  compactly  supported,  where  the  spherical  average  of  a  continuous  function 
/  :  72^^  — >  72  is  defined  by  the  equation 

Af{x)=  f  f{Ux)da{U),  xeR!^.  (2) 

JOd 
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Here  Od  denotes  the  group  of  real  d  x  d  orthogonal  matrices  and  cxd  denotes  the  group 
invariant  probability  measure  on  Od- 


THIN:  A  Collection  of  Routines  for  Fast  Computation 
With  Radial  Basis  Functions 

Rick  Beatson 

University  of  Canterbury 
Christchurch,  New  Zealand 

If  radial  basis  functions  are  used  in  a  direct  way  then  the  demands  in  time  and  storage 
can  be  excessive.  There  are  several  ways  to  reduce  these  demands.  This  talk  concerns  a 
collection  of  routines  implementing  some  of  these  fast  algorithms.  Initially  routines  for  fast 
evaluation,  and  fast  solution  of  the  interpolation  equations  will  be  included.  It  is  hoped  to 
add  further  routines,  such  as  fast  GCV,  at  a  later  date. 


On  a  Problem  of  Approximation  to  the  Analytic  Function 
by  Analytic  Complex  Planar  Splines  in  the  Quasidisk 

V.  Belyi 

Institute  for  Applied  Mathematics  k  Mechanics 
National  Ukrainian  Academy  of  Science 
Donetsk,  Ukraine 

In  fact,  the  problem  of  recovering  of  analytic  function  in  a  quasidisk  by  complex 
planar  splines  is  a  problem  of  spline  approximation.  Here  we  consider  the  approximation 
of  functions  which  are  analytic  in  the  bounded  quasidisk  G  and  continuous  over  its  closure 
G.  Let  San  (/;  be  an  analytic  complex  spline  defined  on  triangular  or  rectangular  grids, 

G 

where  SAsiC)  ^  fbe  complex  planar  spline  of  G.  Opfer  M.  L.  Puri,  and  y{z)  is  the 
quasiconformal  refiection  of  the  complex  plane  with  respect  to  L  =  dG.  Our  main  result 

Wfiz)  -  SaA/-,  ^)IIc(G)  <  MAf-,  hN),  (2) 

where  a;(/;  t)  is  the  modulus  of  continuity  of  /,  Hn  the  mesh  width  and  M  depends  on  G 
and  /  .  The  estimation  (2)  was  known  earlier  for  the  case  of  the  quasismooth  boundary 
dG  only. 

Some  other  results  concerning  problems  of  approximation  by  analytic  splines  are  con¬ 
sidered. 
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On  Riesz-summability  for  Multivariate  Inverse  Fourier  Integrals 

Hubert  Berens 

University  of  Erlangen-Niirnberg 
Erlangen,  Germany 

This  is  joint  work  with  Yuan  Xu.  Motivated  by  Yuan’s  investigations  of  multivariate 
orthogonal  polynomials,  we  were  led  to  study,  what  we  call,  summability  of  inverse 
Fourier  integrals  and  Fourier  seri^. 

Let  /  e  L^(R‘^),  d  G  N,  and  let  /  be  its  Fourier  integral.  We  denote  by 

f  dv,  xeR‘^  and  >  0 

the  £-1  Dirichlet-kernel  and  by 

f  e‘'’*/(v)(fv  =  /*£>«, a(l)(x) 

the  inverse  Fourier  integral  of  /  -  note  that  we  are  integrating  over  the  ^i-ball  in  R*^ 
centered  at  the  origin  with  radius  R>  0.  As  a  central  result  we  proved:  Let  6  >  0.  If 

^  ^  I  ~  '^)  ^  ^  and  R>  0; 

denote  the  Riesz  (R,  S)  means  of  the  inverse  Fourier  integral  of  /,  we  have 

(1)  the  means  define  an  approximate  identity  on  L^(R‘^)  (and  consequently  on  Co(R^) 
and  LP(R‘*),  1  <  p  <  oo));  i.  e., 

V/  G  L^CR*^)  ^li^  •)  =  /  in  norm. 

(2) 

V/  G  L^’CR'')  l4?/5;d(/;  I  <  constd;^  M(/;  x)  a.  e., 

where  M(/;  •)  is  the  maximal  function  of  /;  in  particular, 

V/  €  LP(R'‘)  lim  x)  =  f(x)  a.  e. 

it— >00  *  * 

This  is  in  sharp  contrast  to  the  well-known  Bochner-Riesz  means  which  have  a  critical 
index  (of  convergence),  namely,  6  >  (d—l)/2. 

The  proof  depends  strongly  on  a  closed  form  of  the  Dirichlet  kernel  as  a  divided  difference  of 
a  function  on  11+  with  knots  being  the  square  of  the  coefficients  of  x  =  (xi, . . . ,  xj)  e  R'^. 
There  are  analogous  results  for  Fourier  series. 


Numerical  Transformations  of  Divergent  Neumann  Sequences 
for  EVedholm  Integral  Equations  of  the  Second  Kind 

Barbara  S.  Bertram*  and  Otto  G.  Ruehr 

Michigan  Technological  University 
Houghton,  Michigan 

We  consider  the  numerical  solution  of  second  kind  Fredholm  integral  equations  of  the 
form 

f{x)  =g(x)  + 

where  the  underlying  integral  operator  admits  discrete  spectra  but  has  sup  norm  greater 
than  (or  equal  to)  one.  In  general,  the  application  of  successive  approximations  (Neumann 
series)  is  precluded  because  it  is  divergent.  Abbreviating  (l)as/  =  y  +  Lf,  we  trans¬ 
form  the  Neumann  sequence,  /„  =  X)”=o  ^9^  using  Aitken’s  method  applied  pointwise. 
Somewhat  surprisingly,  we  find  in  several  examples  that  the  transformed  sequence  con¬ 
verges  rapidly  to  the  solution  /.  Since  Aitken’s  method  does  not  usually  sum  divergent 
sequences,  the  numerical  results  imply  that  Neumann  sequences  have  a  simple  structure. 
Indeed,  a  brief  formal  analysis  shows  that 


/ 


k{x,t)f(t),  a<x<b, 


fn-f  -I-  y^AjAj”, 

*>i 

where  Ai  are  the  eigenvalues  of  L  and  A,  depends  on  the  independent  variable  x,  but  not 
on  n. 

Computations  were  done  for  the  following  kernels  on  the  interval  [0,1]:  exp{cx(x  —t)) 
(M.  Kac),  (1  -  xt)-°‘  (M.  L.  Glasser),  exp{axt)  (H.  Brakhage),  and  \x  -  t|”“  (weakly 
singular).  By  varying  a,  the  number  of  eigenvalues  larger  than  1  can  be  controlled.  In  each 
case,  when  o:  was  chosen  to  make  the  Neumann  sequence  diverge,  one  or  more  successive 
Aitken  transforms  effectively  summed  it. 


A  Note  on  Irving  Glicksberg’s  Pseudocompactness  Paper 

Jorg  Blatter*  and  Heinz  Konig 

Universidade  Federal  do  Rio  de  Janeiro 
Rio  de  Janeiro,  Brazil 
and 

Universitat  des  Saarlandes 
Saarbriicken,  Germany 

Dedicated  to  the  Memory  of  Galen  L.  Seever 

Let  X  be  any  non- void  topological  space  and  denote  by  C{X)  the  set  of  all  real- valued 
continuous  functions  on  X  and  by  C*(X)  the  subset  of  aU  bounded  real- valued  continuous 
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functions  on  X.  Irving  Glicksberg  [1],  [2]  solved  the  question  of  when  a  particular  non¬ 
negative  linear  functional  I  on  C{X)  can  be  represented  by  a  Baire-measure  /i  on  A"  in 
the  sense  that 

Hf)  =  f  fdfi  V/  €  CiX) 

Jx 

by  proving  a  pioneering  version  of  what  is  known  today  as  the  Daniell-Stone  Theorem  and 
which  gives  a  necessary  and  sufficient  condition  on  I  in  order  to  have  such  representation. 
He  then  went  on  to  investigate  the  question  of  when  all  non-negative  linear  functionals 
on  C{X)  have  such  representation  and  found  that  a  necessary  and  sufficient  condition 
for  this  is  that  X  be  pseudo-compact,  which  is  to  say  that  C{X)  =  C*{X).  Finally, 
Glicksberg  established,  with  rather  complicated  proofe,  that  pseudo-compactness  of  A  is 
also  equivalent  to  the  validity  of  various  other  classical  theorems  on  C(X),  including  Dini’s 
and  Ascoli’s. 

What  we  do  here  is 

1.  Show  that  pseudo-compactness  of  A  is  equivalent  to  the  validity  of  three  more  classical 
theorems  on  C'(A),  among  these,  surprisingly,  Arzela’s  Theorem,  apparently  much 
stronger  than  Dini’s;  and 

2.  give  completely  elementary  proofs  of  both  Glicksberg’s  results  and  our  own,  of  course. 

References 

[1]  Irving  Glicksberg,  The  representation  of  functionals  by  integrals,  Duke  Math.  J.  19 
(1952),  253-261. 

[2]  Irving  Glicksberg,  Stone-Qech  compactifications  of  products,  Trans.  Amer.  Math.  Soc. 
90,  369-382. 


Convergence  Problem  for  Kergin  Interpolation  II 

Thomas  Bloom*  and  Jean-Paul  Calvi 

University  of  Toronto 
Toronto,  Ontario,  Canada 
and 

Universite  Toulouse 
Toulouse,  France 

Let  E,  F,  G  be  three  compact  sets  in  C”.  We  say  that  (E,  R,  G)  holds  if  for  any 
choice  of  an  interpolating  array  in  F  and  of  an  analytic  function  /  on  G,  the  corresponding 
Kergin  interpolation  polynomials  of  /  exists  and  converge  to  /  on  E.  Given  G,  regular 
C—  convex  and  E  (resp )  we  construct  optimally  F  (resp  E)  so  that  {E,  F,  G)  holds. 
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On  L2  and  Loo- Approximations  Subject  to  Interpolatory  Constraints 

M.  A.  Bokhaxi 

Department  of  Mathematical  Sciences,  K.F.U.P.M. 

Dhahran,  Saudi  Arabia 

Given  n  distinct  points  x-i,X2,- ■  ■  ,Xn  G  [a, 6],  a  function  /  e  C"'[a,6]  and  an  integer 
m  >  (r  +  l)n  —  1,  we  consider  the  problem 

min{||/-  ^11  :  ^  ^  =  f^Hxi),i  =  l,2,...,nJ  =  0,l,...,r}  (*) 

where  1|  •  ||  denotes  either  L2-norm  or  Loo-norm.  Using  a  technique  of  J.  L.  Walsh,  the 
problem  (*)  is  modified  to  an  unconstrained  minimization  problem.  We  note  that  the 
solution  of  (*)  possesses  several  similarities  to  the  respective  classical  cases  of  L2  and  Loo 
approximation  problems. 


Series  Representations  of  Entire  Functions  by 
Laguerre  and  Hermite  Polynomials 

Lyubomir  Boyadjiev 

Graceland  College 
Lamoni,  Iowa 

Eigenfunction  expansions  are  sometimes  useful  in  studying  the  properties  of  the  func¬ 
tions  to  which  they  converge.  Here  we  consider  entire  functions  and  their  expansions  in 
terms  of  Laguerre  and  Hermite  series.  The  class  of  entire  functions  which  have  expansions 
in  terms  of  Hermite  series  is  described  using  the  generating  function  for  Hermite  polyno¬ 
mials  as  a  basic  tool.  The  discussion  about  Laguerre  representations  of  entire  functions 
involve  the  classical  Hankel’s  integral  transformation.  Another  approach  to  this  problem, 
based  on  equiconvergence  theorems  for  series  in  Laguerre  polynomials,  is  also  considered. 


EJquiconvergence  of  Some  Simultaneous  Hermite-Pad4  Interpolants 
Marcel  G.  de  Bruin*  and  A.  Sharma 

Delft  University  of  Technology 
Delft,  The  Netherlands 
and 

University  of  Alberta 
Edmonton,  Canada 

The  object  is  to  extend  some  recent  results  of  E.  B.  Saff,  A.  Sharma  and  R.  S.  Varga 
(R.A.LR.O.  Anal.  Numer.  15  (1981),  no.  4,  371-390)  and  M.  P.  Stojanova  (Constr. 
prox.  4  (1988),  435-445)  on  equiconvergence  of  Hermite-Pade  interpolation  of  meromorphic 
functions  to  the  case  of  simultaneous  interpolation. 
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Let  d,  ui,  1^2,  ■•,^d  be  natural  numbers  and  let  F  =  (Fi, . . . ,  F^)  be  d  functions  mer- 
morphic  in  the  discs  Dp^  (1  <  pi  <  P2  •  •  •  <  Pd>  Pj  real)  respectively,  given  by 


:=  W  •=  £ Oii,kZ^y  limsup  |ai,fc|  <  ^  (i  =  1, . . . ,  d), 

■0*U)  *:-»oo  Pi 


where 

Bi{z)  :=  —  Zij)  Ai,j  = 

j=i  j=i 

Here  it  is  assumed  that  the  poles  of  the  fimctions  are  in  disjoint  sets  (i.e.  Zij  ^  Zk^i 
when  i  ^  k  and  all  j,l).  Let  £  be  a  positive  integer  and  put  n  =  a  +  1,  where  a  = 
+  1^2  +  •  •  •  +  finally  r  >  1  is  a  natural  number.  The  Hermite-interpolant  to  the 
Taylor-sections  on  the  zeros  of  (z”  -  1)''  will  be  denoted  by  fi,t(z)  (note 

that  the  case  ^  =  oo  is  well  defined). 

The  interpolation  problem  with  multiple  nodes  is  now  formulated  for  1  <  ^  <  oo  as 

findd  rational  functions  Ui{z)/B^{z),  degUf^z)  <  nr—i/i  —  l  (1  <  i  <  d),  degB^{z)  < 
n—UQ  —  l,  that  interpolate  the  d  rationals  fi^i{z)/B^{z)  (1  <  i  <  d)  on  the  zeros  of  (z”  — 1)^; 
moreover  the  coefficient  ofz^~'^~^  in  B^  will  be  1;  B^{z)  =  53fc=o°  ^  7n-i/o-i  ~ 

The  main  result  is  as  follows. 

A.  For  n  sufficiently  large,  the  interpolation  problem  stated  above,  has  a  imique 
solution  that  moreover  satisfies 


n— i/Q  —  l  « 

=  Ck,  with  Ckz'"  =  JjFi(z);  1  <  £  <  oo. 

fc=0 


i=l 


B.  Let  W  be  a  compact  subset  of  |z|  <  r,  r  >  0,  that  omits  the  singularities  of  the 
functions  Fi  (1  <  i  <  d),  then 


limsup{max 

n— »oo 


ur^jz)  ufiz) 

B°°{z)  B^{z) 


}'/”<c(i,T), 


where  c(i,r)  =  Bf^~^^^'^^{T/piY  for  r  >  pi  and  c(i,r)  =  for  r  <  pi  with  R  = 

maxi<i<d  1/pi. 

C.  Specifically  we  have  for  |z|  <  piF“(^~^'*'^/'’^: 

B“(z)  B‘{z)  ' 

uniformly  and  geometrically  in  compact  subsets  of  "H  omitting  the  singularities. 
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A  New  Look  at  the  Fekete-Szego  Theorem 

Jean-Paul  Calvi 

Universite  Toulouse  III 
Toulouse,  France 

Let  E  be  a  compact  point  set  in  C,  we  assume  that  its  logarithmic  capacity  is  greater 
than  1  {cap{E)  >  1).  Then  a  classical  [1955]  theorem  of  Fekete  and  Szego  states  that  there 
are  only  finitely  many  monic  polynomials  p{X)  e  Z[i][X]  irreducible  over  Q(i)  and  having 
all  their  roots  in  E,  i.e.,  p“^(0)  C  E. 

Thus  if  K{E)  denotes  the  set  of  all  the  roots  of  these  polynomials,  we  have  cardK (E)  < 
oo.  However,  Fekete  and  Szego  said  nothing  more  on  the  size  of  K  (^E).  In  this  note,  we  are 
interested  in  finding  specific  estimates  for  card  {K{E))  under  some  geometric  hypothesis  on 
E.  To  do  so,  following  the  original  idea  of  Fekete,  we  have  to  find  polynomials  q  €  Z[i][X] 
with  “small”  degree  and  such  that  ||g||f;  >  1.  (The  original  method  does  not  seem  to 
provide  bounds  for  the  degree.) 

We  propose  a  solution  which  makes  use  of  classical  Carleman’s  inequalities  for  orthog¬ 
onal  polynomials  with  respect  to  the  plane  Lebesgue  measure  (on  a  compact  set). 

Here  p(X)  G  Z[i\[X]  means  that  all  the  coefficients  of  p  belong  to  Z[i]  =  {a  +  ib,  a,  be 
Z}. 


Approximation  by  Boolean  Sums  of  Linear  Picugov— LehnhofF  Operators: 
Teljakovskii’s  Type  Estimates  and  Constants  of  Approximation 

Jia-ding  Cao*  and  Heinz  H.  Gonska 
Fudan  University 

Shanghai,  People’s  Republic  of  China 

We  estimate  constant  and  find  best  asymptotic  constant  in  Teljakovskii  s  Theorem 
concerning  approximation  by  algebraic  polynomials.  We  research  approximation  by 
Boolean  sums  of  linear  Picugov-Lehnhoff  operators. 


The  Construction  of  Minimal  Projections  via  Dual  Balls 

B.  L.  Chalmers*  and  F.  T.  Metcalf 

University  of  California 
Riverside,  California 

A  minimal  norm  projection  P  =  Ui  i8>  Vi  from  a  Banach  space  X  onto  V  = 
[ui, . . . ,  Un]  is  obtained  as  follows: 

||ai,---,an||  :=  ||aiui -f  . . . -I- On^nllx 
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defines  a  norm  in  n-space  with  ball  B. 

ui, . . . ,  tin  which  make  P  minimal  are  determined  from  the  requirement  tha 

:=  ||feltil  +  ...  +  6n«n||x* 

defines  a  norm  in  n-space  with  ball  B*,  the  dual  ball  to  B. 


A  Wavelet-like  Unconditional  Basis 

Kuei-Fang  Chang 

Feng  Chia  University 
Tai-Chimg,  Taiwan,  ROC 

We  construct  a  bounded  and  Besselian  basis  in  the  space  for  1  <  p  <  2,  and 

we  construct  a  bounded  and  Hilbertian  basis  in  ^^(R)  for  2  <  p  <  oo.  These  bases 
depend  on  an  orthonormal  wavelet  in  L^(R).  Under  the  Wiener  condition,  a  wavelet  and 
its  bi-orthogonal  wavelet  have  the  same  rate  of  decay  at  infinity.  If  a  pre-wavelet  has  a 
bi-orthogonal  pair  of  Riesz  bases,  then  the  pre-wavelet  provides  an  unconditional  basis  for 
all  LP(R)(1  <  p  <  oo). 


Spline  Wavelets  with  Prescribed  Shape 
Debao  Chen 

Oklahoma  State  University 
Stillwater,  Oklahoma 

In  some  previous  work,  we  studied  the  general  structure  of  cardinal  spline  wavelets, 
and  drew  graphs  of  some  of  them,  to  show  the  variety  of  shapes  included  in  these  families. 
In  this  paper  we  prove  that  we  can  construct  spline  wavelets  with  virtually  any  prescribed 
shape.  As  a  consequence,  we  can  prove  that  for  any  function  /  satisfying  certain  mild 
conditions,  the  set  •  -k)  :  j,k  E  Z},  (which  may  not  be  a  wavelet  basis  or 

wavelet  fi:ame)  is  dense  in  L^(R). 


Fast  Algorithms  for  Matrix- Valued  Nevanlinna-Pick  Interpolation 

Guanrong  Chen 
University  of  Houston 
Houston,  Texas 

In  this  talk,  we  introduce  several  fast  algorithms  for  both  scalar-valued  and  matrbc- 
valued  Nevanlinna-Pick  interpolation.  Given  n  distinct  points  Zi  in  the  unit  disk  |zl  <  1 
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and  n  (m  by  m)  complex  matrices  Wk  satisfying  the  Pick  condition,  the  fast  Nevanlinna- 
Pick  interpolation  algorithms  developed  by  the  author  require  only  0(Tn^  log  n)  complex 
arithmetic  operations  to  evaluate  the  matrix-valued  interpolatory  rational  function  at  any 
particular  value  of  2,  in  contrast  to  the  naive  Nevanlinna  algorithms  which  require  O(m^n) 
arithmetic  operations,  where  they  both  use  O(m^n)  processors. 


Denseness  of  Radial-Basis  Functions  in  L^(R^) 
and  Its  Applications  in  Neural  Networks 

Tianping  Chen*  and  Hong  Chen 
Fudan  University  &  University  of  Notre  Dame 

In  this  paper,  we  discuss  problems  of  approximation  to  functions  in  and  oi>- 

erators  from  to  by  Radial-Basis  Functions.  The  results  obtained  in  this 

paper  solve  the  problem  of  capability  of  RBF  neural  networks,  a  basic  problem  in  Neural 
networks. 

Suppose  g  :  R^.  Many  papers  concerning  approximation  by  functions 

N 

i=l 


have  been  published.  It  is  natural  to  raise  the  following  question:  What  are  the  necessary 
and  sufficient  conditions  imposed  on  the  function  g  such  that  the  family  of  functions 

Uti  CidMX  -  Xi\\Rn)  is  dense  in  L^R^)? 

The  purpose  of  this  paper  is  to  answer  this  question  and  construct  neural  network 
model,  which  can  be  used  to  approximate  nonlinear  continuous  operators  from  L^{R^)  to 
L^{R^).  The  main  results  are  as  follows. 


Theorem  1.  Suppose  g  ■.  R-^  ^  i^^(l  +  (|xl)^^  €  then  the  set  of  linear 

combinations  Ci5(Ai||A'  —  is  dense  in  L^{R^),  where  Ai  >  0, Cj  €  R^,Xi  E 

R^,i  = 


Theorem  2.  Suppose  V  is  a  compact  set  in  L^(R^),g  :  iZ'''  — ►  H^,p(||X||nn)  e  L?'(RP’), 
then  for  any  e  >  0,  there  are  positive  integer  N,  real  numbers  A*  >  0,  Xi  €  iT*,  i  =  1, . . . ,  iV, 
which  are  independent  of  f  E  V,  and  constants  Ci{f),i  =  1, . . . ,  iV,  depending  on  f  eV, 
such  that 


N 

fiX) 

i=l 


<  € 

L2(R") 


(1) 


holds  for  all  f  eV.  Moreover,  all  Ci{f)  are  linear  continuous  functionals  defined  on  V. 

Theorem  3.  Suppose  that  F  E  L'^{R^),  then  the  family  of  all  linear  combinations  X,i=i 
CiF{XiPiiX)  +  Zi)  is  dense  in  L‘^{R^),  where  Ci  E  R\Xi  E  R\Zi  E  R^,Pi  are  rotations  in 
Rr,k  =  \,...,N. 
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Remark  1.  The  assumption  made  in  Theorem  1  is  equivalent  to  ^(||X||;jn)  g 
which  is  also  necessary  for  the  validity  of  Theorem  1. 


Approximation  Using  Positive  Definite  Functions 

E.  W.  Cheney 

University  of  Texas  at  Austin 
Austin,  Texas 

Positive  definite  functions  (not  necessarily  radial)  can  be  used  to  advantage  in  approx¬ 
imating  real-valued  functions  on  IR”  or  on  S^.  In  both  cases,  we  usually  require  strictly 
positive  definite  functions,  and  thus  a  characterization  of  these  functions  is  needed  at  the 
start.  Not  only  can  these  functions  be  used  for  scattered  data  interpolation,  but  they 
form  “fundamental”  families  in  in  C{S"^),  and  in  C{S°°).  Recent  work  in  this 

area  will  be  summarized  in  the  talk.  In  particular,  contributions  by  Xingping  Sun,  Valdir 
Menegatto,  and  Kuei-Fang  Chang  will  be  discussed. 


Continuous  Wavelet  Transform  on  Tangent  Bundles  of  Spheres 

Stephan  Dahlke*  and  Peter  Maass 

RWTH  Aachen 
Aachen,  Germany 

In  this  paper,  the  continuous  wavelet  transforms  on  T^(R)  and  L^(R?)  are  generalized 
to  TS^  and  TS^.  We  consider  specific  groups  that  admit  square-integrable  representations 
in  L^{TS^)  and  L^(TS^),  respectively.  Under  the  aetion  of  these  groups,  points  in  the  fibers 
are  scaled  and  translated  whereas  points  in  the  underlying  manifolds  are  only  translated. 
Modified  versions  of  these  groups  can  be  used  to  define  the  generalized  wavelet  transform. 
It  turns  out  that  this  transform  is  a  multiple  of  an  isometry. 


Multiscale  Analysis,  Approximation,  and  Interpolation  Spaces 

Wolfgang  Dahmen 
RWTH 

Aachen,  Germany 

Striking  advancements  of  mathematical  modelling  and  numerical  simulation  during 
the  past  ten  or  fifteen  years  have  been  to  a  great  extent  due  to  a  significant  change  in 
the  traditional  profile  of  numerical  analysis.  Whereas  classical  numerical  linear  algebra 
schemes  tended  to  isolate  the  final  problem  formulation  (often  as  a  prohibitively  large  linear 
system  of  equations)  firom  its  origin,  a  breakthrough  in  the  complexity  of  solution  processes 
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could  be  achieved  by  exploiting  possibly  much  information  about  the  structure  of  the 
underlying  infinite  dimensional  continuous  problem.  The  objective  of  this  talk  is  to  indicate 
some  contributions  of  approximation  theory  in  this  context.  This  can  be  exemplified  for 
multiscale  basis  transformations  arising  in  connection  with  data  compression  or  multilevel 
solvers  for  differential  or  integral  equations.  The  stability  of  such  transformations  in  a 
general  Hilbert  space  context  can  be  reinterpreted  via  the  concepts  of  biorthogonality  and 
Riesz-bases  in  terms  of  certain  norm  equivalences.  Tools  from  interpolation  theory,  based 
on  classical  ingredients  like  Bernstein  and  Jackson  inequalities,  can  then  be  employed  to 
analyse  the  validity  of  such  norm  equivalences.  We  conclude  with  briefiy  indicating  some 
of  their  consequences  and  applications. 


Approximation  Theory  for  Erdos  Weights 

Steven  B.  Damelin 

University  of  the  Witwatersrand 
Republic  of  South  Africa 

In  recent  years  the  subject  of  weighted  approximation  on  the  real  line  has  received 
considerable  attention  with  the  development  of  a  general  theory  of  orthogonal  polynomials 
for  these  weights.  In  this  talk,  we  survey  the  topic  of  weighted  approximation  for  Erdos 
weights  with  particular  emphasis  on  recent  work  by  D.  S.  Lubinsky  and  the  author. 

Roughly,  an  Erdos  weight  is  the  form 

W{x)  =  exp(— Q(x)) 


where 


Q  :  R  — ^  R 

is  of  faster  than  polynomial  growth  at  infinity.  The  archetypal  example  being 

Q{x)  =  Qk,a{x)  =  expfc  (|x|“) ,  q;  >  1,  A:  >  1, 


where 

expk  =  exp  (exp(exp(. . .))) 

is  the  iterated  exponential.  The  results  to  be  presented  include  a  selection  from  sharp 
necessary  and  sufficient  conditions  for  Lagrange  interpolation,  sharp  necessary  and  suffi¬ 
cient  conditions  for  orthogonal  expansions  and  Jackson— Bernstein  approximation  theorems. 
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Better  Dual  Functions  for  Gabor  Lattices 

Ingrid  Daubechies 

Princeton  University 
Princeton,  New  Jersey 

(Abstract  not  available) 


A  New  Method  to  Determine  Regularity  of  Refinable  Functions 

Ingrid  Daubechies 

Princeton  University 
Princeton,  New  Jersey 

A  function  /  is  refinable  if  it  can  be  written  as  a  linear  combination  of  its  translated 
dilates  f{2x  -  n).  Such  functions  come  up  in  computer  aided  design,  in  subband  filtering, 
and  in  the  construction  of  wavelet  algorithms.  Several  methods  have  been  developed  to 
determine  the  regularity  of  /  firom  the  coefficients  in  the  refinement  equation.  Most  of 
these  methods  are  practical  only  if  the  linear  combination  in  the  refinement  equation  is 
finite.  In  joint  work  with  Albert  Cohen,  we  developed  a  new  method,  borrowed  from 
dynamical  systems  theory,  to  determine  the  regularity  of  refinable  functions;  this  new 
method  works  as  well  for  infinitely  many  as  finitely  many  non-zero  coefficients,  and  it 
extends  to  multidimensional  non-separable  refinement  equations. 


Characterization  of  the  Best  Uniform  Approximations  of  Periodic 
Functions  by  Convex  Classes  Defined  by  Strictly  CVD  Kernels 

Oleg  Davydov 

Dnepropetrovsk  State  University 
Dnepropetrovsk,  Ukraine 

We  say  that  a  kernel  K{x,y)  e  C{T^)  is  strictly  CVD,  if  the  number  of  zeros  of 

/(a:)=  [  K{x,y)h{y)dy 
Jo 

in  a  period  is  less  than  or  equal  to  the  number  of  cyclic  sign  changes  of  h,  for  any  h  € 
L~(0, 27r). 

According  to  S.  Karlin  {Total  Positivity  1,  p.  259),  we  have 

2fc-l 

Sfcdet  f  K{di,y)dy  >0, 

i.7=l 
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t?i  <  •  •  •  <  i92fc-i  <  1?1  +  27r,  /i  <  •  •  ■  <  hk-i  <  h  +  Stt,  n(Ij)  >  0,  j  =  1, ...  ,2k  -  1, 
where  Ij  is  an  interval,  is  +1  or  —1,  A:  =  1, 2,  — 

Set 

M  =  l^fix) :  f{x)  ^  J  K{x,y)h{y)dy,\h{y)\<l  a.  e.,j/ €  [0,27r)| . 

We  present  the  following  characterization  that  is  the  periodic  analog  of  an  result  of 
A.  Pinkus  (Journal  of  Approximation  Theory  (1981),  33  (2)). 

Theorem.  Assume  that  K(x,y)  is  strictly  CVD  and  g  €  C'(T’)  \  A4.  The  best  approxi¬ 
mation  „ 

r(x)=  /  K(x,y)h*(y)dy 
Jo 

to  g  from  A4  is  uniquely  characterized  as  follows. 

There  exists  an  n  =  1, 2, . . .  and  knots  <  *  *  •  <  C2n  <  +  27r  and  t9i  <  •  •  •  <  i92n  < 

+  27r,  such  that 

h*(y)  =  (-1)^+^  a.  e.,  0  <  y  <  6+1,  j  =  1, 2, . . . ,  2n; 
gi'dl)  —  /*(l9l)  =  —(5(1^2)  —  f*  (“d 2))  9i‘'^2n-l)  —  f  (l?2n-l) 

=  -(9(d'2J  -  r (^2n))  =  -(gid'L)  -  n^'D) 


(  I  rij+i  \ 

I  sign  det  jy  K(di,y)dy  1  \\g  —  f*\\c(T) 


for  some  t?2n,  satisfying 


l?2n-l  <  d2n  <  d2n  <  <  l9l+27r. 


Schatten— von  Neumann  Classes  Versus  Besov  Spaces 
of  Integral  Operator’s  Kernels 

Lubomir  T.  Dechevsky 

Sofia  Technical  University 
Sofia,  Bulgaria 

Let  S'p,  0  <  p  <  00,  be  the  Schatten-von  Neumann  ideal  of  compact  operators  bounded 
from  Ho  =  L2(R^)  to  Hi  =  L2(R”)  with  quasinorm  ||T|5pl|  =  (E^^i  SkiTY)"^^,  where 
Sk(T)  are  the  s-numbers  of  T. 

Let  G  R,  be  the  homogeneous  Besov  space  with  equal  metric  indices. 


28 


Our  main  result  states  that  the  factor-space  of  Sp  modulo  composition  with  unitaxies 
on  Hq  and  Hi  is  isomorphic  to  B  =  In  more  detail,  for  every 

T  e  Sp  there  exists  an  integral  operator  Ti: 

I'lfix)  =  f  K{x,y)f{y)dy,  x  e  R"  , 

7R"* 


such  that 

1)  Ti  =  UqTUi  ,  where  Uj  is  a  imitary  on  HjJ  =  0, 1,  and 

2)  ||r|5p|l  =  \\K\B\\. 

Conversely,  every  integral  operator  T  with  kernel  K  G  J5  is  in  Sp,  and  there  exist 
Cl,  C2  :  0  <  Cl  <  C2  <  oo,  depending  on  m,  n  and  p  only,  such  that 


ci\\T\Sp\\  <  \\K\B\\  <  C2\\T\Sp\\  . 


This  is  a  generalization  to  all  p  :  0  <  p  <  oo  of  the  well-known  theorem  about  Hilbert- 
Schmidt  operators  (p  =  2) — in  which  case  ci  =  C2  =  1  and  the  aforesaid  isomorphism 
is  an  isometry.  Best-approximation  properties  of  s-numbers  imply  that  the  present  result 
provides  direct  and  converse  theorems  relating  best  approximation  by  finite-rank  operators 
with  regularity  of  the  operator’s  kernel  in  terms  of  Besov  quasinorms. 

The  limiting  case  p  =  oo  is  of  extreme  importance  in  itself,  since  it  provides  an  explicit 
expression  for  the  uniform  operator  norm  of  a  compact  operator  between  Hq  and  Hi.  In 
particular,  the  norm  of  an  integral  operator  is  equivalent  to  the  norm  of  its  kernel  in 


Generailization  to  the  case  when  Hj,j  =  0, 1,  are  arbitrary  separable  Hilbert  spaces 
is  an  easy  application  of  the  standard  canonic- isometry  argument.  Moreover,  the  most 
general  case  when  Hq  and/or  Hi  is  non-separable  is  easily  reduced  to  the  separable  case 
because  the  image  of  any  compact  operator  is  separable. 

There  is  an  immediate  application  of  the  main  result  in  obtaining  apriori  estimates  for 
the  Fredholm  resolvent  T{\I  —  T)~^  and  for  the  solutions  of  Fredholm  integral  equations 
of  the  second  kind  in  terms  of  the  kernel’s  regularity.  Namely,  the  well-known  Carleman’s 
inequality  for  p  =  2  involving  the  T2-norm  of  the  kernel  can  now  be  extended  for  all 
p  :  0  <  p  <  oo. 

We  note  that  in  the  present  context  the  general  technique  of  proving  a  pair  of  adjusted 
Jackson-  and  Bernstein-type  inequalities  is  insufficient,  because  the  spaces  Xq  and  Xi  in 
which  approximation  is  being  carried  out  are  different  from  the  spaces  Vb  which 

appear  in  the  respective  iC-functional.  The  difficult  part  of  the  proof  is  to  show  that  Xj 
is  isomorphic  to  Yj,j  —  0, 1,  and  to  this  end  some  facts  related  to  unconditional  bases  of 
orthogonal  wavelets  are  crucial. 
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Some  Error  Estimates  on  the  Classes  of  Functions  for 
Interpolation  by  Smooth  Chebyshev  Generalised  Polynomials 

Vassili  Demidovitch 

Moscow  State  University 
Moscow,  Russia 

Let  {^(•)}o®  ^  be  a  ECT-system  defining  the  space  of  smooth  Chebyshev 

generalised  polynomials  (STG-polynomials)  7r[a,  6]  =  span{ei(')}o® • 

It  is  well-known  that  for  any  x(-)  G  [a,  6]  (m  >  0),  determined  by  the  table  of 

their  values  in  the  fixed  collection  of  knots 


An  =  e  la,b]  :  (i'  ^  i")  (ti>  Un)  (n  >  0), 

there  exists  a  unique  interpolation  STG-polynomial 

^„(-)  €  ^n[a,6]  span{ei(-)}o  C  ^[a,6], 

satisfying  the  following  conditions  =  x{ti){i  =  0, n).  To  estimate  the  error  Tn  = 

x{t)  —  we  introduce  the  following  notions: 

V{-;  the  STG-polynomials  (of  Newton’s  type)  with  respect  to  the  corre¬ 

sponding  collections  of  points 

{(^<fc)fcSo  ^  An  :  {k'  ^  k")  =>  {ik>  ^  4»)}  ; 

I?n,[a:(')] — generalized  derivative  of  degree  m  for  x{-)  (“by  virtue  of  the  ECT-system”); 

|a:(-)  G  :  mpc  |£)„,[a;(t)]|  <  Mm|. 

Then  we  can  prove  the  following  result. 

Proposition.  For  any  x(-)  G  a,  6](0  <  m  <  n-l-1)  the  structural  error  estimate 

holds  in  the  form  of  the  inequality 

(1)  |r„(i)|  <  (a<t<b) 

where  =^mint|D„^[^(^)]|  >  0,  and  fm(-)  is  independent  of  x(-)  nonnegative  function 
with  fm{ti)  =0  {i  =  0,n). 

In  particular,  one  can  assmne  that 

/n+l(t)  =  \V{t]to,...,tn)\ 
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and 


(2) 


fn{t)  “  2 


n''(‘  I  ^0>  •  •  •  >  ii—lt  ^<+1  >  •  •  •  >  tn) 


t=0 


1/n+l 


Remark.  It  follows  from  (2)  that  for  the  classical  case 


and 

/„(()  =  2 1{« 

Since  in  this  case  ‘dm  —  rn\,  then  (1)  implies  (when  m  =  n  +  1  and  m  =  n)  the  well-known 
error  estimates  on  the  classes  of  functions  for  interpolation  by  classical  polynomials. 


Nonlinear  Wavelet  Approximation 
R.  A.  DeVore 

University  of  South  Carolina 
Columbia,  South  Carolina 

We  shall  discuss  some  recent  aspects  of  nonlinear  wavelet  approximation  including 
best  basis  selection  and  adaptive  piursuit.  There  is  little  known  in  terms  of  error  estimates 
for  these  highly  nonlinear  methods  of  approximation.  We  shall  discuss  what  is  known  and 
future  directions. 


Best  Approximation  auid  RT-Functionals 
Z.  Ditzian 

University  of  Alberta 
Edmonton,  Alberta,  Canada 

The  best  approximation  from  a  subspace  of  eigenfunctions  of  a  given  differential  op¬ 
erator  is  related  to  the  RT-functionals  given  by  that  operator.  An  approximation  process  is 
introduced  and  its  rate  of  convergence  is  shown  to  be  equivalent  to  a  given  iC-functional 
via  a  strong  converse  inequality.  Application  to  simultaneous  approximation  is  given. 
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Best  Approximation  on  Some  Function  Spaces 

Xiaolei  Dong 

Harbin  Normal  University 
Harbin  City,  PRC 


No  abstract  received. 


Orthogonal  Polynomials  and  the  Construction  of  Scaling  Functions 
with  arbitrary  approximation  order 

G.  C.  Donovan*,  J.  S.  Geronimo,  and  D.  P.  Hardin 

Georgia  Institute  of  Technology 
Atlanta,  Georgia 

We  use  the  theory  of  orthogonal  polynomials  to  construct  continuous,  compactly  sup¬ 
ported,  orthogonal,  spline  scaling  functions  with  arbitrary  order  of  approximation. 


Non-diagonal  Quadratic  Hermite— Pad4  Approximation 
to  the  Exponential  Function 

Kathy  Driver 

University  of  the  Witwatersrand 
Johannesburg,  South  Africa 

Approximation  to  the  exponential  function  of  the  form 

£mn.(z)  ~  P.(z)e-^-  +  (3„{z)e-'  +  R.(z)  =  o  , 

where  Pn,  Qm  and  Rg  are  polynomials  of  degree  at  most  n,  m  and  s  respectively  and  Pn  has 
leading  coefficient  1  is  considered.  Exphcit  formulas  and  exact  asymptotics  are  obtained 
for  certain  ray  sequences  of  {Pn}>  {Qm}  and  {Rg}-  It  is  also  shown  that  these 

Hermite-Pade  Type  I  polynomials  can  be  used  to  asymptotically  minimize  expressions  of 
the  above  form  on  the  unit  disk. 
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On  Interlacing  Properties  of  the  Roots  of 
Orthogonal  and  Euler-Frobenius  Polynomials 

Francois  Dubeau*  and  Jean  Savoie 

Universite  de  Sherbrooke 
Sherbrooke,  Quebec,  Canada 

and 

College  Militaire  Royal  de  Saint-Jean 
Richelain,  Quebec,  Canada 

Ultraspherical  polynomials  and  Euler-Frobenius  polynomials  are  examples  of  sequences 
of  polynomials  generated  recursively  by  the  scheme:  Po{x)  =  {£  a,  nonnegative  integer) 

and  Cn+iPn+i  =  PrAPn)  where  Tr(P)(x)  =  —2rxP{x)  -F-  (1  —  x^)DP{x).  We  investigate 
interlacing  properties  of  the  roots  of  polynomials  Pn  and  Pn+i  or  P„  and  Pti+2-  We  ob¬ 
tain  similar  results  for  Hermite-like  polynomials  obtained  by  the  scheme  Hq{x)  =  x^  and 
Cn^\Hn+i  =  T{Hn)  where  T{H){x)  =  —2xH(x)  +  DH{x).  Finally,  an  increasing  property 
of  the  main  root  of  Euler-FVobenius  polynomials  is  obtained  for  the  most  important  cases. 


Nonstationary  Subdivision  Schemes 
and  Multiresolution  Analysis 

Nira  Dyn 

Tel  Aviv  University 
Tel-Aviv,  Israel 


Nonstationary  subdivision  schemes  consist  of  recursive  refinements  of  an  initial  sparse 
sequence  with  the  use  of  masks  that  may  vary  from  one  scale  to  the  next  finer  one.  The  case 
of  masks  with  supports  that  grow  linearly  from  one  scale  to  the  next  is  analysed.  General 
suflficient  conditions  for  the  convergence  of  such  schemes  to  C°°  compactly  supported  basis 
functions  are  derived.  Under  additional  mild  conditions,  these  limit  functions  allow  to 
define  a  multiresolution  analysis  that  has  the  property  of  spectral  approximation.  Finally, 
these  general  results  are  used  to  construct  C°°  compactly  supported  cardinal  interpolants 
and  also  C°°  compactly  supported  orthonormal  wavelet  bases  that  constitute  Riesz  bases 
for  Sobolev  spaces  of  any  order. 

The  talk  presents  joint  work  with  Albert  Cohen. 
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Asymptotic  Approximation  by  Cardinal  Series 

Ricardo  Estrada 

Texas  A&M  University 
College  Station,  Texas 

Using  the  theory  of  distributional  asymptotic  expansions,  we  obtain  asymptotic  ap¬ 
proximations  of  functions  by  series  of  cardinal  type, 

n=oo 

f{ne)g{x/e  -n),  as  e  0, 

n=“C50 

as  well  as  by  series  of  quasi-interpolants,  for  a  wide  class  of  functions  /  and  g. 


Hermite  Interpolation  with  Radial  Basis  Functions  on  the  Sphere 

Gregory  E.  Fasshauer 

Vanderbilt  University 
Nashville,  Tennessee 

We  show  how  conditionally  negative  definite  functions  on  the  sphere  can  be  used  for 
Hermite  interpolation.  These  results  are  obtained  by  combining  ideas  of  [Menegatto  ’94] 
and  [Narcowich  ’94],  as  well  as  classical  concepts  of  Schoenberg.  We  also  conclude  that 
the  basis  functions  we  use  are  natural  in  the  sense  that  Hermite  interpolation  can  be 
viewed  as  the  limit  of  an  associated  Lagrange  interpolation  problem,  thus  justifying  some 
observations  made  by  [Franke,  Hagen,  and  Nielson  ’93]. 


On  An  Algebraic  Modulus  of  Smoothness 
and  Best  Algebraic  Approximation 


Michael  Felten 

Universitat  Dortmund 
Dortmund,  Germany 


We  consider  the  space  of  continuous  fimctions  on  [—1,1]  as  well  as  the  space 
with  weight  ^  and  introduce  a  modulus  of  smoothness  that  is  based  on  the  algebraic 

addition  a  ©  6  :=  a\/l  —  ^  +  \/l  —  a^b  defined  on  [—1, 1]. 

It  will  be  shown  that  the  difference  operator  of  higher  order  behaves  completely  differ¬ 
ent  near  the  endpoints  ±1  than  the  corresponding  classical  difference  operator  for  periodic 
functions.  This  is  due  to  the  fact  that  ([-1, 1],  ©)  is  not  a  group  because  the  associative 
law  is  missing  around  the  boundaries. 

Steklov  TTipflns  and  an  equivalent  RT-functional  are  given.  Moreover  the  equivalence 
with  the  Butzer-Stens  modulus,  introduced  by  the  Chebychev  translation  will  be  seen. 
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Classification  of  All  Function  Systems  in  Lp 
and  the  Place  of  Wavelets  and  Prewavelets  in  this  Classification 

Vadim  Iv.  Filippov  ' 

Saratov  State  University 
Saratov,  Russia 

The  interest  to  the  wavelets  is  the  natural  development  of  some  part  of  function 
theory,  which  has  connection  with  concrete  function  systems  (trigonometric,  Haar,  etc.). 
It  concerns  applied  using  of  these  systems.  Since  there  exists  the  infinite  number  of  fxmction 
systems,  the  question  on  optimization  of  some  function  systems  (maybe  classical)  appears 
for  some  problems.  To  answer  on  these  questions  we  must  somehow  classify  all  function 
systems  (maybe  the  same  as  living  world  of  Earth  ). 

That  will  malfp  time  quick  for  the  solution  of  the  questions  of  optimizing.  The  prop¬ 
erties  of  function  systems  depend  on  function  spaces  (for  example  there  is  no  basis  in 
LP(0, 1),0  <  p  <  1).  Therefore,  it  is  natural  to  consider  the  function  systems  in  the 
spaces:  a)LP(0, 1),  0  <  p  <  1;  6)LP(0,  1),  1  <  p  <  oo,  (if  we  limit  ourselves  by  the  scale  of 
LP,p  >  0).  But  it  is  more  natural  to  consider  the  systems  in  all  subspaces  of  the  space  5, 
almost  everywhere  finite  measurable  functions  (for  example  the  scale  of  the  Orlicz  spaces) 
The  results  given  below  are  more  general  than  for  wavelets  in  the  paper  [1]  .  It  is  the  first 
step  to  the  classification  of  all  function  systems. 

Definition.  A  system  C  D’,0  <  p  <  oo  is  called  a  representation  system  in 

the  space  IP  if  for  any  f  G  there  exists  a  series  Ckfk  such  that 

n 

lim  |l/-^Cfc/fcl|Lp  =0. 
fc=l 

Let  suppipn  =  {a; :  ^n{x)  ^  0};  Xn  =  inf{x  €  [0, 1]  :  Vc  >  0,  mes{{x,x  +  e)  f]  supp(pn}  ^ 
0};  pn  =  sup{y  G  [0, 1] :  Ve  >  0,  mes{(y  -  e,  y)  H  suppcpn}  7^  0}; 

Denote  d((p„)  =  Vn-Xn,  Qn  =  [xn^Vn]-  We  consider  the  systems  C  L^{0,  l),p  > 
1,  satisfies  the  conditions: 

(1)  d{(pn)  Oj 

(2)  Ji^n{x)dx^0;  ^ 

(3)  VV  G  N  mes{U^jv  Qk}  =  1  let  <7n  =  infx^R\\(l  -  A(pn(x))xo„(x)||p(3^)'>.  Prom 
Lemma  1  [1]  it  follows  that,  if  the  (2)  are  satisfied,  then  an  <  1- 

Theorem  1.  Let  for  a  subsystem  {v’nfc(2^)}  of  system  {<Pti}  the  conditions  (1),  (2)  and 
supn^o^nk  <  1  ore  satisfied.  Then  this  subsystem  >  where  TNT  G  N  is  any  number, 

is  a  representation  system  in  L^(0,  l),q  <  p,  if  and  orily  if  the  condition  (3)  is  fulfilled. 

There  axe  examples  for  non-improving  of  assumptions  (l)-(3)  in  Theorem  1.  If  the 
assumption  (3)  is  not  satisfied,  then  the  system  {(pn}  is  not  a  full  system.  If  the  assumption 
(2)  or  (1)  is  not  satisfied  then  there  exist  the  examples,  when  the  system  {<pn}  is  not  a  full 
one. 


35 


Remark  If  the  ^stem  {(pn)  consists  of  subsystems,  which  satisfy  the  assumptions  in  the 
II  level  or  in  level  and  etc.,  then  we  must  consider  all  these  subsystems. 

[1]  V.  Filippov,  P.  Oswald.  Representation  in  U*  by  series  of  translates  and  dilates  of  one 
function,  J.  Appr.  Theor.  (to  appear). 


On  Best  Linear  Approximation  in  /i(n) 

Martina  Finzel 

Universitat  Erlangen 
Niirnberg,  Germany 

Consider  the  space  li(n);  i.e.,  R”  endowed  with  |  •  |i,  and  let  \Pu  :  R”  — ►  2^  be  the 
metric  projection  of  l\  (n)  onto  an  r-dimensional  linear  subspace  U.  The  subspace  and  its 
metric  complement  :=  decompose  (in  general  not  uniquely)  R”;  i.e., 

R”  =  1/(0)  +  U. 

The  author  explicitly  describes  the  metric  complement:  Let  R  :=  U-^  n  6“  (0)  denote  the 
intersection  of  the  orthogonal  complement  of  U  with  the  closed  dual  unit  ball  of  loo{n);  R  is 
a  polyhedron  of  dimension  n—r.  The  faces  of  the  polyhedron  are  considered  to  be  relatively 
open.  Each  face  of  R  is  uniquely  determined  by  a  point  of  its  relative  interior,  e.g.,  by  its 
center  of  gravity  w,  and  denoted  by  face^w;  obviously,  rel  bdR  = 

Let  Su  :=  D  5'i(0)  denote  the  intersection  of  with  the  /i(n)-imit  sphere  the 
so-called  Schattengrenze  of  17;  in  other  words,  =  cone  {Su]  0).  Each  face  foxxyj  of  R 
determines  a  face  S^,  of  51(0)  belonging  to  the  Schattengrenze;  in  particular, 

Su  =  (J  Sru- 

w^extR 

Introducing  =  cone  {S^,  0)  +  U,  which  is  for  u;  G  ext  R  an  open  cone  with  vertex  in  0, 
leads  to  a  decomposition  of  R”  into  pairwise  disjoint  cones;  i.e., 

R”\17=  ij  and  R”  =  |J 

face^CR  weext  R 


In  particular,  if  the  metric  projection  onto  U  is  unique,  the  decomposition  is  unique: 

Rn  ^  ^(0)  ^  ^  (jong  (5-^.  0)  0  (7  for  all  face,^  C  R. 

In  this  case  the  metric  projection  is  linear  on  each  of  the  cones  w  €  ext  R.  For 
an  arbitrary  U  these  considerations  lead  to  a  selection  of  the  metric  projection  which  is 
characterised  by  the  following  properties: 

•  the  error  6  =  x  —  u*  has  minimal  support  for  the  selected  u*  €  iPu{x)] 
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•  the  selection  is  piecewise  linear,  and  there  is  a  finite  number  of  cones  of  linearity 
decomposing  R”;  in  particular,  the  selection  is  globally  Lipschitz  continuous; 

•  w.r.t.  these  properties  the  selection  is  in  a  way  “closest”  to  the  strict  approximation 
in  looiV')' 


Orthogonal  Polynomials  and  Solving  Linear  Systems  of  Equations 

Bernd  Fischer 

Medical  University  of  Liibeck 
Liibeck,  Germany 

In  this  note  we  discuss  the  connection  between  orthogonal  polynomials  and  the  solu¬ 
tion  of  linear  systems  of  equations  by  polynomial  iteration  methods. 

To  this  end,  let  (•,  •)  denote  an  (appropriate)  inner  product  for  the  space  n„  of  poly¬ 
nomials  of  degree  n.  The  design  of  an  effective  polynomial  iteration  method  leads  to 
approximation  problems  of  the  form 

(pT»,Pn)  =  min{(p,p)  :  p  G  n„,  p(0)  =  1}, 


or  to  problems  like 

(pn,Pn)  =  min{(p,p) :  p  €  Iln,  p(0)  =  1,  p'(0)  =  0}, 

respectively.  We  note  that  the  latter  problem  is  in  particular  interesting  in  connection 
with  singular  systems. 

The  solution  of  the  first  problem  is  given  by  the  Kernel  polynomials  with  respect  to 
(•,•).  Also,  for  the  second  problem  the  unique  solution  is  explicitly  known  and  may  be 
called  Hermite  Kernel  polynomial.  We  will  characterize  the  Hermite  Kernel  polynomial 
and  we  will  devise  (stable)  recurrence  relations  for  the  computation  of  these  polynomials. 
Furthermore,  we  will  point  out  relationships  between  the  two  mentioned  polynomial  classes 
and  we  will  outline  the  connection  to  the  actual  implementation  of  polynomial  iteration 
methods. 


Plane  Preserving  Volume  Deformations  Using 
Scattered  Landmark  Points 

T.  Foley 

Arizona  State  University 
Tempe,  Arizona 

We  address  the  problem  of  deforming  one  three  dimensional  volume  to  another  volume. 
For  our  applications,  landmark  points  in  the  domain  volume  are  arbitrarily  located  and  it 
is  desired  to  construct  a  deformation  to  another  volume  so  that  the  scattered  landmark 
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points  get  mapped  to  another  set  of  corresponding  points  in  the  range  volume.  Applica¬ 
tions  involve  modeling  solids,  volumetric  morphing  and  mapping  one  brain  volume  to  a 
“standard”  brain  volume.  For  the  brain  deformation,  we  modify  the  interpolant  so  that 
the  separating  plane  in  the  first  brain  is  also  mapped  to  the  separating  plane  of  the  second 
brain. 


Computation  of  Pad4  Approximations  of  Transfer  Functions 

via  the  Lanczos  Process 

Roland  W.  Freund 

AT&T  Bell  Laboratories 
Murray  Hill,  New  Jersey 

The  behavior  of  time-invariant  single-input  single-output  linear  dynamical  systems 
is  characterized  by  their  frequency-domain  transfer  functions.  Although  these  functions 
are  rational,  their  numerator  and  denominator  polynomials  are  of  very  high  degree  if 
the  underlying  system  is  large,  and  then  it  becomes  prohibitive  to  compute  the  transfer 
function.  Instead,  for  large  systems,  one  approximates  the  transfer  function  by  a  rational 
function  of  lower  order.  A  common  choice  is  to  employ  Fade  approximations,  which  are 
computed  in  the  usual  fashion  by  first  generating  the  moments,  followed  by  a  solution  of 
a  Ha/nkpl  system.  Unfortunately,  this  approach  is  numerically  unstable  due  to  the  typical 
extreme  ill-conditioning  of  the  Hankel  system,  and  it  can  only  be  used  to  compute  Fade 
approximants  of  very  moderate  order. 

In  this  talk,  we  present  an  alternative  algorithm  for  the  stable  and  efficient  computa¬ 
tion  of  Fade  approximations  of  transfer  functions.  Exploiting  the  connection  between  Fade 
approximation  and  the  Lanczos  process  for  the  successive  tridiagonaUzation  of  large  matri¬ 
ces,  our  algorithm  generates  the  Fade  approximants  directly  from  the  tridiagonal  Lanczos 
matrices,  and  it  thus  avoids  the  ill-conditioned  moment  computations.  We  present  nu¬ 
merical  results  for  a  variety  of  linear  systems  arising  in  the  simulation  of  lairge  electronic 
circmts. 

For  the  general  case  of  multi-input  multi-output  linear  dynamical  systems,  the  associ¬ 
ated  transfer  functions  are  matrix-valued  rational  fimctions.  We  also  describe  an  algorithm 
for  the  computation  of  matrix  Fade  approximants  to  such  matrix-valued  transfer  functions. 
The  algorithm  is  based  on  a  new  Lanczos  process  for  multiple  starting  vectors. 

Fart  of  this  work  is  joint  with  Feter  Feldmann,  AT&T  Bell  Laboratories. 
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Relations  Between  Best  Approximations 
Michael  Ganzburg 

Courant  Institute  of  Mathematical  Sciences 
New  York,  New  York 

Let  E{f,  B,  F)  be  the  best  approximation  of  an  element  /  by  elements  from  B  in  the 
metric  of  a  normed  space  F.  Relations  between  best  approximations  E{f,  B,  F)  by  differ¬ 
ent  subspaces  B  play  an  important  role  in  Approximation  Theory.  Recently  the  author  has 
developed  a  general  approaeh  for  obtaining  such  relations,  including  the  cases  of  approx¬ 
imations  by  algebraic  and  trigonometric  polynomials  and  entire  functions  of  exponential 
type  (Bernstein,  Ganzburg),  of  approximations  by  trigonometric  polynomials  and  spline 
functions  (Velikin),  of  rational  and  spline  approximations  (Petrushev),  and  of  weighted 
polynomial  approximations  with  different  weights. 

These  results  have  different  applications  in  Approximation  Theory  such  as  multi¬ 
dimensional  Jackson  theorems,  multidimensional  Bernstein-,  Nikolskii-,  and  Remez-t)rpe 
inequalities,  sharp  constants  of  multidimensional  approximations  and  others. 

In  this  talk  we  shall  present  a  survey  of  these  problems  and  give  new  relations  between 
best  approximations.  In  particular,  some  results  are  extended  to  the  case  of  quasinormed 
spaces  F.  For  instance,  for  0  <  p  <  1, 

lim  E{f,Vn,m,Lj,(nV*))  =  Eif,Bv,Lp{R^)),  (1) 

n“+oo 

where  R*”  is  the  m-dimensional  Euclidean  space,  P  is  a  centrally  symmetric  convex  com- 
pactum  in  R^,  V*  is  the  polar  of  V,  By  is  the  class  of  entire  functions  of  exponential 
type  whose  spectrum  is  contained  in  V,  'Pn.m  is  the  class  of  polynomials  in  m  variables  of 

degree  n,  Lp{Q)  is  the  space  with  the  quasinorm  ||/llLp(n)  =  (/  0  <  p  <  1. 

n 

For  1  <  p  <  oo  the  relations  like  (1)  have  been  obtained  by  Bernstein,  Raicin  and 
Ganzburg.  This  relation  is  based  on  Remez-  and  Nikolskii-type  inequalities  for  multidi¬ 
mensional  polynomials  in  symmetric  quasinormed  spaces. 


Approximation  by  Separable  Functions 
Manfred  v.  Golitschek 

Institut  fiir  Angewandte  Mathematik  und  Statistik 
Universitat  Wurzburg 
Wurzburg,  Germany 


The  approximation  by  separable  functions  is  one  of  the  topics  in  Approximation  The¬ 
ory  and  Numerical  Analysis  where  Ward  Cheney  has  contributed  many  important  results 
during  the  last  20  years. 
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Separable  Functions  are  defined  as  follows.  Let  iS  C  R  ,  T  C  R  and  D  C  S  x  T  he 
comp>act  sets.  Two  finite-dimensional  subspaces  U  C  £oo{S)  and  V  C  £oo(T)  of  bounded 
functions,  with  basis  (uj)^  or  respectively,  generate  the  linear  space  IFoo(L))  of 

separable  functions 


m  n 

(1)  w{s,t)  =  ^^yk{s)vk(t),  {s,i)  €  D, 

j=l  k=l 

where  Xj  e  £oo{T)  and  y*,  G  ^oo(‘S')  are  arbitrary  bounded  functions.  In  the  applications, 
S,  T  and  D  are  finite  sets  and  U  and  V  are  polynomial  or  spline  spaces. 

Two  problems  will  be  discussed: 

(i)  Given  /  G  £oo{D)  or  /  G  C{D),  does  there  exist  a  best  approximation  of  / 

firom  Woo{D)  in  the  vmiform  norm  ||  •  Hd  on  £>  ? 

(ii)  If  A  is  a  (finite)  subset  of  D  and  is  a  best  approximation  of  /  from  VFoo(A) 

in  the  uniform  norm  ||  •  ||a  on  A,  then  compare  the  distances 

||/-ty^||A,  II/-w^IId,  U-w^Wd- 


Approximations  by  Periodic  Biorthogonal  Wavelets:  Numerical  Algorithms 

Sonia  M.  Gomes* ,  Elsa  Cortina,  and  Irene  Moroz 

IMECC  -  Universidade  de  Campinas 
Campinas  SP,  Brasil 

In  this  paper  we  develop  fast  numerical  algorithms  for  representations  of  1-periodic 
functions  in  biorthogonal  multiresolution  analyses.  The  algorithms  include  decomposition, 
reconstruction,  interpolation,  point  evaluation  and  representation  of  differential  operators 
having  constant  coefficients.  Since  all  these  operations  are  of  convolution  type,  it  is  natmal 
to  make  use  of  fast  Fourier  transform  technique.  We  show  that  in  the  Fourier  domain 
the  decomposition  and  reconstruction  algorithms  have  a  matrix  representation  in  terms  of 
permutations  and  of  2  x  2  block  diagonal  matrices.  For  representation  of  derivatives,  linear 
systems  with  circulant  matrices  must  be  solved  which  can  be  done  directly  using  discrete 
Fourier  transform.  One  important  fact  about  biorthogonal  multiresolution  analyses  is 
that  they  are  compatible  with  derivatives.  We  show  that  this  property  can  be  used  for 
easy  computations  of  stiffness  matrix  elements.  It  also  helps  in  the  study  of  the  Gibbs 
phenomenon  in  approximations  of  functions  having  jump  discontinuities  in  some  of  their 
derivatives.  We  present  some  illustrative  examples  using  biorthogonal  spline  wavelets. 
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Parallel  Wavelet-Galerkin  Methods  by  Means  of  Adapted  Wavelet  Packet  Bases 

C.  Guerrini*  and  E.  Piraccini 

University  of  Bologna 
Bobgna,  Italy 

Wavelet  approximation  methods  for  the  numerical  solution  of  differential  and  pseudod¬ 
ifferential  equations  have  been  recently  investigated  as  a  powerful  tool  for  preconditioning 
and  compressing  the  stiffness  matrices  relative  to  orthogonal  and  biorthogonal  wavelet 
bases  [1],  [2],  [3];  moreover  stability  and  convergence  results  have  been  established  for  a 
very  general  setting  [4]. 

In  recent  studies,  new  compactly  supported  biorthogonal  wavelet  bases  have  been 
constructed,  which  are  adapted  to  some  simple  differential  operators  [5];  on  these  bases 
the  stiffness  matrices  have  a  simple  structure  and  bounded  condition  number.  Looking 
for  an  efficient  parallel  algorithm  for  the  numerical  solution  of  the  Galerkin  approach, 
in  this  paper  we  propose  the  use  of  compactly  supported  wavelet  packet  bases  adapted 
to  quite  general,  constant  coefficient,  elliptic  differential  operators.  We  show  how  these 
bases  can  be  evaluated  by  means  of  simple  convolution  operations  of  finite  support  vectors 
and  we  study,  specially  firom  a  numerical  point  of  view,  the  consequences  of  using  these 
bases  for  the  parallel  solution  of  simple,  constant  coefficient  model  problems.  Performance 
and  efficiency  of  the  parallel  algorithm  are  estimated  from  a  first  experimentation  on  a 
distributed  memory  hypercube  multiprocessor. 

[1]  W.  Dahmen,  S.  Proessdorf,  and  R.  Schneider,  Wavelet  approximation  methods  for 
pseudodifferentiaJ  equations  II:  Matrix  compression  and  fast  solution.  Preprint  No.82, 
RWTH  Aachen,  1992. 

[2]  S.  Jaffard,  Wavelet  method  for  fast  resolution  of  elliptic  problems,  SIAM  J.  Nu- 
mer.Anal.  29  (1992),  965-986. 

[3]  G.  Beylkin,  On  wavelet-based  algorithm  for  solving  differential  equations,  preprint 
1993. 

[4]  W.  Dahmen,  S.  Proessdorf  and  R.  Schneider,  Wavelet  approximation  methods  for 
pseudodifferential  equations  I:  Stability  and  convergence,  lAAS-Preprint  No.  7,  Berlin 
1992. 

[5]  S.  Dahlke,  I.  Weinreich,  Wavelet-Galerkin  methods:  An  adapted  biorthogonal  wavelet 
basis.  Const.  Approx.  9  (1993),  237-262. 


Error  Bounds  for  Linear  Interpolation  on  a  Triangle 

David  Handscomb 

Oxford  University  Computing  Laboratory 
Oxford,  England 

We  attempt  to  determine  best  possible  bounds  for  the  errors  in  function  value  and 
derivative  when  a  function  is  approximated  by  linear  interpolation  between  values  at  the 
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vertices  of  a  triangle  of  known  shape  and  size.  We  consider  both  Ls  and  Loo  bounds,  in 
terms  of  L2  and  Loo  measures  of  smoothness  of  the  function. 


Orthogonal  Multiwavelet  Constructions  in  2D 

G.  C.  Donovan,  J.  S.  Geronimo,  and  D.  P.  Hardin* 

Vanderbilt  University 
Nashville,  Tennessee 

In  the  usual  construction  of  a  wavelet  using  a  multiresolution  analysis  (MRA)  (V^)  is 
a  single  function  (p  called  the  scaling  function  whose  set  of  integer  translates  forms  a  Riesz 
basis  for  Vq.  In  contrast  multiwavelets  are  constructed  from  an  MRA  that  is  generated 
by  a  finite  set  of  scaling  functions  If  there  is  some  set  of  compactly  supported 

scaling  functions  whose  integer  translates  form  an  orthogonal  basis  for  Vq  then  (V^)  is 
called  an  orthogonal  MRA. 

The  authors  and  co-worker  Massopust  constructed  the  first  nontrivial  continuous,  com¬ 
pactly  supported,  orthogonal  multiwavelets  using  univariate  fractal  interpolation  functions. 
These  multiwavelets  have  several  properties  not  possible  in  the  single  scaling  function  case. 
For  instance,  they  are  symmetric,  interpolatory  and  remain  orthogonal  when  restricted  to 
integer  intervals. 

In  recent  work  by  the  authors,  a  general  scheme  for  constructing  multivariate  wavelets 
was  developed.  This  scheme  is  based  on  enlarging  a  given  multiresolution  analysis  (usually 
a  classical  spline  space)  by  adding  certain  functions  to  obtain  an  orthogonal  MRA. 

This  work  is  extended  to  the  multivariate  setting,  in  particular  we  add  continuous 
“fractal  interpolation  surfaces”  that  are  supported  on  triangles  to  the  MRA  generated  by 
the  Cq  piecewise  linear  hat  function  on  a  given  lattice  to  get  an  orthogonal  MRA.  As  in 
the  imivariate  case  the  resulting  scaling  functions  are  symmetric  (with  respect  to  rota¬ 
tions  leaving  the  lattice  invariant),  interpolatory  and  remain  orthogonal  when  restricted 
to  triangles  in  the  associated  triangulation. 


Peano  Kernel  for  Harmonicity  Differences  of  Order  p 

Werner  Haufimann 

Gerhard-Mercator-University  Duisburg 
Duisburg,  Germany 

We  consider  linear  functionals  L  which  vanish  on  poly  harmonic  functions  of  a  fixed 
order  p  >  1  in  a  bounded  domain  I>  C  R*^,  i.e.,  on  the  kernel  of  the  poly  harmonic  operator 
AP,  where  A  denotes  the  Laplacian  differential  operator.  For  such  functionals  we  prove  a 
theorem  of  Peano  type  which  states  that 

L(/)=  f  V(x)A’’f(x)dj:  (/eC^”), 

JD 
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where  V  is  the  corresponding  Peano  kernel. 

This  formula  shows  a  full  analogy  with  the  classical  one-dimensional  case  in  which 


Vo{t)g^^\t)dt 


(9  €  cn, 


and  where  L  is  a  functional  vanishing  on  the  polynomials  of  degree  p  —  1,  i.e.,  on  the  kernel 
of  the  operator  —  . 

It  is  a  well-known  resvdt  that  the  Peano  kernel  Vo  of  the  one-dimensional  finite  differ¬ 
ence  of  order  p  is  a  univariate  B-spline. 

We  introduce  the  concept  of  the  harmonicity  difference  of  order  p  and  study  properties 
of  the  Peano  kernel  associated  with  it.  The  harmonicity  difference  of  order  p  naturally 
arises  from  the  Pizzetti  mean  value  formula  for  polyharmonic  functions  of  order  p. 

The  harmonicity  difference  of  order  p  enjoys  many  of  the  properties  of  the  one¬ 
dimensional  finite  difference.  In  particular,  we  prove  that  its  Peano  kernel  is  (i)  com¬ 
pactly  supported,  (ii)  radially  symmetric  with  high  order  smoothness,  and  (iii)  satisfying 
an  extremal  problem  which  resembles  Holladay’s  theorem  for  one-dimensional  splines. 

Our  paper  is  joint  work  with  O.  Koimchev  (Sofia). 


The  Orthogonality  of  Weighted  Faber  Polynomials 

Matthew  He 

Nova  Southeastern  University 
Ft.  Lauderdale,  Florida 

Weighted  Faber  polynomials  {Fn{z',  g;  E)}^^q  associated  with  a  domain  E  and  a 
weight  function  g  play  a  very  important  part  in  the  study  of  the  asymptotic  properties 
of  orthogonal  polynomials  in  the  complex  domain.  In  the  present  paper,  we  study  the 
orthogonality  of  these  polynomials,  mainly  in  dependence  on  the  weight  function  g  and 
the  structure  of  the  domain  E.  We  introduce  a  class  of  m-fold  orthogonal  polynomials  and 
investigate  their  properties.  Our  results  include  certain  classical  orthogonal  polynomials 
and  ordinary  Faber  polynomials  as  special  cases. 


Spline  Interpolation  and  Its  Wavelet  Analysis 
Tian  Xiao  He 

Illinois  Wesleyan  University 
Bloomington,  Illinois 

In  this  paper,  we  will  discuss  the  spline  interpolants  at  arbitrary  knots  with  the  highest 
possible  approximation  order  and  the  corresponding  wavelet  functions.  Considering  the 
spline  interpolants  as  window  functions,  we  will  also  discuss  the  corresponding  compact 
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window  Fourier  transform.  In  addition,  an  approach  to  the  corresponding  orthogonal 
scaling  function  will  be  also  considered. 


Existence  and  Accuracy  for  Matrix  Refinement  Equations 

Christopher  Heil 

Georgia  Tech 
Atlanta,  Georgia 

Refinement  equations  are  functional  equations  of  the  form  f{x)  =  '^k=o  /(2^  ~  ^)- 

They  play  a  key  role  in  wavelet  theory  and  in  subdivision  schemes  in  approximation  the¬ 
ory.  The  coefficients  Ck  are  usually  real  or  complex  numbers.  Recent  “multiwavelet” 
constructions  require  matrix  Ck  and  vector- valued  /  =  (/i,  •  •  • ,  /r)-  We  consider  the  exis¬ 
tence,  uniqueness,  and  order  of  approximation  of  compactly  supported  solutions  to  matrix 
refinement  equations.  Several  classes  of  solutions  axe  identified,  depending  on  whether 
the  infinite  matrix  product  (FI^i  ■^(2~'^7))  converges,  possibly  in  a  weak  or  constrained 
sense,  where  M('y)  =  5  This  behavior  is  determined  by  the  eigenvalues  of 

M(0).  We  also  show  how  to  calculate  the  order  of  approximation  p,  i.e.,  the  highest  degree 
polynomial  l,x,. . . ,  x^~^  that  can  be  written  exactly  as  a  linear  combination  of  the  integer 
translates  of  /i, . . . ,  /r. 


Stability  and  Linear  Independence 
of  the  Shifts  of  A  Multivariate  Refinable  Function 

Tom  Hogan 

University  of  Wisconsin — Madison 
Madison,  Wisconsin 

Due  to  their  so-called  time-frequency  localization  properties,  wavelets  have  become  a 
powerful  tool  in  signal  analysis  and  image  processing.  Typical  constructions  of  wavelets 
depend  on  the  orthogonality  or  stability  of  the  shifts  of  an  underlying  refinable  function.  In 
this  paper,  these  properties  are  characterized  in  terms  of  the  refinement  mask  for  a  large 
class  of  multivariate  compactly  supported  refinable  distributions.  Linear  independence 
(which  is  related  to  stability)  of  the  shifts  is  also  treated.  Among  the  functions  to  which 
these  results  apply  are  tensor  products  and  box  splines. 
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Optimal  Triangulations  Using  Edge  Swappings 

Dong  Hong 

The  University  of  Texas  at  Austin 
Austin,  Texas 

This  is  joint  work  with  C.  K.  Chm.  The  main  result  of  this  paper  is  an  efficient  method 
for  triangulating  any  finitely  many  scattered  sample  sites  such  that  these  sample  sites  are 
the  only  vertices  of  the  triangulation  and  that  for  any  discrete  data  given  at  these  sample 
sites,  there  is  a  piecewise  quartic  polynomial  on  this  triangulation  that  interpolates  the 
given  data,  with  the  optimal  order  of  approximation. 


Recovering  Band  Limited,  Weakly  Stationary  Process 
from  a  Set  of  Its  Irregularly  Spaced  Samples 

Christian  Houdre 

Georgia  Tech 
Atlanta,  Georgia 

The  problem  of  recovering,  say,  a  band  limited  weakly  stationary  process  from  a  set 
of  its  irregularly  spaced  samples  is  studied.  For  rather  general  sampling  sequences  some 
sufficient  conditions  ensuring  mean  square  or  pathwise  reconstruction  are  obtained.  For 
the  cases  of  regular  samples  with  either  finitely  many  missing  ones  and/or  finitely  many 
irregular  ones  a  necessary  and  sufficient  condition  is  presented.  Some  elements  of  the  proofs 
involve  classical  resxilts  on  non-harmonic  Fourier  series  as  well  as  more  recent  results  on 
frames. 


Convex  Spline  Approximation  in  Lp  Space 

R.  A.  DeVore,  Y.  K.  Hu*,  and  D.  Leviatan 

Georgia  Southern  University 
Statesboro,  Georgia 

Let  /  G  Lp[— 1,  1],  0  <  p  <  00,  be  convex.  Then  for  any  given  partition  T„  :=  {—1  := 
xo  <  xi  <  ...  <  Xn  :=  1}  of  the  interval  [—1,  1],  there  is  a  convex  continuous  piecewise 
quadratic  function  on  Ik  with  approximation  error  us,  either  measured  on  [—1,  1]  or 
on  local  subintervals.  This  function  can  be  smoothed  into  a  or  (cubic)  spline  with 
the  same  degree  of  approximation. 

REMARK:  This  is  our  first  step  of  constructing  convex  polynomial  approximation  in 
Lp  with  degree  of  approximation  0/3 . 
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Generalizations  of  Dykstra’s  Algorithm 
Hein  Hundal 

Pennsylvania  State  University 
State  College,  Pennsylvania 


This  talk  covers  two  generalizations  of  Dykstra’s  Algorithm.  Dykstra’s  algorithm  is 
a  method  for  computing  the  best  approximations  to  a  set  C  which  is  the  intersection  of  a 
finite  number  of  closed  convex  sets  Cj.  The  best  approximation  is  found  by  an  alternating 
projections  scheme.  This  scheme  is  useful  when  the  projection  onto  the  individual  sets  is 
easily  computed,  but  the  projection  onto  C  is  unknown  or  difficult  to  compute.  The  first 
generalization  of  this  algorithm  presented  allows  the  programmer  to  change  the  order  in 
which  the  projections  are  performed.  These  changes  can  be  used  to  accelerate  convergence. 
The  second  generalization  allows  C  to  be  an  intersection  of  an  infinite  sequence  of  closed 
convex  sets.  Since  all  convex  sets  in  a  separable  Hilbert  space  can  be  represented  as 
the  intersection  of  a  countable  number  of  half  spaces,  this  generalization  could  be  useful. 
Examples  of  both  generalizations  and  some  simple  applications  will  be  given. 


Scattered  Data  Reconstruction  and 
Radial  Basis  Function  Interpolation  Methods 

Armin  Iske 

Institut  fur  Numerische  imd  Angewandte  Mathematik 
Universitat  Gottingen 
Gottingen,  Germany 

We  assert  to  reconstruct  a  finite  number  of  scattered  data  Tif,...,TNf  that  are 
generated  by  actions  of  linear  functionals  Ti, . . .  ,Tn  on  a  function  /  :  R*^  — +  R  from  a 
function  space  to  be  specified.  A  possible  specialization  of  our  method  could  be  multivariate 
interpolation.  For  this  case,  radial  basis  functions  provide  a  prominent  tool. 

Here,  we  generalize  radial  basis  function  interpolation  methods  and  assert  to  find 
a  reconstruction  function  Sf  with  respect  to  the  reconstruction  conditions  Tjf  —  TjSf, 
{1  <  j  <  N).  Due  to  the  reconstruction  scheme,  s/  is  assumed  to  be  a  linear  combination 
of  translates  of  a  prescribed  conditionally  positive  definite  function  ip  exceeded  by  a  poly¬ 
nomial  from  which  denotes  the  linear  space  of  all  d  variate  polynomials  of  order  m  or 
less. 

Certainly,  the  reconstruction  functions  as  described  above  form  a  linear  space 
depending  on  ip.  We  assert  to  topologize  by  a  seminorm  |  •  |^  which  is  intimately 
related  to  the  conditional  positive  definiteness  of  ip.  Then,  (F^/,,  |  •  |^)  is  referred  as  native 
function  space  of  the  reconstruction  method. 

This  talk  provides  a  characterization  of  the  native  function  spaces  and  relates  different 
terms  of  conditional  positive  definiteness. 
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Wavelet  Smoothing  Techniques  and  Interactive  Surface  Modeling 

B.  Jawerth 

University  of  South  Carolina 
Columbia,  South  Carolina 

Given  a  set  of  measured  data,  we  discuss  several  different  wavelet  based  techniques  to 
quickly  fit  a  smooth  surface  to  the  points.  Further,  given  a  surface,  we  present  a  method, 
based  on  a  version  of  Nash-Moser  techniques  combined  with  wavelets,  to  edit  the  surface 
to  obtain  a  new  one  with  desired,  user-defined  properties. 


A  Generalization  of  de  Boor’s  Stability  Result 
and  Symmetric  Preconditioning 

Kurt  Jetter 
University  of  Duisburg 
Duisburg,  Germany 

We  provide  a  symmetric  preconditioning  method  based  on  weighted  divided  differences 
which  can  be  applied  in  order  to  solve  ill-conditioned  scattered-data  interpolation  problems 
in  a  stable  way.  Concerning  the  theoretical  background,  an  a  priori  unboxmded  operator 
on  £2(2)  is  preconditioned  so  as  to  get  a  bounded  and  coercive  operator.  The  method  has 
another  interpretation  in  t^ms  of  constructing  certain  Riesz  bases  of  appropriate  closed 
subspaces  of  L2{R^)>  In  extending  Mallat’s  multiresolution  analysis  to  the  scattered  data 
case,  we  construct  nested  sequences  of  spaces  giving  rise  to  orthogonal  decompositions 
of  functions  in  L2{R^)\  in  this  way  the  idea  of  wavelet  decompositions  is  (theoretically) 
carried  over  to  scattered-data  methods. 

This  is  joint  work  with  J.  Stockier,  University  of  Duisburg. 


Refinable  Shift-invariant  Spaces: 

EVom  Splines  to  Wavelets 

Rong-Qing  Jia 
University  of  Alberta 
Edmonton,  Canada 

In  this  talk  we  shall  survey  some  recent  results  concerning  refinable  shift-invariant 
spaces.  A  linear  space  S  of  (measurable)  functions  on  R*  is  said  to  be  shift-invariant  if  it 
is  invariant  under  integer  translation,  that  is, 

feS=^f{^-j)eS  VjGZ^ 

The  space  S  is  said  to  be  refinable  if 

/  e  S  =»  n-/2)  e  S. 
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Refinable  shift-invariant  spaces  form  a  natural  framework  for  the  study  of  splines  and 
wavelets. 

In  this  talk  we  cover  five  topics.  First,  we  review  some  basic  properties  of  shift- 
invariant  spaces  such  as  stability,  linear  independence,  and  orthogonality,  and  discuss 
characterizations  of  shift-invariant  spaces.  Second,  we  give  a  unified  treatment  to  various 
approximation  schemes  such  as  orthogonal  projection,  bi-orthogonal  projection,  cardinal 
interpolation,  and  quasi-interpolation.  We  also  provide  an  insight  into  some  recent  re¬ 
sults  on  characterization  of  the  approximation  order  of  a  shift-invariant  space.  Third,  we 
study  refinement  equations  and  subdivision  schemes.  We  are  especially  concerned  with 
the  existence  and  smoothness  of  the  solution  of  a  refinement  equation,  and  the  rate  of 
convergence  of  the  corresponding  subdivision  scheme.  Foxnth,  we  investigate  multiresolu¬ 
tion  and  construction  of  wavelets.  The  emphasis  is  placed  on  construction  of  multivariate 
wavelets  which  are  not  tensor  products  of  imivariate  wavelets.  Finally,  we  examine  results 
on  approximation  by  wavelets.  In  particular,  the  direct  theorem  (Jackson  estimates)  and 
the  inverse  theorem  (Bernstein  estimates)  are  discussed  in  the  setting  of  Besov  spaces. 


Some  Inequalities  for  Derivatives 

Theodore  Kilgore 

Auburn  University 
Auburn,  Alabama 

Some  inequalities  for  derivatives  will  be  discussed. 


Wavelet  Analysis  of  Speech  Signals 

Mei  Kobayashi*  and  Masaharu  Sakamoto 

IBM  Tokyo  Research  Laboratory 
Tokyo,  Japan 

This  paper  describes  experiments  in  which  speech  signals  were  analyzed  using  several 
types  of  wavelets:  Haar,  Daubechies,  spline,  Bagor,  and  chirp.  We  constructed  wavelet 
sonograms  of  time  versus  resolution  space  and  compared  them  with  those  from  a  Fourier 
method.  Since  wavelet  analysis  showed  promise  for  identifying  some  speech  events,  follow¬ 
up  experiments  were  conducted  to  assess  the  viability  of  wavelet-based  tools  in  speech 
recognition  systems.  We  will  present  our  findings  and  discuss  possible  directions  for  future 
research. 
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A  Variational  Approach  to  Subdivision 

Leif  Kobbelt 

CAGD-Group 
University  of  Karlsruhe 
Karlsruhe,  Germany 

Interpolatory  subdivision  is  a  technique  to  iteratively  produce  smooth  curves  or  siu-- 
faces.  This  is  done  by  inserting  new  points  into  a  given  control  polygon  or  control  net.  As 
more  and  more  points  are  inserted,  the  sequence  of  polygons  or  nets  converges  to  a  smooth 
limit. 

In  the  hteratiure,  subdivision  schemes  are  usually  defined  by  an  explicit  rule  how  to 
compute  the  new  points  in  every  step.  In  our  approach,  however,  we  implicitly  define  the 
new  points  by  choosing  an  energy  functional  which  is  to  be  minimized.  Although  this  in 
general  leads  to  global  dependencies  between  new  and  old  points,  the  resulting  schemes 
can  be  computed  with  linear  time  complexity  due  to  the  banded  structure  of  the  equation 
systems  that  have  to  be  solved. 

We  give  general  conditions  for  the  existence  of  solutions  and  present  some  examples 
for  such  schemes  which  produce  limit  ciuwes  of  arbitrary  smoothness.  There  are  several 
ways  to  modify  this  approach. 


Geometric  Hermite  Interpolation 

K.  Hollig  and  J.  Koch* 

Mathematisches  Institut  A 
Stuttgart,  Germany 

Since  parametric  representations  of  cmves  are  not  unique,  the  approximation  rates  by 
splines  can  be  significantly  improved.  This  surprising  fact  was  first  observed  in  [DeBoor- 
Hoellig-Sabin:1987]  where  a  6-th  order  accurate  cubic  interpolation  scheme  for  planar 
curves  has  been  constructed.  We  describe  a  natural  generalization  of  standard  Hermite 
interpolation  for  space  cmves.  In  addition  to  position  and  tangent  direction  at  the  end¬ 
points  a  third  point  within  the  parameter  interval  is  interpolated.  The  resulting  method  is 
easy  to  implement  and  achieves  the  optimal  approximation  order  5.  This  is  a  substantial 
improvement  over  the  order  4  of  standard  Taylor  approximations.  The  high  accuracy  of 
the  method  is  confirmed  by  examples  which  also  demonstrate  a  good  qualitative  behav¬ 
ior  of  the  approximations.  In  particular  the  curvature  plots  show  that  the  shape  of  the 
interpolated  curves  is  in  general  very  well  preserved. 

1.  C.  de  Boor,  K.  Hollig,  and  M.  Sabin,  High  accuracy  geometric  hermite  interpolation. 
Computer  Aided  Geometric  Design,  4  (692),  269-278,  April  1987. 
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On  the  Laguerre  Polynomials  of  Negative  Index 

Elizabeth  Kochneff 

Eastern  Washington  University 
Cheney,  Washington 

The  Laguerre  polynomials  of  index  a  >  —  1,  {L“(x)},  are  the  unique  polynomials 
orthogonal  on  (0,  cx>)  with  respect  to  the  weight  x°‘e~^  dx  which  satisfy  the  requirements 
that  each  L“  is  a  polynomial  of  exact  degree  n  and  the  coefficient  of  x”  is  (— l)"/n!. 
They  are  defined  for  a;  <  —1  using  their  explicit  representation  for  a  >  —1  and  analytic 
continuation.  We  prove  pointwise  convergence  for  expansions  of  certain  functions  in  terms 
of  the  Laguerre  polynomials  of  index  a  <  —  1.  The  idea  is  that  although  the  kernel 
polynoroial  cannot  be  integrated  against  the  weight  it  behaves  as  an  approximation 

to  the  identity  provided  we  apply  it  to  functions  which  are  well  behaved  near  zero  or  if  we 
subtract  from  the  product  of  the  kernel  polynomial  and  the  function  a  suitable  polynomial. 


Asymptotically  Best  L,-approximation  by  Splines  and 
Asymptotically  Best  Quadrature  Formulae 

Peter  Kohler 

Universitat  Hildesheim 
Hildesheim,  Germany 

By  an  appropriate  choice  of  knots  and  data  points  for  the  quasi-interpolants  of  de 
Boor  and  Fix,  asymptotically  best  spline  approximants  to  a  given  smooth  function  g  can 
be  obtained,  where  the  underlying  norm  is  the  L,-norm  with  1  <  g  <  oo.  This  method  can 
be  used  to  construct  Peano  kernels  for  asymptotically  best  quadrature  formulae.  Instead 
of  using  these  Peano  kernels,  the  quadrature  formulae  resulting  in  this  way  can  also  be 
obtained  simply  by  replacing  the  integrand  /  by  some  piecewice  polynomial  interpolant 
(this  works  for  integrals  with  arbitrary  integrable  weight  functions). 


On  A:— monotone  Polynomial  and  Spline  Approximation 
in  Lp,  0  <  p  <  00  Metric 

Kirill  Kopotun 

University  of  Alberta 
Edmonton,  Canada 

Let  n„  denote  the  set  of  all  algebraic  polynomials  of  degree  <  n,  «5 (r,  N)  be  the  space 
of  all  splines  of  order  r  with  knots  ^rn{f,  ^)p  denote  the  usual  integral  modulus  of 

smoothness,  and  let  A*'  be  the  set  of  all  functions  /  such  that  ^ 

for  all  X  and  h  such  that  [x,x+kh]  C  [0, 1].  In  particular,  A^  is  the  set  of  all  nondecreasing 
functions.  Also,  /  e  C^[0, 1]  Cl  A^  4=^  f^'^\x)  >0,xe  [0, 1]. 
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The  following  theorem,  which  is  a  generalization  (in  the  case  1  <  p  <  oo)  and  a 
complement  (0  <  p  <  1)  of  a  well  known  result  of  A.  S.  Shvedov,  is  discussed. 

Theorem  1.  Let  k  £  N  and  0  <  p  <  oo  be  fixed,  m  £  N\J  {0},  and  let  u  £  N  be  such 
that  max{A:  +  2  —  m,k}  <  u  <  k-h  p~^.  Then  for  any  n  £  N  and  A  >  0  there  exists  a 

function  f  £  C7®®[0, 1],  f^^\x)  >0,x£  [0, 1]  such  that  for  every  Pn  €  Iln,  Pn^\o)  >0  the 
follotving  inequality  holds 

Corollary,  (monotone  polynomial  and  spline  approximation).  The  estimates 

<  Ca-2(/M)p 


and 

4‘2(/)p  <  c«2(/M)p 

are  not  true  in  general  for  0  <  p  <  oo  and  f  £  1]  fl  A^,  where  En\f)p  := 

lly  -fnllp  and  Sr,n(^f')p  I—  iufag5(r,n)nA^  11.^  ®llp* 

It  is  interesting,  that  in  the  case  p  =  oo  these  estimates  are  true.  In  fact,  the  following 
is  valid: 

if  /  e  1]  n  A^,  then  <  Cn-^u}m(f,n-^)oo,m  >  1  (I.  A.  Shevchuk, 

and  X.  M.  Yu  and  Y.  P.  Ma),  and  i^r,n(/)oo  ^  C!n~^u}r-i{f',n~^)ooj‘>'  ^  2  (D.  Leviatan 
and  H.  N.  Mhaskar). 

Theorem  1  shows  that  the  estimates  in  the  following  theorem  are  exact  in  the  sense 
of  the  orders  of  moduli  of  smoothness. 

Theorem  A.  (monotone  spline  approximation).  Let  L^[0,  1],  1  <  p  <  oo  denote  the  space 
of  functions  which  are  j -fold  integrals  of  Lp[0,  1]  functions.  Then 

4^iif)p<Cn~^uj{f',n~\, 

if  f  E  I'pfO,  1]  n  A^  and  r  >2  (C.  K.  Chui,  P.  W.  Smith  and  J.  D.  Ward)  and 

C*(/)p  < 

if  f  £  Pp[0,  1]  n  A^  and  r  >3  (D.  Leviatan  and  H.  N.  Mhaskar). 


51 


Generalized  Exponential  Box  Splines 
Thomas  Kunkle 

College  of  Charleston 
Charleston,  South  Carolina 

Up  to  a  constant  factor,  the  s-variate  exponential  box  spline  of  Ron,  Dahmen,  and 
Micchelli  is  the  convolution  of  several  distributions  M  of  the  form 

{M,  <p)-  I  M(t  I  ro, . . . ,  Xn)<piKt)  dt 

Jxo 

where  ^  is  an  s-variate  test  function,  ac  is  a  vector  in  R'*,  and  M{t\xo, . . .  ,x„)  is  the 
univariate  5-spline  with  equidistant  knots  {xo,  • .  •  jXn}-  In  this  talk,  we  will  discuss  the 
compactly  supported  Generalized  Exponential  Box  Spline  (or  GEB-spline)  obtained  by 
letting  {xo, . . .  ,x„}  have  variable  stepsize. 

The  recurrence  relations  for  GEB-splines  closely  resemble  those  of  the  exponential  box 
spline  and  univariate  5-spline. 

Taking  for  a  basis  a  set  of  shifts  of  finitely  many  GEB-splines,  we  obtain  a  spline  space 
with  a  triangulation  of  periodic  stepsize  in  each  of  its  directions. 

Results  concerning  a  partition  of  unity,  the  exponential-polynomials  contained  in  this 
space,  and  the  linear  independence  of  these  shifts  all  generalize  corresponding  results  for 
the  exponential  box  spline. 


Bivariate  Splines  on  FVS-triangulation 
Ming-Jun  Lai 
University  of  Georgia 
Athens,  Georgia 

We  survey  the  study  of  bivariate  splines  on  a  FVS-triangulation.  FVS-triangulation  is 
a  triangulation  obtained  firom  a  quadrangulation  by  adding  two  diagonals  of  each  quadri¬ 
lateral.  Such  a  triangulation  was  first  considered  by  Praeijs  de  Veubeke  and  Sander  in 
1974.  They  studied  piecewise  cubic  splines  on  this  special  triangulation  to  construct 
finite  elements.  Their  study  was  generalized  by  Laghchim-lahlou  and  Sablonniere  in  1989. 
They  constructed  finite  elements  in  spline  spaces  of  smoothness  r  and  degree  3r  for 
odd  integer  r  and  in  spline  spaces  S^Sr+l  of  smoothness  r  and  degree  3r-l- 1  for  even  integer 
r  on  a  FVS-triangulation.  In  Lai  1993,  an  application  of  Si  to  CAGD  was  considered. 
Lai  showed  the  dimension  of  Si  on  such  a  triangulation  is  smaller  than  other  bivariate 
spline  spaces.  In  Lai  and  Schumaker  1994,  they  studied  the  construction  of  locally  sup¬ 
ported  basis  and  approximation  order  of  Sg  on  such  a  triangulation.  They  demonstrated 
that  the  dimension  of  this  spline  space  is  smaller  than  many  spline  spaces.  In  Lai  1994, 
Lai  studied  the  approximation  order  of  bivariate  spline  spaces  52^^)  on  a  FVS-triangulation 
for  d(r)  =  3r  if  r  is  odd  and  d(r)  =  3r  -  1  if  /  >  2  is  even.  He  showed  those  spline  spaces 
with  lower  degree  d(r)  <  3r  + 2  have  full  approximation  order. 
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Composite  Wavelet  Transforms 
Joe  Lakey 

Texas  A&M  University 
College  Station,  Texas 

Meaningful  unifications  of  the  theory  of  affine  firames  and  Gabor  frames  are  obtained 
by  means  of  group  extensions  that  contain  both  affine  and  Heisenberg  subgroups.  Exten¬ 
sions  of  the  continuous  wavelet  and  Gabor  transforms  £ire  discussed  and  two  methods  for 
designing  corresponding  firames  adaptable  to  analyzing  specific  signals  are  introduced. 


Monotone  Rational  Approximation 

E.  H.  Kaufman,  Jr.,  D.  J.  Leeming*,  and  G.  D.  Taylor 

University  of  Victoria 
Victoria,  B.C.,  Canada 

We  consider  the  problem  of  uniform  approximation  on  a  finite  set  of  real  numbers 
by  rational  functions  whose  first  derivative  is  required  to  be  nonnegative  at  each  of  the 
points.  We  discuss  computational  results  obtained  by  applying  the  CONMAX  nonlinear 
optimization  code  (developed  by  the  authors)  to  this  problem.  Some  theorems  concerning 
best  approximations  are  presented,  and  some  possible  extensions  of  the  work  are  mentioned. 


Multivariate  Polynomial  Interpolation: 

Explicitation  of  Lagrange  and  Hermite  Polynomials 

Alain  Le  Mehaute 

Telecom  Bretagne 
Brest,  Rrance 

We  present  some  explicit  algorithm  for  the  expUcitation  of  Multivariate  Hermite  or 
Lagrange  Polynomials  interpolating  data  in  or  R**.  Techniques  from  real  algebraic 
geometry  axe  well  suited  for  Lagrange  interpolation  problems  in  the  plane,  providing  algo¬ 
rithms  to  decide  whether  a  data  set  is  unisolvent  or  not,  and  also  a  recursive  algorithm  that 
enables  us  to  construct  the  associated  Lagrange  basis.  For  a  Newton  approach,  it  is  more 
relevant  to  use  a  multiresolution  algorithm  which  is  available  in  any  dimension.  Finally, 
we  consider  the  Abel  interpolation  problem  in  R*^,  i.e.  the  interpolation  problem  where  a 
point  Aa  is  associated  to  the  derivative  5“/  for  any  a  such  that  |a|  <  k.  Some  of  the  points 
a  can  coalesce,  and  the  most  obvious  important  case  of  Abel  interpolation  corresponds  to 
the  Taylor  interpolation  at  a  single  point.  The  importance  of  these  interpolation  problems 
comes  from  the  fact  that  an  interpolation  scheme  is  always  regular  (i.e.,  solvable  whatever 
the  location  of  the  data  points  is)  if  and  only  if  the  interpolation  problem  is  an  Abel’s  one. 
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We  provide  an  explicit  formula  for  the  solution  which  can  be  considered  as  a  Multipoint 
Taylor  Formula. 


Approximation  FVom  Shift-Invariant  Spaces  by  Certain  Integral  Operators 

Junjiang  Lei 

Southern  Illinois  University 
Carbondale,  Illinois 

We  consider  approximation  from  shift- invariant  spaces  by  using  an  integral  operator. 
We  characterize  the  approximation  order  provided  by  such  an  integral  operator.  Applica¬ 
tions  of  this  characterization  are  also  given. 


On  a  Density  Question  for  Neural  Networks 

Burkhard  Lenze 

Fachhochschule  Dortmund 
Dortmund,  Germany 

One  of  the  most  interesting  questions  connecting  neural  networks  and  approximation 
theory  concerns  the  principal  capability  of  three-layer  feedforward  networks  to  simulate 
a  given  function  up  to  any  desired  accuracy.  The  most  complete  and  satisfactory  result 
for  standard  McCulloch-Pitts  nets  was  recently  given  by  Leshno  et  al.  [2] .  In  this  talk  we 
consider  a  related  question  where  the  McCulloch-Pitts  neurons  are  substituted  by  so-called 
hyperbolic  sigma-pi  neurons.  We  will  prove  that  three-layer  feedforward  neural  networks 
with  hyperbolic  sigma-pi  neurons  in  the  hidden  layer  can  approximate  any  continuous  mul¬ 
tivariate  function  up  to  any  desired  accuracy,  more  precisely,  that  they  induce  a  translation 
and  dilation  invariant  subspace  of  CCR”)  in  the  topology  of  uniform  convergence  on  com¬ 
pact  subsets.  For  those  interested  in  the  general  neural  network  background  we  refer  to 
[1],  [2],  and  [3];  generalizations  of  our  results  by  means  of  completely  different  techniques 
may  be  found  in  [4]. 
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Generalised  sA:-spline  Interpolants  on  Compact  Abelian  Groups 

Jeremy  Levesley 

University  of  Leicester 
Leicester,  United  Kingdom 

Here  we  investigate  interpolation  and  approximation  of  continuous  functions  on  a 
compact  Abelian  group,  using  shifts  of  some  fixed  continuous  function.  We  also  consider 
some  extremal  properties  of  such  interpolants.  Finally,  for  a  wide  range  of  smooth  functions 
on  the  torus,  we  obtain  error  estimates  for  generalised  sfc-spline  interpolants. 


Convex  Polynomial  Approximation  in  the  Lp  Norm 

R.  A.  DeVore,  Y.  K.  Hu,  and  D.  Leviatan* 

Tel  Aviv  University 
Tel  Aviv,  Israel 

If  F  €  Lp[-l,l],  0  <  p  <  oo,  is  convex,  we  would  like  to  estimate  the  degree  of 
polynomial  approximation  to  /  under  the  constraint  that  the  polynomials  are  convex. 
Estimates  so  far  involved  the  second  modulus  of  smoothness  of  /  in  the  Lp-norm.  While  it 
was  known  that  estimates  involving  the  fourth  modulus  of  smoothness  of  /  axe  impossible 
with  an  absolute  constant,  it  was  an  open  question  whether  estimates  involving  the  third 
are  valid.  We  show  that  for  such  an  /  and  for  each  n  >  2,  there  exists  convex  polynomials 
Pn  €  n„  such  that 

where  C  =  C(p)  when  p  — »  0. 
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Distribution  of  Zeros  and  Poles  of  Rational  Functions 
that  are  Asymtopically  Optimal  for  the  Zolotarev  Problem 

Eli  Levin*  and  E.  B.  Saff 

The  Open  University  of  Israel 
Tel  Aviv,  Israel 
and 

The  University  of  Florida 
Tampa,  Florida 

Given  two  disjoint  compact  subsets  Ex,  E2  of  the  complex  plane,  the  generalized  third 
Zolotarev  problem  concerns  estimates  for  the  quantity 


Zmn  • —  i^f  ^mn(^) 


where  Rm«  denotes  the  class  of  all  rational  functions  with  numerator  degree  m  and  de¬ 
nominator  degree  n.  Recently  we  showed  that  for  any  *'ray  sequence  of  integers 
that  is  m/n  A,m  -h  n  oo^zUiT'^^  approaches  a  limit  L{\).  A  sequence  {r^„}  is 
called  asymptotically  optimal  if 


z„„  -  L(X). 


Under  some  geometric  assumptions  on  Ei,E2^  we  describe  the  asymptotic  behaviour  of 
zeros  and  poles  of  such  rmn’s.  As  a  consequence,  we  describe  (for  the  case  A  <  1)  the 
behaviour  of  rational  functions  that  are  asymptotically  optimal  for  the  fourth  Zolotarev 
problem  (approximation  of  a  signum-type  function). 


Multivariate  Interpolating  AVavelets  and  Duals 

Chim  Li 

Center  for  Approximation  Theory 
Texas  A&M  University 
College  Station,  Texas 

This  is  joint  work  with  C.  K.  Chui.  We  apply  our  general  frame  work  of  multivariate 
wavelets  with  duals  to  the  interpolatory  setting  and  construct  multivariate  interpolating 
wavelets  as  well  as  their  distributional  duals.  Examples  in  terms  of  multivariate  box  splines 
are  given. 
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Minimization  of  Multivariate  Quadratic  Splines 

W.  Li 

Old  Dominion  University 
Norfolk,  Virginia 

We  study  multivariate  quadratic  splines  that  have  simple  matrix  representations.  The 
problem  of  minimizing  such  a  spline  function  is  closely  related  to  quadratic  programming 
problems,  linear  programming  problems,  Huber’s  M-estimator,  and  best  -approximation 
problems.  There  are  many  different  issues  related  to  minimization  of  quadratic  splines, 
including  computational  complexity,  combinatorial  analysis,  algorithms,  error  bounds,  and 
applications.  We  will  give  a  brief  overview  of  this  interesting  topic. 


Laguerre  Series  Method  for  Finding  the  Best  Hankel  Approximation 

of  Infinite  Dimensional  Systems 

Xin  Li 

University  of  Nevada 
Las  Vegas,  Nevada 

Orthonormal  Laguerre  basis  for  L^(0,  oo)  is  applied  to  compute  the  best  Hankel 
approximation  of  certain  classes  of  continuous-time  infinite  dimensional  systems  in  time 
domain.  And  moreover,  the  modulus  of  continuity  of  best  Hankel  approximation  operator 
is  derived. 


Characterization  of  the  Order  of  Polynomial- Reproduction 
for  Multi-Scaling  Functions 

Jian-ao  Lian 

Prairie  View  A&M  University 
Prairie  View,  Texas 

If  the  integer  translates  of  a  compactly  supported  uni-scaling  function  y?  form  a  Riesz 
basis  for  a  subspace  Vq  of  L^(R),  it  is  well-known  that  ip  has  polynomial-reproduction  order 
m  if  and  only  if  the  two-scale  symbol  P  of  y?  is  divisible  by  the  factor  (1  -f  z)*”  but  not  by 
However,  this  result  does  not  extend  to  the  multi-scaling  function  setting.  The 
objective  of  this  paper  is  to  establish  certain  necessary  and  sufficient  conditions  in  terms 
of  the  eigenvalues  of  two  finite  matrices  for  which  a  multi-scaling  function  has  polynomial- 
reproduction  order  m.  Some  applications  will  be  given. 
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Hybrid  Cubic  Bezier  Triangle  Patches  on  Spheres  and  Sphere-like  Surfaces 

X.-L.  Liu*  and  Larry  L.  Schumaker 

Vanderbilt  University 
Nashville,  Tennessee 

We  discuss  a  new  method  for  fitting  scattered  data  on  the  sphere  based  on  the 
spherical  barycentric  coordinates  and  Bernstein-Bezier  polynomials  introduced  by  Alfeld, 
Neamtu  and  Schumaker  [ms,  1994].  The  method  extends  Foley  and  Opitz’s  hybrid  cubic 
patches  [Thomas  A.  Foley  and  Karsten  Opitz,  Mathematical  Methods  in  Computer  Aided 
Geometric  Design  n,  pp.  275-286]  to  spheres  and  general  sphere-like  surfaces.  The  hybrid 
form  thus  defined  gives  a  local  representation  of  cubic  triangular  patches,  and  can  be 
evaluated  with  the  deCasteljau  algorithm.  A  interpolant  to  scattered  data  with  cubic 
precision  is  obtained  which  is  significantly  smoother  than  Clough-Tocher  macro  elements 
with  linear  cross  boundary  derivatives. 


Bdzier  Wavelets 
Suresh  K.  Lodha 

University  of  California,  Santa  Cruz 
Santa  Cruz,  California 

Classical  approach  to  wavelets  based  on  Mallat’s  multiresolution  analysis  for  functions 
on  the  real  line  has  been  extended  to  functions  on  bounded  intervals.  These  wavelet-like 
bases  or  multiwavelet  bases  for  bounded  intervals  are  much  easier  to  construct  and  are  not 
unique.  A  unique  construction  can  be  achieved  by  imposing  additional  requirements.  For 
example,  additional  vanishing  moments  criteria  have  been  used  to  create  Legendre  wavelets 
that  have  been  useful  in  sparse  representation  of  integral  operators.  More  recently  other 
wavelets  have  been  created  that  do  not  necessarily  have  the  translation-invariant  or  scale- 
invariant  properties  and  have  found  useful  applications  in  automatic  level-of-detail  control 
for  editing  and  rendering  curves  and  surfaces.  Orthogonal  properties  of  wavelets  have 
also  been  given  up  to  construct  biorthogonal  wavelets  on  bounded  intervals,  where  the 
dual  scaling  functions  are  interpolating  and  have  found  applications  in  compressing  eleva¬ 
tion  contours.  Other  applications  of  wavelets  on  bounded  intervals  include  fast  methods 
for  solving  problems  in  global  illumination,  image  editing  and  compression,  and  surface 
reconstruction  firom  contours. 

There  is  a  need  to  construct  wavelets  on  boimded  intervals  whose  coefficients  reflect 
certain  geometric  properties  that  are  inherent  in  scaling  functions  or  are  suitable  to  the  a,p- 
plications  at  hand.  Different  basis  or  scaling  functions  however  possess  different  geometric 
properties.  For  example,  coefficients  in  the  Bezier  representation  of  curves  carry  geometric 
properties  that  allow  a  designer  to  control  the  shape  of  the  curve.  In  contrast,  coefficients 
in  the  Lagrange  representation  of  curves  carry  geometric  property  of  interpolation.  Certain 
desirable  geometric  properties  of  wavelets  from  the  point  of  view  of  geometry  applications 
are  symmetry,  smoothness,  and  small  support. 
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Here  we  investigate  the  construction  of  wavelets  for  curves  represented  in  the  Bezier 
bases.  The  de  Casteljau  algorithm  can  be  viewed  as  the  reconstruction  relations  for  the 
Bezier  scaling  fimctions.  In  this  work  we  study  the  decomposition  relations  for  the  Bezier 
basis  functions.  In  particular,  we  construct  Bezier  wavelets  whose  coefficients  carry  some 
precise  geometric  information.  We  then  discuss  an  application  of  this  representation  to 
smoothness  preserving  approximation  of  curves. 


Optimal  Polynomial  Orthogonal  Systems 

George  Lorentz 

University  of  Texas,  Austin 
Austin,  Texas 


(Abstract  not  available) 


Trigonometric  Polynomial  Bases  of  C[0, 1] 

R.  A.  Lorentz 

Gesellschaft  fur  Mathematik  und  Datenverarbeitung 
St.  Augustin,  Germany 

The  problem  of  finding  a  basis  of  CfO,  1]  which  consists  of  trigonometric  polynomials 
of  lowest  possible  degrees  has  long  been  open.  Privalov  gave  strict  lower  and  upper  bounds 
on  the  growth  of  the  degrees  of  such  a  (growth-)  optimal  basis.  More  recently,  Sahakian 
and  the  author  used  wavelet  techniques  to  construct  trigonometric  bases  that  are  not  only 
optimal  with  respect  to  degree  but  are  also  orthogonal.  This  result  as  well  as  extensions 
will  be  presented. 


Weighted  Approximation  and  Orthogonal 
Polynomials  on  (—00,00) 

D.  S.  Lubinsky 

Witwatersrand  University 
Witwatersrand,  South  Africa 

It  was  S.  N.  Bernstein  who  began  the  theory  of  weighted  approximation  on  (-00, 00), 
with  the  following  question  posed  about  1911:  Let  W  :  R,  — »  [0, 1].  Under  what  conditions 
on  W  is  it  true  that  for  every  continuous  /  :  R  — *•  R,  we  can  find  polynomials  {Pn}  with 

lim  l|(/-Pn)W^llL^(R)=0. 
n— ►oo  * 
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This  qualitative  problem  was  solved  by  Pollard,  Achieser  and  Mergelyan  in  the  1950’s. 
It  was  M.  Dzrbasjan  who  really  began  the  quantitative  side  of  the  theory  in  the  early 
1950’s,  but  he  soon  turned  to  other  topics.  G.  FVeud,  assisted  by  P.  Nevai,  developed  the 
theory  in  the  1970’s  side  by  side  with  work  on  orthogonal  polynomials  for  weights  such  as 
exp(— |x|“),a  >  0.  Fundamental  advances  in  the  subject  were  made  with  the  systematic 
use  of  potential  theory  in  the  1980’s  until  the  present.  We  survey  part  of  the  history  and 
the  most  recent  results,  including  a  selection  from: 

•  Jackson  and  Bernstein  Theorems 

•  Markov-Bernstein  Inequalities 

•  Bounds  on  Orthogonal  Polynomials 

•  Lagrange  Interpolation 

•  Orthogonal  Expansions 

•  A  Unified  Treatment  of  FVeud,  Erdos,  and  non-symmetric  weights 


A  Preliminary  Approach  to  a  Meromorphic  Approximation  in 
A  Strictly  Levi  Pseudoconvex  Domain  with  a  C^-boundary 

C.  H.  Lutterodt*  and  P.  W.  Darko 

Howard  University 
Washington,  D.C. 

Let  y  be  a  complex  manifold  and  suppose  D  is  a  strictly  Levi  pseudoconvex  in  Y .  By 
a  powerful  classical  Hfting  theorem,  we  explore  the  possibilities  of  conveying  an  Oka- Weil 
type  meromorphic  approximation  from  a  Stein  manifold  X  to  D. 


A  Link  Between  Bernstein  Polynomials  and 
Durrmeyer  Polynomials  with  Jacobi  Weights 

Detlef  H.  Mache 

University  of  Dortmund 
Dortmund,  Germany 

We  will  study  a  family  of  linear  operators  P„  defined  by 

n 

fc=0 


with  the  positive  linear  functionals  Ta,k  •  C!{I)  — R  for  A:  =  0, . . . ,  n 


TcAf) 


B(ck  -I-  a  -I- 1,  c(n  -  k)  +  b  +  1)' 
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Here  one  used  for  c  the  special  sequence  in  n  and  a,  with  n€N,  0<a<oo  defined  as 
c:=Cn  :=  [n“]. 

This  family  generalizes  the  classical  Bernstein  operators  and  in  terms  of  the  modulus 
of  continuity  and  the  Ditzian  Totik  modulus  of  smoothness  we  give  among  other  results  a 

•  Voronovskaja  Type  Result, 

•  Global  and  Local  Characterization  Result, 

•  Equivalence  Results  and  Linear  Combination. 

This  fact  of  covering  a  great  class  of  operators  has  been  introduced  in  [6]  (and  was 
generalized  in  [3]).  The  above  polynomial  sequence  {Pn)neN  includes  for  example  a  num¬ 
ber  of  interesting  approximation  processes:  Kantorovic  polynomials  Knf,  the  Bernstein 
polynomials  Bnf,  the  Bemstein-Durrmeyer-Lupas  polynomials  M„/  (J.  J.  Durrmeyer  [2], 
A.  Lupas  [4],  see  also  [1]),  and  the  de  la  Valle-Poussin  type  polynomials  [5]. 
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a  la  theorie  des  moments.  Thfee  de  3e  cycle,  Faculte  des  Sciences  de  I’Universite  de 
Paris,  1967. 

[3]  Gavrea,  I.  and  Mache,  D.  H.,  Generalization  of  Bernstein- Type- Approximation  meth¬ 
ods,  manuscript. 

[4]  Lupas,  A.,  Die  Folge  der  Beta  Operatoren,  Dissertation,  Stuttgart,  1972,  75  pages. 
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An  Efficient  Lower  Bound  for  the  Generalized  Spectral  Radius 

of  a  Set  of  Matrices 

Mohsen  Maesumi 

Lamar  University 

Beaumont,  Texas 

Generalized  spectral  radius  of  a  set  of  matrices  is  a  relatively  new  concept  in  Linear 
Algebra  and  an  important  tool  in  regularity  analysis  of  Wavelets.  Let  E  be  a  collection 
of  m  matrices  and  define  £„  as  the  set  of  all  m”  products  of  length  n  of  the  elements 
of  E.  Then  for  each  A  G  Cn  define  its  average  spectral  radius  by  p{A)  =  [p(A)]^/”. 
Let  Pn(S)  =  maxyi€C„  p{^)  define  the  generalized  spectral  radius  of  E  as  p(E)  = 
limsup„_ooPn(S).  The  cost  of  calculating  p(E)  generally  increases  exponentially  with  the 
required  accuracy.  Here  we  reduce  the  cost  of  Pn(S)  from  m”  matrix  calculations  to  m”/n. 
This  is  accomplished  by  utilizing  the  invariance  of  average  spectral  radius  under  (1)  cyclic 
permutation  of  elements  of  matrix  product,  and  (2)  exponentiation  of  the  product.  The 
equivalence  classes  introduced  by  the  first  relationship  can  be  counted  using  MacMahon’s 
formula,  and  within  the  resulting  representative  classes  the  equivalence  class  introduced  by 
the  second  relationship  can  be  counted  using  Dedekind-Liouville  inversion  principle.  We 
simplify  the  resulting  formula  and  establish  an  inequality  that  limits  the  operation  cost  for 
Pn(S)  to  m”/n. 


Evaluation  of  Kergin  Interpolants 

Ulrike  Maier 

Universitat  Dortmund 
Dortmund,  Germany 

It  is  well  known  that  the  Kergin-operator  C''^(R’')  — >■  is  uniquely  interpolating  in 
p  + 1  nodes,  only.  This  means  that  it  is  interpolating  in  far  less  nodes  than  is  corresponding 
with  the  dimension  of  PJ^. 

In  contrast  to  Kergin  interpolants  to  the  monomials,  Kergin  interpolants  to  arbitrary 
functions  of  C^(R'')  cannot  be  computed  directly  though  there  exists  an  explicit  represen¬ 
tation  of  the  operator  given  by  Micchelli  and  Milman.  The  problem  in  their  representation 
is  that  integrals  of  multiple  directional  derivatives  appear. 

We  give  another  representation  which  is  using  partial  derivatives,  instead,  and  which 
is  reducing  the  problem  to  Kergin  interpolation  of  the  monomials.  Consequently,  Kergin 
interpolation  to  the  monomials  is  investigated.  The  coefficients  of  the  corresponding  inter¬ 
polants  are  presented  by  an  explicit  expression  by  which  we  are  able  to  obtain  recurrence 
relations  for  them. 

Using  results  of  Micchelli  concerning  simplex  splines,  it  is,  after  all,  possible  to 
give  a  computationally  evaluable  representation  of  Kergin  interpolants  for  functions  from 

C'^(RO. 
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On  the  Rate  of  the  Uniform  Rational  Approximation  on  Quasismooth  Arcs 

Victor  V.  Maimeskul 

Institute  for  Applied  Math.  &  Mech. 

National  Ukrainian  Academy  of  Science 
Donetsk,  Ukraine 

Let  L  be  a  bounded  Jordan  arc  in  the  complex  plane  C.  Denote  by  0  <  o  <  1, 

the  class  of  functions  /,  defined  in  L,  which  satisfy  the  Holder  condition 

i/(^)-/(c)i<c(/)iz-cr. 

For  the  best  uniform  approximation  of  continuous  function  /  on  L  by  rational  functions 
(polynomials)  of  degree  at  most  n  we  use  the  notation  Rn{f,  L)  {En{f,  L)). 

If  L=  [0, 1]  is  the  real  interval,  then 

R^(f,L)=0{n-‘‘)  (1) 

for  any  function  /  G  H°‘{L).  This  result  follows  immediately  from  the  same  estimate  for 
En{f,  L)  and  can  not  be  improved,  in  general,  if  a  1. 

Now  let  L  =  [0, 1]  U  [0,  i].  It  is  known  that  /  6  H°‘{L)  implies 

E„(/,i)=0  (2) 

Note  that  the  factor  1/2  here  is  closely  connected  with  the  angle  of  L  at  the  origin,  and 
the  estimate  is  precise  for  the  class  H°‘{L).  As  the  consequence  of  (2),  the  same  relation 
is  valid  for  Rn(f,L),  but  what  about  the  sharpness  of  it  ? 

We  investigate  the  rate  of  the  decrease  of  Rn{f,  L)  for  functions  /  G  11°“  {L),  0  <  a  < 
1,  on  quasismooth  arcs  (arcs  without  “cusps”).  It  turns  out  that  the  non-zero  angles  of 
L  (if  any)  are  not  essential  for  the  rational  approximation  unlike  to  polynomial  one,  and 
estimate  like  (1)  holds  for  an  arbitrary  quasismooth  arc. 


Chziractcrization  of  Best  Approximants  in  Spaces 

Miguel  Marano*  and  Jose  M.  Quesada 

Universidad  de  Jaen 
Jaen,  Spain 

Let  ip  :  [0, 00)  — +  [0, 00)  be  a  convex  function,  y?  ^  0,  (^(0)  =  0,  and  let  (fl.  A,  n) 
be  a  measure  space.  Consider  the  set  =  {/  :  J^ip{a\f\)  <  00  for  all  o  >  0},  where 
/  is  a  real- (complex-)  valued  A-measmrable  function.  Then  C,p  is  a  real- (complex-)  Orlicz 
linear  space.  The  classical  Lp  spaces,  1  <  p  <  00,  are  included  in  this  situation.  For  a 
fixed  /  in  consider  the  set  M^{f)  of  elements  go  that  minimize  fp  tp{\f  —  5I)  among 
the  functions  g  of  some  subclass  M  of  C^.  M,f{f)  is  nonvoid  in  many  contexts.  If  M 
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is  convex,  then  a  characterization  of  go  can  be  obtained  from  the  fact  that  the  function 
G  :  [0, 1]  [0,  oo),G'(A)  =  (p{\f  -  (1  -  Xjgo  -  Ah|)  =  ip{\f  -go-  \{h  -  ^o)|),  is  a 

bounded  convex  function  for  aU  /i  in  jM.  From  this  characterization  we  can  get  a  description 
of  go  for  several  examples  of  /A.  In  this  way,  known  results  in  Lp  approximation  can  be 
extended  to  the  present  approximation.  Moreover,  these  results  can  now  be  explained 
in  terms  of  the  behavior  of  tp.  Analogous  conclusions  are  valid  when  approximating  with 
the  Luxemburg  norm  with  which  any  Orlicz  space  can  be  endowed. 


Quantizing  Multiwavelet  Coefficients 

D.  P.  Hardin,  J.  A.  Marasovich*,  and  D.  Roach 

Vanderbilt  University 
Nashville,  Tennessee 

We  compare  various  quantization  schemes  for  multiwavelet  expansions  using  a  variety 
of  multiwavelets  and  wavelets. 


An  Efficient  Variation  on  the  Oslo  Algorithm 

Avraham  A.  Melkman 

Ben  Gurion  University  of  the  Negev 
Beer  Sheva,  Israel 

Algorithms  for  the  computation  of  the  coefficients  involved  in  knot  insertion  are 
presented  which  are  more  robust  and  efficient  than  the  Oslo  algorithm,  or  its  improvement 
by  Lyche  and  M0orken.  They  are  based  on  a  recurrence  relation  which  is  slightly  different 
from  the  one  for  discrete  B-splines.  This  relation  has  the  advantage  that  it  is  easy  to 
determine  a  priori  which  elements  in  the  recurrence  chain  are  non-zero  using  only  tests  on 
the  indices  of  the  knots. 


TO  BE  ANNOUNCED 
Charles  Micchelli 

IBM  T.  J.  Watson  Research  Center 
Yorktown  Heights,  New  York 
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An  Adapted  Wavelet  Packet  Basis  for  the  Parallel  Solution 
of  Image  Restoration  Problems 

L.  Bacchelli  Montefusco*  and  D.  Lazzaro 

University  of  Bologna 
Bolonga,  Italy 

It  is  well  known  that  the  image  restoration  problem  is  described  by  a  Fredholm 
integral  equation  of  the  first  kind  and  its  approximate  solution  leads  to  a  very  large  ill 
posed  linear  system.  Iterative  regularization  methods  and  parallel  algorithms  are  therefore 
needed  [1],[2],[4].  On  the  contrary,  a  different  strategy  to  solve  such  large-scale  problems 
is  based  on  a  suitable  choice  of  the  approximation  space  basis.  In  fact,  a  proper  choice  can 
be  used  to  convert  the  linear  systems  matrices  to  sparse  form,  as  it  happens,  for  example, 
by  using  orthonormal  wavelet  and  wavelet  packet  bases  [3], (6]. 

In  this  paper  we  look  for  a  least  squares  solution  of  the  restoration  problem,  which 
belongs  to  the  imiform  knots  cubic  spline  space.  We  have  therefore  considered  a  mul- 
tiresolution  non-orthogonaJ  wavelet  packet  decomposition  of  this  space  [5],  which  leads  to 
an  equal  sized  block  form  of  the  original  matrix  [6].  Furthermore,  if  the  kernel  of  the 
integral  equation  (point  spread  function)  is,  as  it  usually  occurs,  of  convolution  type  and 
with  finite  support,  it  is  possible,  under  mild  hypothesis,  to  construct  a  new  compactly 
supported  wavelet  packet  basis  which  is  adapted  to  the  considered  kernel.  On  this  basis 
the  least  squares  matrix  presents  a  nice,  diagonal  block  form,  which  is  a  good  starting 
point  for  the  reafization  of  parallel  algorithms  well  suited  for  coarse  grained  distributed 
memory  parallel  machines. 

Timing  and  efficiency  results  from  a  first  experimentation  carried  out  on  a  i860  hy¬ 
percube  multiprocessor  are  given. 

[1]  A.  N.  Tikhonov,  Solution  of  incorrectly  formulated  problems  and  the  regularization 
method,  Dkl,  AkadNauk  SSSR,  151  (1963),  pp.  501-504. 

[2]  H.  S.  Hou,  H.  C.  Andrews,  Least  squares  image  restoration  using  spline  basis  function, 
IEEE  Tirans.,  Comput.,  vol.  C-26,  n.9  (1977),  pp.  856-873. 
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L°°  Markov-Bernstein  Inequalities  for  Erdos  Weights 

T.  Z.  Mthembu 

University  of  the  North-Qwaqwa 
Phuthaditjhaba,  South  Africa 

Let  W{x)  :=  exp(—Q(x)),x  €  K,  where  Q{x)  is  even,  sufficiently  smooth,  and  is  of 
faster  than  polynomial  growth  at  infinity.  Let  On  denote  the  Mhaskar-Rahmanov- 
Saff  number  for  Q,  and  6n  :=  {jiT{zn))'^ >  1,  where  T  is  an  increasing  function.  For 
|a;|  <  an(H-  2L6n),  we  define 


:=  max  < 


'l-|x| 


On 


+  2L6n 


r(*n)\/5 


+2L«„ 


We  prove  that  for  n  >  1,  and  polynomials  P  of  degree  at  most  n,  there  is  a  C  >  0  such 
that 

||(PW0'»»Wlli.-(s)  <  rr  • 


Climate  Surfaces  for  Texas  Using  Thin  Plate  Splines 

Ranjan  S.  Muttiah*,  Michael  F.  Hutchinson,  and  Paul  T.  Dyke 

Blackland  Research  Center,  Center  for  Resources  and  Environmental  Studies, 

and  Blackland  Research  Center 
Texas  A&M  University,  Australian  National  University 
and  Texas  A&M  University 
Temple,  Texas  and  Canberra 

Thin  plate  splines  and  the  digital  elevation  model  (DEM)  are  used  to  fit  monthly 
average  precipitation  and  temperature  surfaces  for  Texas.  The  elevation  of  the  gauging 
station  is  used  as  an  auxihary  variable  for  interpolation  of  the  precipitation  surfaces.  The 
surfaces  compare  very  favorably  with  existing  climate  contour  maps.  When  surfaces  are  fit 
for  precipitation  without  the  DEM,  large  residual  errors  result.  The  ANUSPLIN  software 
that  was  used  is  easily  learnt  by  a  novice,  allows  a  user  to  easily  locate  stations  that  have 
problems  in  their  datasets,  and  runs  within  a  few  minutes  for  a  grid  of  1000  x  1000  cells  of 
Texas.  This  makes  it  particularly  attractive  for  Geographic  Information  Systems  (GIS). 
The  effects  of  removing  observations  from  input  to  the  spline  model  are  explored. 
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Nonstationary  Spherical  Wavelets  for  Scattered  Data 

F.  J.  Narcowich*  and  J.  D.  Ward 

Texas  A&M  University 
College  Station,  Texas 

We  construct  a  class  of  nonstationary,  orthogonal,  wavelets  generated  by  what  we 
call  spherical  basis  functions  (SBFs),  which  comprise  a  subclass  of  Schoenberg’s  positive 
definite  functions  on  the  m-sphere.  We  will  discuss  decomposition,  reconstruction,  local¬ 
ization,  and  algorithms  for  analyzing  scattered  data  on  the  m-sphere. 


Reproducing  Kernel  Hilbert  Spaces  in  Computational  Harmonic  Analysis 

M.  Z.  Hashed 

University  of  Delaware 
Newark,  New  Jersey 

Reproducing  kernel  Hilbert  spaces  (RKHS)  play  an  important  role  in  several  areas  of 
approximation  theory,  numerical  analysis,  and  estimation  theory.  In  recent  years,  the  use 
of  RKHS  has  led  to  interesting  interactions  between  harmonic  analysis,  signal  processing, 
and  computational  mathematics.  In  this  talk,  I  will  indicate  some  results  of  the  interac¬ 
tion,  with  emphasis  on  sampling  expansions,  integral  equations,  and  numerical  conformal 
mappings. 


Sparse  Approximate  Multiquadric  Interpolation 

B.  K.  Natarajan*  and  R.  E.  Carlson 
Hewlett  Packard  Laboratories  and  Lawrence  Livermore  Labs 


Multiquadric  interpolation  is  a  technique  for  interpolating  scattered  samples  of  mul¬ 
tivariate  functions,  typically  to  support  data  visualization.  We  present  an  algorithm  for 
computing  sparse  but  approximate  multiquadric  interpolants.  Such  interpolants  are  useful 
since  (1)  the  cost  of  evaluating  the  interpolant  scales  directly  with  the  number  of  coeffi¬ 
cients,  and  (2)  the  principle  of  Occam’s  Razor  suggests  that  the  interpolant  with  fewer 
coefficients  better  approximates  the  underlying  function. 
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Spherical  Splines 

P.  Alfeld,  M.  Neamtu*,  and  L.  Schumaker 

Vanderbilt  University 
Nashville,  Tennessee 

We  present  a  natural  way  to  define  spaces  of  splines  on  the  sphere  and  on  sphere-like 
surfaces,  which  is  based  on  a  theory  of  spherical  Bernstein-Bezier  polynomials.  Such  spaces 
exhibit  most  of  the  important  properties  of  the  classical  BB-polynomials  on  planar  trian¬ 
gulations.  The  methods  discussed  have  immediate  applications  to  a  variety  of  important 
practical  problems  including  interpolation,  data  fitting,  and  modelling  of  surfaces  defined 
on  the  sphere. 


Perturbation  of  Orthogonal  Polynomials  on 
An  Arc  of  the  Unit  Circle 

Paul  Nevai*  and  Walter  Van  Assche 

Ohio  State  University,  Columbus,  Ohio  and 
Katholieke  Universiteit  Leuven,  Leuven,  Belgium 

Orthogonal  polynomials  on  the  unit  circle  axe  completely  determined  by  their  reflec¬ 
tion  coefficients  through  the  Szego  recurrences.  We  assume  that  the  reflection  coefficients 
converge  to  some  complex  number  a  with  0  <  |a|  <  1.  The  polynomials  then  live  essen¬ 
tially  on  the  arc  {e*^  :  a  <  ^  <  2tt  -  a}  where  cos  f  y/l  -  |a|2  with  a  €  (0,7r).  We 
analyze  the  orthogonal  polynomials  by  comparing  them  with  the  orthogonal  polynomi¬ 
als  with  constant  reflection  coefficients,  which  were  studied  earlier  by  Ya.  L.  Geronimus 
and  N.  1.  Akhiezer.  In  particular,  we  show  that  under  certain  assumptions  on  the  rate 
of  convergence  of  the  reflection  coefficients  the  orthogonality  measure  will  be  absolutely 
continuous  on  the  arc  with  a  continuous  weight  function. 


Non-Stationary  Wavelet  Bases  for  Sobolev  Spaces 

Igor  Ya.  Novikov 

Voronezh  State  University 
Voronezh,  Russia 

An  approach  for  the  construction  of  non-stationary  orthonormal  wavelets  is  pre¬ 
sented.  The  method  allows  us  to  obtain  non-stationary  orthonormal  infinitely  differentiable 
compactly  supported  wavelets.  It  is  proved  that  these  wavelets  form  a  Riesz  basis  for  all 
Sobolev  spaces  simultaneously.  The  estimates  of  the  uncertainty  bounds  of  non-stationary 
wavelets  will  be  discussed. 


68 


On  the  Peculiarities  of  Canonical  Injection  Operators 
for  Rearrangement  Invariant  Spaces 

S.  Ya.  Novikov 

Samara  State  University 
Samara,  Russia 


(Abstract  not  avaialble) 


Bivariate  Spline  Interpolation  and  Approximation  Order 

Gunther  Niirnberger 

University  of  Mannheim 
Mannheim,  Germany 

We  describe  a  general  method  for  constructing  point  sets  which  admit  unique  La¬ 
grange  and  Hermite  interpolation  by  spaces  of  bivariate  splines  of  arbitrary  degree  and 
smoothness.  The  splines  are  defined  on  uniform  type  triangulations  of  simply  connected 
domains.  It  is  shown  that  Hermite  interpolation  by  splines  of  smoothness  one  yields 
(nearly)  optimal  approximation  order. 

Numerical  examples  using  2000  Lagrange  interpolation  points  are  given.  For  comput¬ 
ing  the  interpolating  splines,  only  small  systems  have  to  be  solved.  Interpolation  methods 
for  spline  spaces  of  special  degree  are  known  in  the  literature. 


On  the  Number  of  Orthogonal  Wavelet  Multiplier  Matrices 

Gerhard  Opfer 

University  of  Hamburg 
Hamburg,  Germany 

Scaling  equations  have  usually  the  form 

^  CkV’CAx  -  k) 
k€Z" 

where  A  is  a  square,  nonsingulax  matrix  of  order  n  with  integer  entries  and  therefore, 
det(A)  is  a  non  vanishing  integer.  In  addition,  it  is  required  that  A  has  a  certain  stretching 
property  expressed  by  |A(A)|  >  1  which  should  mean  that  all  eigenvalues  of  A  are  located 
outside  the  closed  unit  disk.  Since  the  determinant  is  the  product  of  the  eigenvalues,  it 
follows  that  |det(A)|  >  2.  In  the  applications  the  matrices  considered  are  orthogonal  in 
the  sense  that  AA^  =  cl  where  I  is  the  identity  matrix  and  c  a  certain  constant.  One 
is  usually  interested  in  matrices  A  where  the  absolute  value  of  the  determinant  is  small, 
because  |det(A)|  —  1  is  the  number  of  wavelets  necessary  to  span  L2(R”).  Our  lecture 
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will  mainly  be  concerned  with  the  number  of  matrices  with  the  prescribed  properties  when 
either  |  det(A)|  is  prescribed  or  is  bounded  by  a  given  constant.  There  will  be  results  for 
general  n  and  some  explicit  counting  results  forn  =  3  and  n  =  4. 


Piecewise  Linear  Prewavelets  of  Small  Support 

Uwe  Kotycka  and  Peter  Oswald* 

Priedrich-Schiller  Universitat  Jena 
Jena,  Germany 

For  preconditioning  elliptic  problems  in  Sobolev  spaces  over  domains,  multi¬ 
level  frames  composed  of  finite  element  nodal  basis  functions  have  gained  popularity  over 
the  last  five  years.  For  certain  applications,  the  construction  of  L2-semiorthogonal  pre¬ 
wavelets  related  to  these  nodal  basis  frames  is  desirable.  We  give  first  results  in  the  special 
case  of  piecewise  linear  finite  element  functions  on  sequences  of  triangulations  in  the  two- 
dimensional  case.  One  topic  considered  is  the  minimization  of  the  size  of  the  underlying 
“refinement”  masks. 


On  Continuity  of  the  Prox  Map  and  the  Restricted  Center  Map 

D.  V.  Pai*  and  B.  M.  Deshpande 

Indian  Institute  of  Technology,  Bombay 
Powai,  Bombay,  INDIA 

In  this  paper  we  review  some  of  the  recently  introduced  hypertopologies  on  collections 
of  nonempty  closed  subsets  of  a  metrizable  space  in  the  context  of  continuity  properties 
of  the  prox  map  and  the  restricted  center  map.  We  identify  here  suitable  families  A,  B  of 
nonempty  closed  sets  equipped  with  natural  topologies,  as  coordinate  spaces  for  the  gap 
functional  and  the  restricted  radius  functional,  with  a  view  to  explore  bivariate  continuity 
of  these  functionals.  This  leads  us  to  sharpenings  of  many  known  results  as  well  as  to 
some  new  results  for  upper  semicontinuity  of  the  prox  map  and  the  restricted  center  map. 
This,  in  turn,  also  leads  us  to  improvements  of  certain  best  approximation  results  for 
convex-valued  multifunctions  and  fixed  point  theorems  for  such  multifunctions. 


An  Equivalence  Relation  Between  Multiresolution  Analyses  of  ^^(R) 

Manos  Papadakis*,  Theodoros  Stavropoulos,  and  N.  Kalouptsidis 

Hellenic  Military  Academy 
Athens,  Hellas  (Greece) 

Classification  problems  are  of  major  importance  in  various  fields  of  mathematics.  We 
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present  a  natural  notion  of  equivalence  called  unitary  equivalence  between  multiresolution 
analyses  of  L^(R). 

Let  M  =  {Vj  :  j  €  Z}  be  a  multiresolution  analysis  of  I'^(R)  (MRA)  and  (p  in  I-^(R). 
We  say  that  v?  is  a  scaling  function  for  M  if  {(po,n  :  n  e  Z}  is  an  orthonormal  basis  of  Fq- 

Definition  1.  Two  MRAs  AA  and  on  L^(R)  are  called  unitarily  equivalent  if  there 
exist  scaling  functions  p  and  x  for  M  and  Af  respectively  and  a  unitary  operator  U  defined 
on  X^(R)  such  that  U(pm,n  =  Xm,n,  m,n  e  Z. 

Employing  operator  theoretic  techiques,  we  can  prove  the  following  characterization 
of  the  unitary  equivalence. 

Theorem  2.  Let  Ai  and  ff  he  two  MRAs.  Then  they  are  unitarily  equivalent  iff  there 
exist  scaling  functions  ip  for  M.  and  x  for  Af  such  that  x(t)  =  h{t)<p{t)  a.  e.  ,  where  h  is 
a  Borel  measurable  function  defined  on  R  such  that  |h(t)|  =  1  and  h{2t)  =  h{t)  a.  e. 

The  equivalence  relation  we  introduced  is  naturally  associated  with  the  well-known 
Decomposition  and  Reconstruction  algorithms  of  Mallat  defined  for  every  MRA. 

Theorem  3.  Let  Ai  and  M  he  two  MRAs  of  L^(R).  Then  Ai  and  Af  are  unitarily 
equivalent  iff  they  give  the  same  Decomposition  and  Reconstruction  algorithms. 

Theorem  3  allows  us  to  obtain  a  partial  generalisation  of  the  uniqueness  results  of  I. 
Daubechies  and  J.  C.  Lagarias  concerning  the  solutions  of  two-scale  difference  equations. 
Theorem  2  provides  a  tool  to  examine  whether  any  two  given  MRAs  are  unitarily  equiv¬ 
alent.  The  conditions  for  the  fimction  h  of  Theorem  2  seem  rather  restrictive.  There  are 
certain  cases  where  the  equivalence  classes  are  singletons.  Focusing  to  the  functions  h  which 
are  constant  in  each  of  the  intervals  (— oo,0)  and  [0, -l-oo),  we  obtain  equivalences  imple¬ 
mented  by  certain  linear  combinations  of  the  identity  operator  and  the  Hilbert  Transform. 
Finally  we  prove  that  MRAs  defined  by  R-splines  of  different  degree  are  non-equivalent. 


On  the  Best  Approximation  of  Abstract  Integrable  Functions 

Irina  Peterburgsky 
Suffolk  University 
Boston,  Massachusetts 

Let  E  be  a  conjugate  complex  Banach  space,  L  be  a  Banach  space  of  Bochner  inte¬ 
grable  functions  from  the  unit  circumference  of  a  complex  plane  to  E,  and  A”  be  a  subspace 
of  L  that  consists  of  functions  with  zero  Fourier  coefficients  of  non-positive  indices.  For 
any  g  belonging  to  L,  a  problem  about 

inf  11^  -  h\\L,  h  from  K,  (1) 

h 

is  considered.  The  existence  of  the  best  approximation  in  (1)  has  been  established  based 
on  specific  properties  of  Banach-space-valued  functions  obtained  in  the  paper. 
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Polynomial  Interpolation  and  Peano  Kernels 
Knut  Petras 

Institut  fiir  Angewandte  Mathematik 
Braunschweig,  Germany 

A  description  of  the  Peano  kernel  method  for  the  estimation  of  errors  of  polynomial 
interpolation  is  given.  The  dependence  of  the  error  bounds,  which  use  the  first  derivative, 
on  the  location  of  the  evaluation  point  is  characterized.  The  results  are  applied  to  poly¬ 
nomial  interpolation  involving  Chebyshev  nodes.  Sharp  bounds  for  the  error  estimation 
with  the  first  derivative  as  well  as  asymptotic  estimates  involving  higher  derivatives  are 
derived. 

Similar  results  as  in  the  algebraic  case  are  also  proved  for  trigonometric  interpolation 
involving  equidistant  nodes. 


Approximation  Properties  of  Neural  Networks 

Ron  DeVore,  Konstantin  Oskolkov,  Pencho  Petrushev*,  and  Vladimir  Temlyakov 

University  of  South  Carolina 
Columbia,  South  Carolina 

P.  Petrushev 

Bulgarian  Academy  of  Sciences 
Sofia,  Bulgaria 

Neural  networks  are  being  proposed  as  an  efficient  approximation  method  in  several 
applications.  However,  their  approximation  properties  are  not  well  understood.  This  talk 
will  present  recent  results  on  the  degree  of  approximation  possible  from  linear  and  nonlinear 
spaces  generated  by  sigmoidal  functions  and  other  ridge  functions.  In  particular,  we  shall 
show  that  it  is  possible  to  obtain  optimal  approximation  order  for  classical  function  spaces 
firom  certain  linear  spaces  generated  by  sigmoidal  functions. 


Nonorthogonal  Compactly  Supported  Wavelets 
Alexander  P.  Petukhov 

St.  Petersburg  Branch  of  Steklov  Math.  Institute 
St.  Petersburg,  Russia 

We  consider  here  only  the  1-periodic  case,  but  most  of  the  facts  are  valid  for  the 
nonperiodic  one.  Let  T  =  [0, 1] ,  the  sequence  of  spaces  Vo  C  Vi  •  •  •  C  C  .  •  •  generate 
a  Multiresolution  approximation  (MRA)  of  L^(T),  the  corresponding  translates  of  the 
function  form  a  basis  (in  general  nonorthogonal)  of  Vj,  Vj+i  =  IVo  ©  Wi  ©  •  •  •  0  IVj,  j  = 
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0, 1, 2, . . . ,  be  an  orthogonal  sum  of  wavelet  spaces, 


-  =  /' 


<p^(x)  ■  -  i/2N)dx,  i  =  0, 1, . . . ,  2iV  -  1,  where  N  =  2^. 


Proposition.  Suppose 

2N-1 

^{x  +  l/2Ar)  =  Y,  +  n/2Ny 

n=0 


Then  the  system  {rj)^{x  -  k/N)}^^Q  of  functions  forms  a  basis  of  the  space  . 

In  particular,  it  follows  from  the  Proposition  that  if  supp  ifJ  <  p/N  for  all  j  (it 
corresponds  to  the  case  of  a  finite  support  for  L^R)),  then  supp  ^  <  2p/N .  Hence,  it  is 
easy  to  construct  algorithms  of  decomposition  and  reconstruction  for  elements  of  which 
have  the  complexity  0{pN)  (the  same  as  for  the  case  of  an  orthonormal  basis). 

Unfortunately,  the  mentioned  algorithm  has  dealt  with  solving  of  systems  of  linear 
equations,  so  it  cannot  be  used  in  the  real  time  mode.  However,  if  we  do  not  require 
orthogonality  of  the  spaces  and  we  only  suppose  that  V^+i  =  +  •  •  •  +  Wj, 

where  the  sum  is  direct,  then  there  exists  a  decomposition  of  Vj  of  such  type  for  which 
a  decomposition  algorithm  can  be  realized  in  real  time  as  a  discrete  convolution  with  a 
finite  window.  If  at  the  same  time  two-scale  relations  between  the  spaces  Vj  take  place, 
then  an  algorithm  of  reconstruction  can  be  also  realized  a  discrete  convolution  with  a  finite 
window. 

Besides,  we  establish  a  natural  one-to-one  correspondence  between  all  possible  MRA 
of  L^(T)  and  functions  p,  €  L^(T)  for  which  p{2k  for  any  k  e7i. 


Perturbation  of  Orthogonal  Polynomials 
On  an  Arc  of  the  Unit  Circle 

Leonid  Golinskii,  Paul  Nevai,  Walter  Van  Assche,  and  Prank  Pinter* 

Ohio  State  University 
Columbus,  Ohio 

Katholieke  Universiteit  Leuven 
Leuven,  Belgium 

The  perturbation  analysis  of  orthogonal  polynomials  on  an  arc  of  the  unit  circle 
was  initiated  by  L.  Golinskii,  P.  Nevai,  and  W.  Van  Assche,  and,  simultaneously  by  F. 
Peherstorfer  and  R.  Steinbauer.  They  utilized  the  one-step  matrix  recurrence  relation  to 
derive  asymptotic  results  about  the  perturbed  system  of  orthogonal  polynomials  and  to 
describe  the  new  measure  of  orthogonality. 

This  is  a  joint  work  with  L.  Golinskii,  P.  Nevai,  and  W.  Van  Assche.  The  purpose  of 
this  talk  is  to  obtain  further  asymptotic  results  and  conclusions  about  the  measure  using 
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different  techniques.  We  use  the  less  obvious  three-term  recurrences  along  with  former 
techniques  developed  by  P.  Nevai,  T.  S.  Chihara  and  A.  Mate  for  orthogonal  polynomials 
on  the  real  line;  as  well  as  incorporate  the  decomposition  of  linear  recurrences  to  place  our 
results  into  much  more  general  context. 


Real  and  Complex  Interpolation  of  Operator  Spaces 

Gilles  Pisier 

Texas  A&M  University,  College  Station,  Texas,  and 
Universite  Paris,  Paris  France 

The  talk  will  discuss  several  results  on  the  real  and  complex  interpolation  method 
applied  to  couples  of  operator  spaces  of  the  form  {B{Xo),B{Xi))^  where  (Xo,Xi)  is  a 
couple  of  Banach  spaces  and  B(X)  denotes  the  space  of  all  bounded  linear  operators  on 
X.  The  main  example  is  the  case  Xq  =  Tooj  Xi  =  Li.  In  that  case  we  are  able  to 
describe  the  real  and  complex  interpolation  spaces  and  to  compute  an  equivalent  of  the 
fft-functional  of  Peetre.  In  the  last  quarter  of  the  talk  we  will  discuss  the  non-commutative 
case,  i.e.  the  case  when  ATo  (resp.  Xi)  is  a  “non-commutative”  Loo-space  (resp.  Li-space), 
i.e.  a  von  Neumann  algebra  (resp.  the  predual  of  a  von  Neumann  algebra). 


Multi  wavelets  With  Short  Support 

Gerlind  Plonka 

University  of  Rostock 
Rostock,  Germany 

We  want  to  generalize  the  concept  of  cardinal  wavelets  in  the  following  way:  We 
introduce  the  shift — ^invariant  sample  samples  Vj  of  level  j  spanned  by  a  finite  set  of 
linearly  independent  (orthogonal)  scaling  functions  (pk  (A:  =  0, . . . ,  r  —  1) 

Vj  :=  span{(p*;(2-^  •  — f);  A:  =  0, . . . ,  r  —  1,  f  €  Z} 


with  r  G  N. 

The  necessary  and  sufficient  conditions  for  obtaining  a  generalized  multiresolution 
analysis  {Vjjjez  will  be  analyzed.  In  particular,  a  matrix  refinement  equation  for  the 
scaling  functions  ipk  (A:  =  0, . . . ,  r  —  1)  must  be  satisfied.  It  is  possible  to  construct  scaling 
functions  that  have  extra  properties  which  can  not  be  achieved  with  only  one  scaling 
function  (cf.  Daubechies).  The  new  scaling  functions  can  be  required  to  be  symmetric 
(linear  phase)  and  with  short  support  while  their  translates  form  an  orthogonal  family. 
Furthermore,  we  want  to  obtain  a  certain  approximation  order  p  of  Vj.  Some  results 
by  Donovan-Geronimo-Hardin-Massopust  and  Strang-Strela  for  the  case  r  =  2  can  be 
generalized. 
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Then  our  next  step  is  to  construct  associated  wavelets  ipk  (A:  =  0, . . . ,  r  —  1)  by  finding 
an  orthonormal  basis  of  the  shift-invariant  wavelet  space 


Wj  :=  Vj+i  ©  V,- 


in  the  form 

Wj  =  span{^fc(2-'  •  —^)  :  A:  =  0, . . .  ,r  —  1,  £  e  Z}. 

Our  main  goal  is  the  construction  of  new  wavelets  with  shortest  support.  The  very  short 
support  of  wavelets  opens  new  possibilities  for  applications,  for  instance  in  numerical  so¬ 
lution  of  operator  equations. 


Interpolation  and  Incidence  Geometry 
Burkard  Bolster 

University  of  Canterbury 
Christchurch,  New  Zealand 

Around  1950  people  on  both  sides  of  the  pure/applied  divide  in  mathematics  started 
to  work  on  very  similar  problems.  On  the  one  hand,  people  interested  in  interpolation 
theory,  approximation,  and  convexity  started  to  explore  non-linear  systems  of  interpolating 
functions.  On  the  other  hand,  incidence  geometers  started  to  investigate  topological  affine 
planes,  projective  planes,  Laguerre  planes,  etc.  Over  the  past  forty  years  both  areas  have 
been  developed  considerably,  though  there  has  been  little  or  no  communication  between 
the  two  groups. 

I  will  describe  some  of  the  fundamental  objects  of  interest  to  both  groups  from  both 
possible  points  of  view.  I  will  also  talk  about  some  recent  progress  that  has  been  made  by 
bringing  together  results  from  both  areas. 


Fast  Poisson  Solvers  in  Shift-Invariant  Spaces 

Gisela  Poplau 

University  of  Rostock 
Rostock,  Germany 


Using  periodic  shift-invariant  spaces,  we  present  a  new  approach  to  fast  Poisson 
solvers.  Our  consideration  generalizes  the  known  concept  of  fast  Poisson  solvers.  A  con¬ 
venient  tool  to  construct  new  fast  Poisson  solvers  is  given.  This  method  is  demonstrated 
in  periodic  shift-invariant  spaces  spanned  by  translates  of  box  splines. 

Further,  we  obtain  error  estimates  for  the  approximate  solution  of  the  Poisson  equa¬ 
tion  in  periodic  shift-invariant  spaces.  Strang-Fix  conditions  are  essential  for  these  error 
estimates.  The  close  relation  to  interpolation  in  periodic  shift-invariant  spaces  is  investi¬ 
gated. 
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The  Fourier  technique  provides  very  efficient  niunerical  algorithms  based  on  fast 
Fourier  transforms.  Numerical  tests  using  box  splines  will  be  presented. 


Polynomial  Bases  for  Spaces  of  Continuous  Functions 

Jurgen  Prestin 
Universitat  Rostock 
Rostock,  Germany 

It  is  the  aim  of  the  talk  to  investigate  orthogonal  polynomial  bases  for  the  spaces  C2-n 
or  C[-l,  1].  Of  special  interest  are  bases  {pn}  where  p„  is  a  polynomial  of  small  degree. 
Recently,  R.  Lorentz  and  A.  A.  Sahakian  constructed  such  an  orthogonal  basis  consisting 
of  trigonometric  polynomials  of  optimal  degree. 

Here  we  consider  finite  dimensional  nested  spaces  of  trigonometric  polynomials  con¬ 
structed  from  de  la  Vallee  Poussin  means  of  the  Dirichlet  kernel.  Following  an  approach  of 
A.  A.  Privalov  (1991)  we  investigate  the  corresponding  Multiresolution  Analysis.  The  scal¬ 
ing  functions  and  wavelets  are  given  explicitly  as  trigonometric  fundamental  interpolants 
and  decomposition  and  reconstruction  algorithms  can  be  described  in  simple  matrix  no¬ 
tation.  The  circulant  structure  of  all  relevant  matrices  allows  the  use  of  Fast-Fourier- 
TVansform  techniques  for  the  actual  implementation.  Thus  we  achieve  almost  optimal 
complexity  compared  to  other  wavelet  approaches  derived  from  implicit  two-scale  rela¬ 
tions,  while  dealing  with  a  fully  computable  trigonometric  multiresolution  analysis  with 
explicit  algebraic  formulas. 

Furthermore  we  describe  corresponding  wavelet  packets  which  yield  refined  frequency 
localization  properties  and  can  be  used  for  the  direct  construction  of  certain  orthogonal 
polynomial  bases  of  C2Tr  with  optimal  degree.  For  the  corresponding  partial  sum  operators 
we  obtain  estimates  with  explicit  constants. 

The  special  structure  of  the  underlying  de  la  Vallee  Poussin  means  allows  to  transform 
most  of  these  results  to  the  algebraic  case.  In  particular  we  obtain  algebraic  polynomial 
wavelet  bases  for  Cf— 1,1],  where  certain  interpolation  conditions  for  the  wavelets  are 
satisfied.  Here  the  orthogonality  is  with  respect  to  the  Chebyshev  weight. 


76 


Extensions  of  a  Theorem  of  Ky  Fan 

Joao  B.  Prolla*,  Ary  O.  Chiacchio,  and  Maria  Sueli  M.  Roversi 

IMECC-UNICAMP 
Campinas,  Brazil 

Let  M  be  a  compact  and  convex  non-empty  subset  of  a  normed  space  E.  Let  /,  g 
and  h  be  three  functions  from  M  into  E  such  that  |lp(a;)  —  /(x)||  <  ||/i(x)  -  /(x)||  for  all 
X  G  M.  We  prove  that  when  /  and  g  are  continuous  and  h  is  almost  quasi-convex,  there 
exists  xo  E  M  such  that  ||5(xo)  -  /(xo)||  <  \\z  -  /(xo)|l  for  all  z  €  S{h{My,g{xo)),  where 
for  any  set  X  C  E  and  v  €  E,  S(M;  v)  is  the  set  of  all  z  G  E  of  the  form  2  =  At  -f-  (1  -  X)v 
ioT  t  E  X  and  A  >  1.  The  same  result  holds  for  the  weah  topology  if  we  assume  /(M) 
relatively  compact.  As  corollaries  we  get  coincidence  and  fixed  point  results,  some  of  them 
from  the  literature  and  several  new  ones.  The  main  tool  used  is  a  lemma  due  to  Ky  Fan. 


Recognition  of  Biological  Shapes  by  Means  of  A  Wavelet  Analysis 

M.  Cotronei  zind  L.  Puccio* 

University  of  Messina 
Messina,  Italy 

Most  of  the  research  efforts  in  quantitative  morphology  deals  with  shape  determina¬ 
tion  of  the  objects  (e.g.,  biomedical  forms)  under  investigation.  Recently,  we  have  faced 
this  kind  of  problem  with  an  analytical  approach  using  an  algorithm  based  on  spline  ap¬ 
proximation.  In  this  work  we  present  a  method  for  the  assessment  of  form  descriptors  by 
means  of  wavelet  and  wavelet-packet  analysis.  We  use  a  multivariate  unsupervised  learning 
scheme  (cluster  analysis)  as  objective  estimator  of  the  recognition  process.  The  classifier 
objects,  taken  as  input  data,  are  some  fundamental  geometrical  shapes. 


Trigonometric  Wavelets  for  Hermite  Interpolation 

Ewald  Quak 

Center  for  Approximation  Theory 
Texas  A&M  University 
College  Station,  Texas 

In  this  talk,  a  multiresolution  analysis  of  nested  subspaces  of  trigonometric  poly¬ 
nomials  is  investigated.  The  pair  of  scaling  functions  which  span  the  sample  spaces  are 
fundamental  functions  for  Hermite  interpolation  on  a  dyadic  partition  of  nodes  on  the 
interval  [0, 27r).  Two  wavelet  functions  that  generate  the  corresponding  orthogonal  com¬ 
plementary  subspaces  are  constructed  so  as  to  possess  the  same  fundamental  interpolatory 
properties  as  the  scaling  functions.  Together  with  the  corresponding  dual  functions,  the 
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interpolatory  properties  of  the  scaling  functions  and  wavelets  are  used  to  find  the  specific 
decomposition  and  reconstruction  sequences  of  this  trigonometric  multiresolution  analysis. 


A  Degree  Estimate  for  Curvature  Continuous  Subdivision  Surfaces 

Ulrich  Reif 

University  of  Stuttgart 
Stuttgart,  Germany 

Subdivision  schemes  like  the  Catmull-Clark  algorithm  admit  the  generation  of  smooth 
surfaces  from  control  meshes  of  arbitrary  topological  type.  Recently,  this  more  or  less  em¬ 
pirical  fact  was  verified  by  the  speaker.  A  short  survey  of  these  results  providing  sufficient 
conditions  for  convergence  to  tangent  plane  continuous  limit  surfaces  is  given.  When  seek¬ 
ing  algorithms  converging  to  curvature  continuous  surfaces,  which  are  unknown  so  far,  the 
following  estimate  has  to  be  taken  into  account:  Consider  an  algorithm  generating  limit 
surfaces  which  consist  of  piecewise  polynomial  patches  of  bi-degree  d  joining  parametri¬ 
cally  smooth  of  order  k.  Further,  assume  that  the  algorithm  is  cyclic,  i.e.,  all  sides  of  an 
7i-sided  configuration  are  treated  equally,  and  that  it  is  not  degenerated,  i.e.,  there  are  no 
enforced  flat  spots.  Then -the  limit  surface  can  be  curvature  continuous  only  if  d  >  2A:  -f  2. 
In  particular,  the  minimal  degree  is  dmin  =  2A:inin  +  2  =  6.  This  result  applies  even  to 
non-stationary  and  non-compactly  supported  schemes. 


Piecewise  Convex/Concave  Function  Estimation 
and  Model  Selection 

Kurt  S.  Riedel 

Courant  Institute  of  Mathematical  Sciences 
New  York,  New  York 

We  consider  estimation  problems  where  the  function  is  known  to  consist  of  a  small 
number  of  pieces  on  which  it  is  strictly  concave  or  convex.  To  estimate  the  function  we 
need  to  identify  these  regions  and  perform  a  constrained  fit.  The  function  estimation 
problem  consists  of  three  steps:  1)  Given  the  data,  fit  the  best  piecewise  convex  function 
to  the  data;  2)  Given  the  number  of  intervals,  estimate  the  change  point  locations;  and  3) 
Determine  the  number  of  change  point  locations.  We  refer  to  the  last  step  as  the  model 
selection  problem. 

We  consider  a  function  g{t)  measured  at  N  locations  {ti,i  =  1 . . .  N}.  The  measure¬ 
ments  are  noisy:  yi  =  ^(ti)+Ct  where  the  Cj  are  independent  random  variables  with  variance 
[cj]  =  (T^ .  Initially,  we  seek  an  estimate  only  at  the  measurement  points:  gi  =  g{ti)-  If 
gi{t)  is  convex,  it  satisfies  the  constraints 


aig{ti-i)  +  {ai  +  ai-i)g{ti)  +  ai-ig{tt+i)  >  0 


(1) 
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where  a*  =  U^i  —  to.  We  begin  by  fitting  g(t)  assuming  that  g{t)  is  convex.  This  gives 
N  —  2  linear  constraints  on  the  least  squares  problem: 

Min  II  y  -  g  IP  subject  to  Ag  >  0  (2) 

where  A  is  the  tridiagonal  matrix  with  entries  o*.  We  denote  the  mean  square  error  of 
the  constrained  fit  by  Gq.  We  then  try  all  possible  fits  with  one  change  from  convex  to 
concave. 

In  general,  we  allow  K  changes  and  consider  choices  of  possible  change  points. 

We  denote  the  candidate  set  of  change  points  by  Si{w) . . .  iS'fc(iy).  We  then  minimize  the 
least  squares  error  subject  to  the  constraints 

_l'yi(«')[a.^._i  _|.  (a-  4.  +  ai-  l^<+i]  >  0  , 

so  where  'yi{w)  =  1  if  S2t{w)  <ti  <  821+1  and  7i(iu)  =  0  if  S2e+i{w)  <ti<  821+2  for  each 
set  of  change  points.  This  is  a  quadratic  programming  problem  and  there  is  an  unique 
solution.  The  minimization  problem  may  be  solved  exactly  using  “active  set”  methods  in 
quadratic  programming  or  may  be  solved  approximately  using  projected  successive  over 
relaxation  (PSOR)  iterations.  The  PSOR  method  uses  the  band  structure  in  A{w)  while 
the  active  set  programs  require  modification  to  talce  advantage  of  the  tridiagonal  structure. 
In  both  cases,  the  following  dual  formulation  is  easier  to  implement: 

Ao  =  argminA  ||  A^A  +  y  \\^  subject  to  A  >  0 

where  is  the  {N  —  2)  x  N  constraint  matrbc  whose  i-th  row  is  given  by  (2).  The 
principal  advantage  of  the  dual  formulation  is  that  the  positivity  constraints  are  easier  to 
implement  than  the  three  point  difference  inequality  constraints  of  (2).  In  principle,  all 
2  (•^^)  cases  need  to  be  evaluated  to  find  the  best  fit.  Since  this  is  computationally  costly, 
we  consider  suboptimal  methods  which  examine  only  a  subset  of  the  possible  change  points. 
These  candidate  change  points  are  determined  by  1)  initializing  with  a  kernel  smoother, 
2)  examining  long  runs  of  data  points  where  the  constraints  are  simultaneously  active. 

Given  the  best  fit  with  K  changes  firom  convex  to  concave  and  mean  square  error  cr| 
and  the  best  fit  with  K  +  1  change  points,  we  need  to  decide  which  model  is  preferable. 
We  examine  a  variety  of  different  model  selection  criteria  including  information  criteria 
and  Bayesian  criteria. 


Smooth  Refinable  Functions  Provide  Good  Approximation  Orders 

Amos  Ron 

University  of  Wisconsin 
Madison,  Wisconsin 

We  apply  the  general  theory  of  approximation  orders  of  principal  shift-invariant 
spaces  of  de  Boor,  DeVore,  Ron,  TAMS,  1994,  to  the  special  case  when  the  generator 


79 


<f>  of  the  underlying  space  satisfies  an  n-scale  relation  (i.e.,  is  a  “father  wavelet”).  We 
show  that  the  approximation  orders  provided  by  such  principal  shift-invariant  spaces  are 
bounded  from  below  by  the  smoothness  class  of  the  scaling  function,  as  well  as  by  the 
decay  rate  of  its  Fourier  transform.  Specifically,  it  is  shown  that,  under  some  mild  technical 
conditions  on  the  scaling  function  <f>,  approximation  order  k  is  provided  if  either  (f)  lies  in 
the  Sobolev  space  W* or  its  Fourier  transform  ${w)  decays  near  oo  like  o(|u;|*’“^). 

The  first  class  of  results  (those  that  are  concerned  with  the  condition  <j)  E 
improve  previously  known  results  of  Yves  Meyer  and  of  Cavareta-Dahmen-Micchelli. 

The  techniques  apply  to  non-stationary  ladders  as  well,  as  well  as  to  ladders  generated 
by  several  functions. 

If  you  are  aware  of  further  known  results  on  this  topic,  please  share  this  information 
with  the  author. 


Rational  Approximation  with  Varying  Weights  I 

P.  Borwein,  E.  A.  Rakhmanov,  and  E.  B.  Sail* 

Institute  for  Constructive  Mathematics 
University  of  South  Florida 
Tampa,  Florida 

We  investigate  two  problems  concerning  uniform  approximation  by  weighted  rationals 
where  Tn  =  Pn/Qn  is  a  rational  function  of  order  n.  Namely,  for  w{x)  := 
we  prove  that  uniform  convergence  to  1  of  u;”r„  is  not  possible  on  any  interval  [0,  a]  with 
a  >  27r.  For  w(x)  :=  x^,  9  >  1,  we  show  that  uniform  convergence  to  1  of  «;”rn  is 
not  possible  on  any  interval  [6,1]  with  6  <  tan^(7r(0  —  1)/4B).  (The  latter  result  can  be 
interpreted  as  a  rational  analogue  of  results  concerning  “incomplete  polynomials”.)  More 
generally,  for  a  >  0,  ^  >  0,  a  -I-  >  0,  we  investigate  for  w{x)  =  e"'  and  wix)  =  x®, 

the  possibility  of  approximation  by  {iy”pn/9n}^i>  where  degp„  <  Q:n,degq„  <  fin.  The 
analysis  utilizes  potential  theoretic  methods.  These  are  essentially  sharp  results  though 
this  will  not  be  established  in  this  paper. 


On  A  Method  for  the  Numerical  Solution  of 
Singular  Integral  Equations  of  Cauchy  Type 

Lama  Gori,  Philip  Rabinowitz,  and  Elisabetta  Santi* 

University  of  La  Sapienza,  Rome,  Italy 
The  Weizmann  Institute  of  Science,  Rehovoth,  Israel 
University  of  L’Aquila,  Italy 

It  is  well-known  that  the  Cauchy  Singular  Integral  Equations 

(1)  a{x)v{x)  +  ^^  f  [  k{x,t)v(t)dt  =  f{x),x  e  (-1,1), 

J -it- X  J-i 


80 


axe  involved  in  several  problems  of  a  wide  variety  of  mathematical,  physical,  and  engineer¬ 
ing  areas. 

In  this  paper,  we  consider  the  problem  of  numerically  solving  singular  integral  equa¬ 
tions  of  type  (1),  with  a  and  b  constant. 

The  approach  is  based  on  a  discretization  of  the  problem,  obtained  by  the  use  of 
suitable  quadrature  rules,  related  to  weight  functions  of  Jacobi  type. 

The  existence  and  uniqueness  of  the  solution  of  the  numerical  problem  are  proved; 
these  results  are  accomplished  by  obtaining  the  expression  of  the  inverse  of  a  certain  matrix 
related  to  the  above  problem. 


A  Newton  Method  for  Polynomial  Interpolation  of  Minimal  Degree 

Thomas  Sauer 

Mathematical  Institute,  University  of  Erlangen 
Niirnberg,  Germany 

Let  Xi, . . . ,  xjv  be  a  sequence  of  points  in  R*^.  The  task  of  polynomial  interpolation  of 
minimal  degree  is  to  find  a  minimal  n  such  that  there  exists  a  subspace  P  of  the  space 
of  polynomials  in  d  variables  of  total  degree  at  most  n,  where  the  Lagrange  interpolation 
problem  is  uniquely  solvable;  i.e.,  for  any  /  :  — ♦  R  there  exists  a  imique  p  €  P  such 

that 

f{xi)=p{xi),  i  = 

A  i>articular  example  of  minimal  degree  interpolation  is  the  concept  of  least  interpolation, 
introduced  by  de  Boor  and  Ron. 

The  purpose  of  the  talk  is  to  present  a  method  to  construct  the  space  P  and  derive 
a  Newton  method  based  on  a  finite  difference  which  also  leads  to  a  remainder  formula  of 
interpolation  in  terms  of  directional  derivatives  and  simplex  splines. 


Optimal  Error  Bounds  for  Spline  Interpolation  on  a  Uniform  Partition 

F.  Dubeau  and  J.  Savoie* 

Universite  de  Sherbrooke 
Sherbrooke,  Quebec,  Canada 
and 

College  Militaire  Royal  de  Saint-Jean 
Richelain,  Quebec,  Canada 

An  n-degree  spline  s  defined  over  the  uniform  partition,  =  a-t-  Zh,  of  the  real  line 
R  of  mesh  size  /i  is  a  function  s  E  C”~^(R)  such  that  s  restricted  to  [a  +  Ih,  a  +  {I  +  l)/i] 
is  an  algebraic  polynomial  of  degree  at  most  n  for  any  I  eZ. 

To  simplify  we  use  the  notation  xi+t  =  a  +  (I  +  t)h.  For  a  function  /  defined  on  R 
and  t  €  R,  we  write  fi+t  =  fi^i+t)-  A  spline  s  is  said  to  be  the  interpolating  spline  of  /  if. 
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for  a  given  v  E  [0, 1],  we  have  sj+v  =  fi+v  for  all  I  E  Z.  Also,  e{x)  =  f{x)  —  s{x)  is  called 
the  error  function. 

Under  some  conditions  on  /  and  v,  we  show  the  two  following  inequalities. 


PJ+ul 

<  .4j(.i,t;)A"+‘-‘||/<"+'>||o. 

(1) 

and 

lie  W II 

<  4{t;)ft”+‘-‘||/<’*+‘>||oo 

(2) 

for  A:  =  0, . . . ,  n,  where 

AI^{u,v)=  f  and  A^(v)  ^  sup  A'^(u,v) 

J-oo  «€[0,1] 


are  optimal  bounds,  and  N^{u,v,‘d)  is  a  Peano  kernel. 

In  a  previous  work,  we  have  studied  inequalities  similar  to  (1)  and  (2),  but  the 
optimality  of  the  explicit  bounds  obtained  was  only  proved  for  k  =  0.  Since  A^{v)  has  a 
TniniTmim  at  V  =  Vn  where  Un  =  0  if  n  is  odd,  and  Un  =  ^  if  n  is  even,  we  only  consider 
the  case  where  v  =  Vn- 

Our  main  goal  is  to  develop  useful  techniques  to  obtain  the  numerical  values  A^{u,  v^) 
and  Using  Euler-Frobenius  polynomials,  we  obtain  new  representations  for  the  ker¬ 

nels  N^{u,  V,  t?).  Using  a  Budan-Fourier  Theorem  for  HB-splines,  it  follows  that  {u,  u„,  ■d) 
is  an  HB-spline  having  a  root  at  each  integer  (usually  simple)  and  possibly  another  one 
called  'ipiu).  When  ^(it)  does  not  exists,  we  obtain  that  all  roots  are  simple  and  deduce 
the  value  of  A^(u,Vn).  Those  values  are  related  to  Euler  splines.  When  V’Cu)  exists,  the 
problem  is  much  more  difficult.  We  obtain  results  for  n  =  2  and  n  =  3  in  that  case. 


Multivariate  Interpolation  and  Approximation 

Robert  Schaback 

University  of  Gottingen 
Gottingen,  Germany 

This  survey  talk  will  touch  the  following  topics: 

Short  introduction  into  the  theory  of  multivariate  interpolation  and  approximation 
by  finitely  many  (irregular)  translates  of  a  (not  necessarily  radial)  basis  function. 

Native  spaces  of  functions  associated  to  conditionally  positive  definite  functions,  and 
relations  between  such  spaces. 

Operators  on  radial  functions  and  construction  of  compactly  supported  positive  defi¬ 
nite  radial  functions. 

Error  bounds  and  condition  numbers  for  interpolation. 
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Uncertainty  Relation:  Why  are  good  error  bounds  always  tied  to  bad  condition 
numbers? 

Optimality  questions  for  approximation  orders:  Quasi-Optimality,  Superconvergence, 
and  Saturation. 


Approximation  on  SO(3)  and  by  Splines  and  Wavelets 

Eberhard  Schmitt 

Institut  fur  Numerische  imd  Angewandte  Mathematik 
Universitat  Gottingen 
Gottingen,  Germany 

The  compact  3-dimensional  Lie  group  SO  (3)  of  all  orthogonal  3x3  matrices  with 
determinant  +1  which  describe  the  rotations  in  appears  in  many  very  different  real 
world  applications  like  robotics  or  molecular  design.  It  can  be  parameterized  by  only  one 
chart  (C7,t)  with  C  :=  [-7r,7r]  x  [-7r,7r]  x  [-7r/2,7r/2]  using  Eulerian  angles  r  :  € 

C  — »•  SO(3).  The  coordinate  representation  f  =  F  or  :  C  — ►  R  of  a  function  F  :  S0(3) 
— >  R  on  SO(3)  has  to  fulfill  certain  obvious  periodicity  conditions  on  the  boundary  of  C. 
For  d  =  ±7r  additional  conditions  apply  due  to  the  singular  behaviour  of  r.  Furthermore 
for  differentiable  functions  F  on  SO  (3)  we  have  the  following  cheiracterization: 

Theorem.  F  is  if  and  only  if  f  fulfills  the  periodicity  conditions  and  there  are  real- 
valued  21:— periodic  functions  F±  as  well  as  2n-periodic  continuous  functions  A±,  B± 
such  that  for  — tt  <  (p  <Tr  and  — tt  ^ijj  : 

(1)  f{(p,  ±7r/2)  =  F±{(p±ip) 

(2)  ±7r/2)  =  A±{(p±  ijj)  cos  (p B±{(p  ±  ip)  sin  (p, 

(3)  (/v  =  -A±((p±ip)sm(p-\- B±{(p±ip)c(^(p. 

This  situation  is  not  very  well  suited  for  the  construction  of  splines.  But  after  trans¬ 
forming  coordinates  (p  =  x  +  y  and  ip  =  x  —  y  we  can  use  the  tensor  product 

s{x,y,i9)  =  '^CijkTi{x)Tj{y)Nk{‘d) 

ijk 

of  polynomial  (Nk)  and  trigonometric  (Tj)  R-Splines  to  calculate  an  approximation  s  to  a 
given  function  /.  The  coefficients  Cijk  then  are  obtained  as  the  solution  of  a  large  system 
of  linear  equations  consisting  of  the  peridocity  and  differentiability  constraints  and  the 
normal  equations. 

The  same  approach  can  be  used  in  a  slightly  modified  form  for  approximation  on 
S^.  In  this  case  we  will  also  consider  a  tensor  product  approach  to  the  construction  of  a 
multiresolution  setup  and  give  a  characterization  of  the  corresponding  wavelet  spaces.  By 
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introducing  the  usual  quotient  map  from  to  S0(3)  these  function  spaces  can  also  be 
used  for  approximation  on  SO  (3). 

We  will  present  an  application  of  these  methods  to  an  approximation  problem  arising 
in  molecular  modelling  and  compare  the  results  with  other  methods,  for  instance  approx¬ 
imation  by  positive  definite  functions. 


Compression  through  Rotations 
Bl.  Sendov 

Bulgarian  Academy  of  Sciences 
Sofia,  Bulgaria 

The  use  of  approximation  theory  in  the  signal  and  image  compression  explores  the 
appropriate  choice  of  the  basis  in  the  respective  functional  space.  The  most  successful  in 
this  direction  are  the  wavelet  bases,  which  are  in  some  sense  fractal  structures  and  adapt 
themselves  to  the  local  peculiarities. 

In  this  paper  we  define  a  rotation  U{X)‘,  |A|  <  1  in  the  Hilbert  space  L2  with  an  uncon¬ 
ditional  basis  $  =  Let  Hk,p  C  Ho,o  =  -^2  be  the  span  over  {<p2''t+p(i)}^o;  ^  = 

0, 1, 2, . . .  ,p  =  0, 1, 2, . . . ,  2*"  -  1.  The  rotation  U (A)  has  the  property  that  for  a  given 
f  e  H,  f{x)  =  there  exists  a  Ao  €  (-1,1),  such  that  for  the  ro¬ 

tated  basis  C/(A)$  =  {y)i(A;a;)}go  one  of  the  values  A(/;A)  =  X)»~o(‘^2i(A;/))^  and 
B{f]X)  =  E^~o(v’2^+l(A;/))^  (A(/;A)-hH(/;A)  =  \ff)  is  maximal  and  the  other  is 
Tninimal.  By  consecutive  rotations  of  the  spaces  Hp^k  we  obtain  a  representation  of  the 
function  /  which  is  appropriate  for  compression. 

This  approach  explores  the  fundamental  idea  of  the  wavelet  bases  and  the  potential 
of  the  fi:actal  transforms  for  compression. 

Computer  experiments  of  signal  and  image  compression  are  given. 


On  the  Unconditional  Convergence  of  the  Multivariate  £>’”-Splines 

Aleksei  Shadrin 

Siberian  Branch,  Russian  Academy  of  Sciences 
Novosibirsk,  Russia 


One  way  for  defining  the  multivariate  splines  is  the  variational  approach  which  leads 
to  the  £)’”-splines  of  Atteia: 

s  =  s(/,m,  A,Q)  =  argmin{  WD'^gy  :  g  G  W^{0),g\^  =  /|^}  . 

The  I>"*-splines  inherit  a  lot  of  remarkable  properties  of  the  univariate  polynomial  splines 
of  degree  2m-  1.  In  particular,  for  /  e  the  £)^-spline  interpolation  process  converges 
unconditionally  with  respect  to  the  mesh  structure.  In  Sobolev  spaces  different  from 
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W^,  convergence  of  the  £>”*-splines  can  be  provided  on  the  meshes  with  bounded  global 
mesh  ratio,  but  if  no  restrictions  are  assumed,  then  divergence  examples  (in  some  W^) 
already  exist. 

Here  we  study  the  problem  of  finding  the  necessary  and  sufficient  conditions  on  I 
and  p  providing  unconditional  convergence  of  £>'"-splines  in  Wj,.  This  problem  originated 
from  Carl  de  Boor’s  conjecture  that  unconditional  convergence  of  univariate  interpolating 
splines  of  degree  2m  —  1  takes  place  for  any  /  €  W^,  where  1  <  p  <  oo.  This  conjecture 
was  verified  for  m  =  1, 2, 3,  and  recently  we  had  shown  that  it  is  true  for  any  m  whenever 
p  is  sufficiently  close  to  2. 

We  show  that  in  contrast  to  the  univariate  case,  the  unconditional  convergence  of  the 
multivariate  i?"‘-splines  in  Wp  takes  place  if  and  only  if  p  =  2. 


Adaptive  Kernel  Estimation  of  the  Spectral  Density 
with  Boundary  Kernel  Analysis 

Alex  Sidorenko*  and  Kurt  S.  Riedel 

Courant  Institute  of  Mathematical  Sciences 
New  York,  New  York 

A  data-adaptive  non-parametric  procedure  for  estimating  the  log-spectral  density, 
d(f)  =  ln[iS'(/)],  is  proposed.  Our  scheme  has  the  following  steps: 

0)  Compute  the  Fourier  transform,  y{f)  on  a  grid  of  size  2N  -f  2  and  9i{f)  =  ln[lj/(/  •+ 
~  yif  —  ‘^)P/2(1V  +  1)]  on  a  grid  of  size  N  +  1. 

la)  Kernel  smooth  6i(/)  with  a  kernel  of  order  (0,4)  for  a  number  of  different  bandwidths, 
he,  and  evaluate  the  average  square  residual  (ASR)  as  a  function  of  he 

N 

n=l 

where  0{fn\he)  is  the  kernel  estimate  of  9{f)  using  bandwidth  he. 

lb)  Estimate  the  optimal  (0, 4)  global  halfwidth  using  a  goodness  of  fit  method.  Relate 
this  to  the  optimal  (2, 4)  using  the  halfwidth  quotient  relation. 

2)  Evaluate  the  multitaper  estimate  of  S{f)  on  a  grid  of  size  2N  +  2.  If  the  computa¬ 
tional  effort  is  not  important,  set  K  =  otherwise  choose  K  according  to  your 

computational  budget. 

3)  Estimate  6"{f)  by  smoothing  the  multitaper  estimate  with  global  halfwidth  h2,4. 

4)  Estimate  6{f)  by  substituting  0"{f)  into  the  optimal  halfwidth  expression  correspond¬ 
ing  to  the  minimum  of  (14). 

Relative  convergence  rates  are  given.  Step  2  uses  a  family  of  orthonormal  spectral 
windows  followed  by  kernel  smoothing.  The  multiple  spectral  window  procedure  pres¬ 
moothes  the  empirical  spectrum  and  reduces  the  impleasant  effects  of  the  long  tail  of  the 

log-chi  square  distribution.  We  use  the  sinusoidal  tapers:  Vn^  =  sin(^^).  These 


85 


tapers  are  an  orthonormal  family  and  their  sensitivity  to  bias  is  asymptotically  minimal. 
Multiwindowing  reduces  the  expected  mean  square  error  by  (x)  ®  over  simply  smoothing 
the  log  tapered  periodogram.  The  optimal  number  of  tapers  is 

Near  the  spectral  peaks  and  at  /  =  0  or  /  =  1/2  we  modify  the  kernel  shape  to 
take  into  account  the  boundary.  Kernels  which  minimize  the  expected  mean  square  error 
are  derived.  These  kernels  are  equivalent  to  using  a  linear  weighting  function  in  the  local 
polynomial  regression.  It  is  shown  that  any  kernel  estimator  that  satisfies  the  moment 
conditions  up  to  order  m  is  equivalent  to  a  local  polynomial  regression  of  order  m  with 
some  non-negative  weight  function  if  and  only  if  the  kernel  has  at  most  m  sign  changes. 
Fast  algorithms  are  proposed  for  computing  the  kernel  estimate  in  the  boundary  region 
and  for  evaluating  the  Rice  criterion. 

To  initialize  a  plug-in  kernel  estimator  (Step  lb),  we  use  the  fitted  average  square 
residual  criterion  (ASR)  i.e.  we  fit  the  to  the  parametric  model:  ASR{h)  =  do/ Nh  -I-  (i\h^ 
where  h  is  the  kernel  halfwidth.  This  parametric  fit  results  in  more  accurate,  less  variable 
estimates  of  the  starting  halfwidth. 


Jackson-Type  Theorems  on  the  Real  Line 
Bruno  Brosowski  and  Antonio  R.  da  Silva* 

Universitat  Frankfurt,  FVankfurt,  Germany 
and 

Federal  University  of  Rio  de  Janeiro 
Rio  de  Janeiro,  Brazil 

The  famous  theorems  of  Jackson  about  the  degree  of  approximation  by  trigonometric 
and  algebraic  polynomials  have  been  object  of  many  generalizations  and  applications. 
Typical  for  these  results  is  that  the  approximated  functions  are  bounded  and  so  they  allow 
the  construction  of  suitable  sequences  of  linear  operators  (so  called  Jackson  operators), 
which  lead  to  the  desired  estimates. 

In  our  approach  we  deal  with  functions  defined  on  the  whole  line,  not  necessarily 
bounded.  In  this  case  the  classical  approach  does  not  work  (a  result  of  Albinus).  We 
consider  then  a  convenient  metric  and  estimate  the  rate  of  convergence  of  a  approximating 
sequence,  when  structural  properties  of  the  involved  functions  are  known,  e.g.,  the  growth 
of  their  derivatives.  The  method  may  be  applied,  mutatis  mutandis,  to  entire  functions. 
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On  Fan’s  Best  Approximation  and  Applications 

S.  P.  Singh 

Memorial  University 
St.  John’s,  Newfoundland,  Canada 

Ky  Fan  proved  the  following  theorem  known  as  the  best  approximation  theorem. 

Let  C  be  a  compact  convex  subset  of  a  Banach  space  X  and  f  :C  —*  X  a.  continuous 
function.  Then  there  is  a  y  €  C  such  that 

b-fv\\  =  d{fy,C). 

=  inf  {||/j/ -  ar||  :  I  €  C} , 

This  theorem  has  applications  in  fixed  point  theory  and  variational  inequalities.  An  ex¬ 
tension  of  this  theorem  and  applications  will  be  presented  in  this  talk. 


On  the  Pointwise  Convergence  of  Fourier  Series 
with  respect  to  Periodized  Wavelets 

M.  Skopina 

St.  Petersburg  State  University 
St.  Petersburg,  Russia 

Let  3,  miilt Resolution  analysis  with  the  father-  and  mother-wavelets  re¬ 

spectively  (p  and  Ip,  satisfying  condition  (A):  ly’(r)|,  \'iP(t)\  <  J^drl),  where  uvsa  monotone 
decreasing  sununable  on  [0,-|-oo)  function.  The  tensor  product  of  m  analyses  {V^  }^_oo  ^ 
an  m-dimensional  multiresolution  analysis  in  L2(R"').  The  periodized  wavelet 

basis  of  {Vj}jL-oo  ^  orthonormal  basis  in  L2([0, 1]”*).  The  Fomrier  series  of  functions 
/  €  ^([0, 1]”*)  with  respect  to  this  basis  are  considered  (we  will  call  them  wavelet  Fomier 
series).  We  study  the  convergence  at  the  Lebesgue  points  and  at  the  strong  Lebesgue 
points,  i.e.  points  x  such  that 

f  \dt  =  0,  sup f  |/(x -f-  <)  -  f{x)\ dt  <  oo, 

I  I 

here  /  is  an  m-dimensional  interval,  symmetric  respectively  to  the  origin,  d{I)  is  the 
diameter  oil,  fil  is  the  measure  of  J.  It  is  proved  that  only  under  condition  (A)  the  wavelet 
Fomder  series  of  /  G  L([0, 1]’”)  converges  at  each  strong  Lebesgue  point.  In  particular, 
convergence  holds  almost  everywhere  for  /  G  X/logL([0, 1]”’).  If  instead  of  (A)  we  assume 
condition  (B):  |^(r)|,  |V’(t)|  ^  where  z/  is  a  monotone  decreasing  on  [0,  -l-oo)  function 

OO 

such  that  /  dr  <  oo,  then  the  wavelet  Fourier  series  of  /  G  L([0, 1]’")  converges 

0 

at  each  Lebesgue  point.  In  particular,  convergence  holds  almost  everywhere.  Condition 
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(B)  is  close  to  shaxp  because  we  have  an  example  of  multiresolution  analysis  with  father- 
wavelet  <p(t)  decaying  as  such  that  there  exists  a  function  /  €  i>([0, 1]*”)  for  which 
the  wavelet  Fourier  series  diverges  at  some  Lebesgue  point. 


Approximation  on  Compact  Sets  by  RBFs 
F.  J.  Narcowich,  P.  Smith*,  and  J.  Ward 

Texas  A&M  University,  College  Station,  Texas,  and 
WindWard  Technologies,  Houston,  Texas 

This  paper  examines  both  the  theoretical  and  practical  aspects  of  using  RBFs  as 
approximants  over  compact  sets  in  R".  In  particular,  we  point  out  the  importance  of 
choosing  the  centers  and  the  inherent  difficulties  therein.  Several  real-world  problems  are 
presented  to  emphasize  the  importance  of  choosing  the  centers  correctly. 


Multivariate  Periodic  Interpolation  and  Wavelets 
Frauke  Sprengel 

University  of  Rostock 
Rostock,  Germany 

The  aim  of  the  talk  is  to  describe  wavelet  decompositions  for  periodic  functions  in 
L2(T‘^).  Therefore  we  consider  nested  spaces  of  multivariate  periodic  functions  forming  a 
non-stationary  multiresolution  analysis. 

The  scaling  functions  of  these  spaces  are  defined  as  fundamental  interpolants  on  a 
regular  grid  or  a  sparse  grid,  respectively.  The  corresponding  wavelets  and  wavelet  spaces 
are  constructed.  Using  the  interpolatory  conditions  we  obtain  the  decomposition  and 
reconstruction  relations. 

As  an  example  we  use  certain  de  la  Vall^  Poussin  means. 

Our  investigations  generalize  the  results  for  the  univariate  case  by  J.  Prestin,  E.  Quak, 
and  K.  Selig  (1993/94)  which  were  based  on  the  paper  by  A.  A.  Privalov  (1991)  on  orthog¬ 
onal  bases  of  the  space  C2^. 


An  Algorithm  for  Minimal  Norm 
Solutions  of  Linear  Inequalities 

V.  P.  Sreedharan 
Michigan  State  University 
East  Lansing,  Michigan 

Consider  a  system  of  m  real  linear  inequalities  in  n  real  variables.  When  this  system 
is  feasible,  that  is  when  this  system  is  consistent,  it  is  natural  to  single  out  an  n- vector  x  of 
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least  norm  satisfying  these  inequalities.  In  this  paper  a  method  is  given  which  accomplishes 
this  task  for  the  problem  Ax  >  b,  ||a;||  (min),  where  the  norm  on  R”  is  the  ^-norm, 
1  <  p  <  oo.  We  show  how  this  problem  can  be  converted  to  a  non-negative  least  distance 
problem:  j  >  0,  ||c-i?u||,(min),  where  q  as  usual  is  the  index  conjugate  to  p,  viz.  p+q  —  pq. 
Algorithms  we  developed  in  some  earlier  papers  are  now  applicable  to  this  least  distance 
problem.  If  r  denotes  the  residual  of  this  problem,  then  one  can  show  that  the  case 
r  =  0  corresponds  to  the  infeasibility  of  the  system  of  inequalities.  On  the  other  hand,  if 
r  ^  0,  then  r'  is  p-dual  as  defined  in  our  earlier  papers  can  be  computed.  Prom  this  it  is 
a  small  step  to  the  solution  of  our  present  problem. 


Asymptotic  Distribution  of  Poles  and  Zeros 
of  Best  Rational  Approximants  to  |a;|“  on  [—1, 1] 

Herbert  Stahl 

TFH-Berlin/FB  2 
Berlin,  Germany 


Like  in  the  polynomial  case  so  also  in  rational  approximation  there  is  a  fundamental 
interest  in  the  approximants  of  the  functions  |a:|“,Q;  >  0,  on  [—1, 1].  Recently  the  strong 
error  formula 


.  TT 
sm  —a 
2 


^-iTy/an 


has  been  proved  for  the  minimalrational  approximation  error 


5nn  (l*r,  1-1, 1])  :=  llkl*  -  , 


where  ||||[-i,i]  denotes  the  uniform  norm  on  [—1, 1]. 

Let  r*  be  the  best  rational  approximant  of  degree  n.  It  is  shown  that  all  poles  and 
zeros  of  r*  lie  on  the  imaginary  axis,  further  poles  and  zeros  interlace.  As  a  consequence  of 
this  configuration  we  have  overconvergence  throughout  both  the  right  and  left  halfplane. 

In  the  talk  we  shall  investigate  the  asymptotic  distribution  of  these  poles  and  ze¬ 
ros.  In  addition  the  asymptotic  distribution  of  extreme  points  of  the  error  function  on 
[—1, 1]  will  be  studied.  These  questions  are  important  for  an  understanding  of  the  specific 
approximation  speed. 

Typically,  asymptotic  formulae  can  be  proved  in  two  degrees  of  precision:  The  more 
precise  formulae  allow  to  determine  the  location  of  individual  objects  (poles,  zeros,  or 
extreme  points)  as  n  — ►  oo.  The  less  precise  version  gives  a  description  only  in  terms  of 
asymptotic  densities. 

In  the  investigation  rational  interpolation,  weighted  orthogonal  polynomials,  potential 
theory,  and  some  results  from  the  theory  of  elliptical  integrals  play  a  major  role. 
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Circular  Bernstein-B4zier  Curves  and  Degree  Raising 

Sonya  Stanley 

Vanderbilt  University 
Nashville,  Tennessee 


The  theory  of  circular  Bernstein-Bezier  polynomials  was  discussed  in  a  recent  paper 
of  Alfeld,  Neamtu,  and  Schumaker.  Given  coefficients  co,  ci,  C2,  • . . ,  Cn,  the  nth  degree 
circular  Bemstein-Bezier  polynomial  is  given  by 


//>\  f  n\  /  sin 2(^2  */sin2(^  ^1)  V  «  ^  i 


The  authors  present  a  degree  raising  formula  for  circular  Bernstein-Bezier  curves,  but 
unlike  the  standard  case,  the  control  polygon  of  a  circular  Bernstein-Bezier  curve  does  not 
converge  to  the  curve  under  degree  raising.  We  show  that  with  proper  normalization,  the 
control  polygon  does  converge  to  the  curve  d&  n  00. 


On  Multivariate  Attenuation  Factors 

Gabriele  Steidl 
TH  Darmstadt 
Darmstadt,  Germany 

We  introduce  the  method  of  discrete  attenuation  factors  for  the  approximate  compu¬ 
tation  of  multivariate  discrete  Fourier  transforms. 

We  consider  attenuation  factors  related  with  multivariate  discrete  Bernoulli  functions 
and  deduce  a  best  approximation  property  of  the  corresponding  method  of  attenuation 
factors. 

Attenuation  factors  for  the  computation  of  multivariate  Fourier  coefficients  were  ex¬ 
amined  in  [3].  Choosing  a  unique  approach  to  the  discrete  and  non-discrete  settings  based 
on  the  Fourier  analysis  of  locally  compact  abelian  groups,  we  emphasize  the  close  relation 
between  both  cases  and  interpret  results  in  literature  from  a  more  general  point  of  view. 
Especially,  we  generalize  the  considerations  in  [1]  and  [2]. 

[1]  M.  G.  Ber,  W.  N.  Malozyomov,  The  best  formulae  for  the  approximate  computation 
of  discrete  Fourier  transforms,  Comput.  Math  Math.  Phys.  32  (12)  (1992),  1533—1544. 

[2]  H.  Ehlich,  Untersuchungen  zur  numerischen  Fourieranalyse,  Math.  Z.  91  (1966),  380- 
420. 

[3]  M.  H.  Gutknecht,  Attenuation  factors  in  multivariate  approximation,  Numer.  Math. 
51  (1987),  615-€29. 
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Wavelet  Analysis  of  Scattered  Data  on  the  Real  Line 

Joachim  Stodder 

University  of  Duisburg 
Duisburg,  Germany 

This  is  joint  work  with  C.  K.  Chui,  K.  Jetter,  and  J.  D.  Ward. 

We  provide  a  general  framework  for  the  construction  of  wavelets  that  are  generated 
by  translates  of  a  given  slowly  growing  function  h  where  the  translates  are  defined  by  a 
sequence  of  arbitrarily  spaced  points.  Our  approach  is  operator  theoretic  in  nature  and 
relies  only  on  specific  properties  of  the  function  h  (e.g.,  the  order  of  the  singularity  of  its 
Fourier  transform  at  the  origin)  and  not  on  a  detailed  knowledge  of  the  function  itself. 

The  methods  used  here  provide  an  exact  identification  of  the  L2-spaces  generated  by 
linear  combinations  of  shifts  of  h,  and  they  allow  the  description  of  Riesz  bases  by  certain 
coefficient  sequences  in  ^2(Z).  Furthermore,  assuming  boundedness  of  the  global  mesh 
ratio  of  the  sequence  of  points,  the  orthogonal  decomposition  of  two  embedded  L2-spaces 
is  presented,  and  Riesz  bases  for  the  corresponding  wavelet  space  are  constructed. 

Oiu:  investigation  covers  certain  classes  of  functions,  including  splines  with  simple 
knots  and  the  Hardy  multiquadric.  We  also  present  an  example  of  a  regularized  Ha.ar 
wavelet  for  scattered  data. 


Strictly  Positive  Definite  Functions 
on  the  Unit  Sphere  of  the  Complex  Hilbert  Space 

Dedicated  to  E.  Ward  Cheney 
in  Honor  of  His  65th  Birthday 

Xingping  Sun 

Southwest  Missouri  State  University 
Springfield,  Missouri 

Let  denote  the  unit  sphere  in  the  complex  Hilbert  space  H.  A  continuous  / 
defined  on  the  unit  disc  D:={z:|2|<l}is  said  to  be  positive  definite  on  S°°  if  for 
every  natural  number  n  and  any  n  points  zi,...,Zn  on  the  matrix  A  with  entries 
Aij  =  f{{zi,Zj))  is  nonnegative  definite,  where  (•,•)  denotes  the  inner  product  on  H.  If 
the  matrix  A  is  positive  definite  whenever  the  n  points  zi, . . . ,  are  distinct,  then  /  is 
said  to  be  strictly  positive  definite  on  S°°.  In  this  paper,  we  examine  several  sufficient 
and/or  necessary  conditions  that  indicate  if  a  given  function  /  is  strictly  positive  definite. 
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Malvar  Wavelets  on  Hexagons 

Xiang-Gen  Xia  and  Bruce  W.  Suter* 

Air  Force  Institute  of  Technology 
Wright-Patterson  AFB,  Ohio 

Malvar  wavelets  were  originally  developed  from  the  lapped  orthogonal  transforms  to 
eliminate  the  blocking  effects  in  transform  coding,  which  can  be  thought  of  as  a  family  of 
local  sinusoid  bases.  A  general  description  of  window  functions  for  Malvar  wavelets  was 
given  by  Coifman  and  Meyer.  With  this  description  of  window  functions,  Malvar  wavelets 
were  generalized  to  more  general  forms,  for  example,  continuous-time  case  and  discrete¬ 
time  case.  Most  recently,  a  family  of  two-dimensional  nonseparable  Malvar  wavelets  on 
rectangular  regions  were  obtained  by  Xia  and  Suter.  These  two-dimensional  nonseparable 
Malvar  wavelets  have  potential  applications  in  many  fields,  including  the  simulation  of 
turbulence  flow.  In  addition,  a  discrete  implementation  of  spatial-varying  filter  banks  was 
also  obtained. 

In  this  work,  we  want  to  construct  orthonormal  bases  or  Malvar  wavelets  defined  on  a 
hexagon  A,  using  local  bases  defined  on  its  three  rhombuses  Aj,j  =  1, 2, 3.  In  the  following, 
small  bold  English  letters,  such  as  x,y,  always  denote  two-dimensional  vectors  in  R^.  Let 
/j^/fc(x),  A:  €  Z,  be  an  orthonormal  basis  for  L^iAj)  for  each  j  =  1, 2, 3,  where  L^{A)  denotes 
all  square  integrable  functions  on  Aj .  A  trivial  way  to  construct  an  orthonormal  basis  for 
L^(A)  is  that  we  simply  use  the  truncation  windows  which  is  1  for  x  G  Aj  and 

0  otherwise,  then  fj,k{x)XAj{^)J  =  1,2, 3,  A:  €  Z,  form  an  orthonormal  basis  for  L^{A). 
However,  the  basis  elements  in  this  construction  may  have  discontinuities.  This  may  cause 
long  terms  in  signal  expansion.  Therefore,  one  desires  other  window  functions  Wj  (x)  rather 
than  the  truncations  xAj  (x)  so  that  the  windowed  local  bases  form  an  orthonormal  basis 
for  L^(A)  and  moreover,  they  are  continuous  in  A.  To  construct  such  window  functions 
and  bases  is  the  main  goal  of  this  work.  The  method  developed  in  this  work  can  be  easily 
generalized  to  Malvar  wavelets  defined  on  many  hexagons. 


A  New  Construction  of  Wavelets:  The  Lifting  Trick 

Wim  Sweldens 

University  of  South  Carolina,  Columbia,  South  Carolina 
and  Katholieke  Universiteit  Leuven 

In  this  talk  we  present  a  new  construction  of  wavelets  that  are  not  necessarily  trans¬ 
lates  or  dilates  of  one  particular  function.  The  advantage  of  such  wavelets  is  that  they 
can  be  tailored  to  specific  situations.  For  example,  one  can  construct  wavelets  that  are 
biorthogonaJ  with  respect  to  a  weighted  inner  product,  or  wavelets  that  are  defined  on 
general  subsets  of  a  Euclidean  space  (such  as  curves  or  surfaces).  Since  the  wavelets  are 
not  the  translates  and  dilates  of  one  function,  the  Fourier  transform  can  no  longer  be  used 
in  their  construction.  We  present  a  very  simple,  but  quite  powerful  alternative,  which  we 
namp  the  “lifting  trick” .  This  construction  leads  to  wavelets  that  still  have  the  powerful 
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properties  of  classical  wavelets,  such  as  biorthogonality,  compact  support,  a  fast  transform, 
vanishing  moments,  interpolation,  and  smoothness.  We  will  present  one  example,  namely 
the  construction  of  wavelets  on  a  sphere,  in  greater  detail. 


Interpolating  Wavelets  on  [—1, 1]  and  the  Sphere 
Manfred  Tasche 

University  of  Rostock 
Rostock,  Germany 

The  construction  of  wavelets  on  an  interval  was  recently  considered  by  Meyer  1992, 
Chui-Quak  1992,  Cohen-Daubechies-ViaJ  1993  and  some  others.  In  the  following  we 
describe  a  new  approach  to  polynomial  wavelets  on  [—1,  1]  introduced  by  Kilgore-Prestin 
1994.  Here  we  apply  Chebyshev  transforms,  related  shifts  (see  Butzer-Stens  1977)  and 
discrete  cosine  transforms.  It  is  essential  that  we  do  not  use  interior  and  edge  scaling 
functions  and  wavelets. 

Our  sampling  space  Vj  (J  €  No)  is  spanned  by  Chebyshev-shifts  of  a  Lagrange  function 
(fj  with  respect  to  the  Chebyshev  nodes  cos  {2~^kTr)  (fc  =  0, . . . ,  2^).  Note  that  Vj  coincides 
with  the  set  of  all  polynomials  of  degree  <  2^.  The  chain  of  all  {j  E  No)  forms 
a  multiresolution  of  the  weighted  Hilbert  space  L^[— 1,  1]  (with  w{x)  :=  (1  — 

(x  G  (—1,  1)).  The  wavelet  space  Wj  (j  G  No)  is  the  orthogonal  complement  of  Vj  in  V)+i 
and  is  spanned  by  Chebyshev-shifts  of  the  wavelet  V’j+i*  The  shift-invariant  spaces  Vj  and 
Wj  can  be  easily  characterized  by  Chebyshev  transforms  and  discrete  cosine  transforms. 

Using  tensor  products  of  interpolating  trigonometric  wavelets  and  polynomial  wavelets 
introduced  above,  we  obtain  new  wavelets  on  the  sphere.  The  scaling  functions  of  level  j 
(j  G  No)  are  chosen  as  Lagrange  fundamental  polynomials  with  respect  to  grid  {2~^k‘K, 
cos(2~'^77r))  (A:,  /  =  0, . . . ,  2^).  For  the  decomposition  and  reconstruction  of  a  given  func¬ 
tion,  we  describe  very  efficient  and  numerically  stable  algorithms  based  on  fast  Fourier 
transform  and  discrete  cosine  transform.  Numerical  tests  will  be  presented. 


Complex  Rational  Uniform  Approximation  Applied  to  the  Design  of 
Iterative  Methods  for  Linear  Systems  of  Equations 

Jean-Pierre  Thiran  and  Suzanne  Thiry* 

Facultes  Universitaires  N.D.  de  la  Paix 
Namur,  Belgium 

For  solving  the  real  nonsingular  system  of  linear  equations  x  =  Tx  +  c,  we  consider 
stationary  (A:,  /)-step  methods  defined  by  the  iterative  scheme 

i  k 

PO^m  ~  T(^  i9t)  "I"  C  ^  ]  ^m—jPj 

i=l  i=l 
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where  po  ^  0,  ^JLq  Pi  =  E!=i  9t  =  1- 

One  way  of  determining  the  parameters  pi  and  Qj,  is  to  find  a  rational  approximation 
of  the  conformal  mapping  from  the  exterior  of  the  unit  disk  onto  the  exterior  of  a  region 
that  encloses  the  spectrum  of  T  [1].  In  the  case  of  a  polygonal  domain  fi,  Li  [2]  employed 
a  least-square  criterion  for  approximating  the  Schwarz-Christoffel  mapping  computed  with 
the  aid  of  SCPACK  [3]. 

The  present  contribution  is  intended  to  provide  some  numerical  comparison  based  on 
the  use  of  imiform  approximation  procedures. 

References 

[1]  Gutknecht,  M.  H.,  Stationary  and  almost  stationary  iterative  (A:,/)-step  methods  for 
linear  and  nonlinear  ^sterns  of  equations,  Numer.Math.  56  (1989),  179-213. 

[2]  Li  X.,  An  adaptive  method  for  solving  nonsymmetric  linear  systems  involving  appli¬ 
cations  of  SCPACK,  J.  Comp.  App.  Math.  44  (1992),  351-370. 

[3]  TVefethen  L.  N.,  Numerical  computation  of  the  Schwarz-Christoffel  transformation, 
Siam.  J.  Sci.  Slot.  Comput.  1  (1)  (1980),  82-102. 


Quaternions  on  Wavelet  Problems 
Leonardo  Traversoni 

Universidad  Autonoma  Metropolitana,  Iztapalapa 
Mexico  D.F. 

There  are  several  problems  and  applications  where  wavelets  are  used  where  inclusion 
of  Quaternions  could  be  of  great  help  simplifying  the  formulation  of  the  problem. 

This  paper  shows  first  how  Quaternion  Wavelets  can  be  built  and  some  of  its  properties 
and  later  some  applications  first  in  its  classical  use  of  representing  rotations  and  afterwards 
in  not  so  classical  applications  dealing  with  high  dimensions  and  providing  a  tool  for 
additional  steps  of  compression  of  information. 

One  of  the  big  advantages  of  using  Quaternions  is  their  easy  algebraic  representation 
that  fits  specially  for  computational  uses. 


Data  Fitting  by  Integer  Quasi-Convex  Functions 

Vasant  A.  Ubhaya 

North  Dakota  State  University 
Fargo,  North  Dakota 

Given  n  data  points  /  =  (/i,/2, •  •  •  ,/n),  the  problem  is  to  find  a  best  quasi-convex 
fit  g  =  ,9n)  with  integer  gi  which  minimizes  a  suitable  norm  \\f  -  gW-  The 

norms  under  consideration  axe  least  squares,  £i  and  uniform.  Algorithms  of  polynomial 
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complexity  (in  n)  axe  developed  for  computation  of  best  fits.  These  problems  arise  in 
data-fitting  and  statistical  estimation. 


Divergence-FVee  Wavelet— Bases 

Kaxsten  Urban 

RWTH  Aachen 
Aachen,  Germany 

In  the  last  few  years  the  existence  of  divergence-free  wavelets,  that  was  shown  by  Bat¬ 
tle  and  Federbush,  was  used  for  the  proof  of  existence-theorems  of  classical  solutions  of  the 
incompressible  Navier-Stokes-Equations.  For  numerical  purposes  these  vector  fields  have 
some  disadvantages:  they  axe  in  general  not  compactly  supported  and  the  computation 
of  these  wavelets  seems  to  be  difficult.  In  1992  Lemaxie-Rieusset  constructed  compactly 
supported  divergence-free  wavelets,  whose  computation  was  standard.  His  construction 
was  based  on  a  tensor-product  approach,  that  seems  not  be  helpful  for  applications  in  fluid 
dynamics. 

In  this  talk  we  show  the  basic  ideas  of  a  new  construction  of  vector-wavelets  with 
vanishing  divergence  without  using  tensor  products.  It  will  be  pointed  out,  that  these 
wavelets  give  rise  to  a  stable  basis  of  the  space  of  all  divergence-free  vector  fields  and  that 
they  axe  not  hard  to  compute.  We  give  some  examples  based  on  box-splines  and  show 
applications  in  the  numerical  treatment  of  incompressible  Navier-Stokes-Equations. 


Step-by-Step  Algorithms  in  the  Multidimensional  Spline  Approximation 

Vladimir  A.  Vasilenko 

Siberian  Branch,  Russian  Academy  of  Sciences 
Novosibirsk,  Russia 

The  multidimensional  spline  approximation,  especially  at  the  huge  scattered  mesh  in¬ 
stead  of  regular  grid,  is  the  expensive  computational  procedure  both  in  analytical  and  in 
finite  element  methods.  In  this  communication  we  present  an  effective  method  for  the  solv¬ 
ing  of  system  for  the  variational  FE  analog  of  spline  based  on  special  tensor  decomposition 
of  multi-index  matrices  and  vectors.  It  leads  to  non-expensive  step-by-step  improvement 
algorithm  to  provide  an  a  priori  given  error  bound.  Simultaneously,  the  fluent  linear  in¬ 
creasing  of  data  volume  is  provided  in  the  steps,  and  effective  data  compression  is  achieved 
for  the  low  accuracy  requirements. 
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Construction  of  Fundamental  Interpolants  and  Wavelets 

Johan  M.  de  Villiers 

University  of  Stellenbosch 
Stellenbosch,  South  Africa 

This  is  joint  work  with  C.  K.  Chui.  Our  objective  is  to  introduce  a  general  scheme 
for  the  construction  of  interpolatory  approximation  formulas  and  compactly  supported 
wavelets  by  using  spline  functions  with  arbitrary  (nonuniform)  knots.  Both  constructions 
use  “optimally  local”  interpolatory  fundamental  spline  functions  which  are  not  required  to 
possess  any  approximation  property. 


Representing  the  Error  in  Multivariate  Polynomial  Interpolation 

Shayne  Waldron 

University  of  Wisconsin-Madison 
Madison,  Wisconsin 


I  give  a  quick  tour  of  my  recent  results  concerning  constructive  instances  of  the  er¬ 
ror  bounds  that  numerical  analysts  use  to  conclude  that  a  scheme  based  on  polynomial 
interpolation  has  the  highest  possible  order  of  accuracy  that  its  polynomial  reproduction 
allows. 

Highlights  in  the  univariate  case  include  the  bound 

|/(x)  -  Hemi  < 

where  H  —  ©/  is  the  Hermite  interpolant  to  /  at  the  multiset  0  of  n  points,  (j{x)  := 
-  9),  and  diam  (resp.  conv)  is  the  diameter  (resp.  convex  hull)  of  a  (multi)set 
of  points.  This  inequality  significantly  improves  upon  Beesack’s  inequality,  upon  which 
almost  all  the  error  bounds  for  Hermite  interpolation  given  over  the  last  30  years  have 
been  based.  Some  interesting  pictures  of  the  Birkhoff  splines  which  describe  the  error  in 
Hermite  interpolation  are  presented. 

Highlights  in  the  multivariate  case  include  the  multivariate  form  of  Hardy’s  inequality, 
that  for  m  —  n/p  >  0 


/  llLp(n)  < 


r(m  -  n/p)  ,, 
r(m)r(#©  +  m  -  n/p) 


V/  €  Lp(f2), 


where  ©  is  a  finite  sequence  of  points  in  M”,  and  /  Jq  /  is  (up  to  a  scalar  multiple)  the 
linear  fimctional  of  integration  against  the  simplex  spline  M(-|©)  (introduced  by  Micchelli, 
1980).  This  inequality  has  great  utility  and  can  be  applied  to  the  many  integral  error 
formulae  that  exist,  including  those  for  Kergin/Hakopian  interpolation,  and  recent  work  of 
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Sauer  and  Xu.  There  is  some  discussion  on  formulae  for  the  error  in  linear  interpolation  on 
a  triangle,  including  presentation  of  a  formula  which  exhibits  the  geometry  of  the  problem 
in  a  pleasing  way. 


Construction  of  Orthonormal  Wavelets  with  Multiplicity  r  with 
the  Same  Support  to  Increase  the  Number  of  Vanishing  Moments  r-fold 

Jianrong  Wang 

HARC,  The  Woodlands,  Texas 

This  is  a  joint  work  with  C.  K.  Chui.  Let  r  >  2  be  any  preassigned  integer.  To  any 
compactly  supported  orthonormal  scaling  function  whose  integer-translates  contain  locally 
all  polynomials  of  order  m  (or  degree  m  —  1)  in  their  linear  span,  we  construct  r  —  1  addi¬ 
tional  scaling  functions  with  the  same  support,  so  that  the  integer  translates  of  all  of  the 
r  scaling  functions  constitute  an  orthonormal  family,  and  that  together,  they  provide  local 
approximation  of  order  rm.  In  addition,  we  construct  the  r  corresponding  orthonormal 
wavelets  with  the  same  support.  Consequently,  these  wavelets  have  rm  (instead  of  just  m) 
vanishing  moments. 


Cubic  Spline- Wavelet  Packets  of  Sobolev  Spaces 
and  Adaptive  Approximation 

Jianzhong  Wang 

Sam  Houston  State  University 
Himtsville,  Texas 

Let  Hq{I)  denote  homogenuous  Sobolev  space  : 

=  {m  e  m  =  /'(O)  =  f(L)  =  /'(£)  =  0}, 

with  the  inner  product  <  f,g  >=  fjf"{x)g"(x)dx,  where  I  =  [0,  L].  (For  the  sake  of 
simplicity,  we  assume  that  L  is  an  integer  greater  than  3.) 

Let  Vj  be  the  subspace  of  all  cubic  splines  in  Hq{I)  with  knots  {^}.  It  is  well- 
known  that  forms  an  MRA  for  Hq(I)  equipped  with  norm  |||/|||  =<  /,/  >= 

fj  f"{x)g''{x)  dx.  Let  fj  be  the  cubic  spline  in  Vj  interpolating  /.  Then  fj  is  the  optimal 
approximation  of  /  in  V  under  the  norm  III  •  |||.  The  function  /,-  is  uniquely  determined 

by  the  vector  f*  :=  {/(^)}fcif  ^ 

We  construct  two  compactly  supported  cubic  splines  in  Vi,  ‘ip{x)  and  tpbix)  with 
suppif)  =  [0,3]  and  suppx()b  =  [0,2].  The  j-level  dilation  and  translations  of  ip  and  j-level 
dilation  and  reflections  of  ipb  form  the  basis  of  Wj,  which  is  the  orthogonal  complement  of 
Vj  in  Vj^i.  Then  Hq{I)  =  Vq  Wj.  By  the  way,  Vq  and  Wj,0  <  j  <  oo,  also  form  a 
decomposition  of  the  space  Cq.  Notice  that  any  function  gj  in  Wj  is  uniquely  determined 
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by  :=  {93(20^)}%^  ^j(^)  =  0,  Wk.  Because  of  this  property,  this  kind  of  basis 

has  two  advantages.  (1)  It  can  be  used  to  decompose  sampling  data  directly.  (2)  Distinct 
from  the  usual  wavelet  basis  in  L^,  it  provides  an  algorithm  which  is  easy  to  perform  by 
using  parallel  algorithm.  That  means,  if  the  sampling  data  fj  :=  are  given, 

we  can  decompose  it  into  all  levels  in  parallel. 

Now  we  write  gj{x)  —  The  following  theorem  tells  us  that  the  coeffi¬ 

cients  of  the  decomposition  can  reflect  the  local  singularity  of  the  fimction. 

Theorem  1.  Let  0  <  a  <  4.  Suppose  the  function  /  is  Holder  continuous  with  exponent 
a  at  xo  G  I.  Then  for  any  k  E  Z,  j  £  Z"*"  such  that  2~^k  E  {xq  —  6f2,xo  +  <5/2), 

dj,k  =  0(2-“^),  j  -  00. 

Hence,  for  big  j,  only  the  coefiicients  corresponding  to  the  singularities  are  relatively 
large.  Then  the  terms  with  small  coefficients  can  be  deleted  without  causing  an  error 
beyond  tolerance.  The  following  theorem  tells  us  that  these  terms  can  be  essentially 
determined  by  the  corresponding  sampling  values  of  gj.  Hence  we  can  predetermine  the 
terms  to  be  deleted  . 


Theorem  2.  Let  M  =  maxi\f{x)\  and  gj(x)  be  its  wavelet  component  in  Wj.  Suppose 
that  for  c  >  0,  — 1  <  fci  <  A:2  <  Uj  —  2, 


(2k -I 
9j\  2j+i  ' 


<  e, 


for  ki  <  k  <  k^. 


Then  the  function 

-I<fc<nj-2,fc3(fci+/,fc2-l] 

satisfies 

where  I  =  l{e)  =  [0.87  -f  log  j]. 

If  the  function  /  has  only  isolated  singularities,  the  quantity  of  nonvanishing  terms 
in  each  j-level  are  almost  the  same  while  the  dimension  of  Wj  is  exponentially  increasing. 
This  fact  is  crucial  for  adaptive  approximation,  data  compression,  numerically  solving  PDE 
and  other  applications. 

The  coded  data  by  using  wave  decomposition  described  above  have  poor  frequency 
resolution  when  the  level  is  high.  In  order  to  get  finer  frequency  resolution  in  high  level,  we 
construct  the  spline- wavelet  packet  in  Sobolev  space  Hq{I)  as  following.  Let  di,d2,  -  •  •  ,dn 
be  a  sequence  of  Os  and  Is,  ^  =  1  —  dn,  and  q:„  =  Define  C  Ws 

forn  >  0  such  that  any  E  satisfies 


dTd-2  -d„(  \  _  Q 

s  V2n+s-|-l'' 


kEZ\{2^+^Z  +  an}. 


98 


(We  denote  forn  =  0.)  Then  •"‘in+i  ^  --in+i 

and  Vj  =  V^_a  ®iZi  ©(dida  -dj)^/-**  ***•  Each  ■‘^’*,  n  =  1,2,  •••,  has  the  same 

dimension  2*L. 

Using  the  tensor  product  we  can  obtain  the  multidimensional  spline-wavelet  packet. 
For  any  function  decomposed  by  such  a  packet,  the  nonvanishing  coefficients  of  its  com¬ 
ponents  in  the  distinct  subspace  <S>  have  almost  the  same  time  (or 

space)  locations  in  their  corresponding  coefficient  matrices.  The  procedure  to  obtain  the 
coefficients  for  the  spline-wavelet  packet  basis  is  simpler  than  that  for  wavelet  basis.  It 
is  particularly  useful  for  multi-dimensional  problems  such  as  multivariable  approximation, 
image  processing  and  nonlinear  PDE  etc. 
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Wavelets  Associated  with  Periodic  Basis  Functions 

Francis  J.  Narcowich  and  Joseph  D.  Ward* 

Department  of  Mathematics 
College  Station,  Texas 

In  this  talk,  we  investigate  a  class  of  nonstationary,  orthogonal,  periodic  scaling  func¬ 
tions  and  wavelets  generated  by  continuously  differentiable  periodic  functions  with  positive 
Fourier  coefficients;  these  functions,  which  were  introduced  long  ago  by  Schoenberg,  will 
be  called  periodic  basis  functions  (PBFs). 

The  scaling  fimctions  and  wavelets  presented  here  have  a  number  of  attractive  features. 
Their  translates  are  orthonormal,  thereby  eliminating  the  need  for  dual  functionals.  Such 
wavelets  also  show  greater  stability,  because,  for  orthonormal  bases,  the  Riesz  constants  are 
unity.  The  algorithms  for  analyzing  with  them  are  ultimately  simple.  The  decomposition 
and  reconstruction  coefficients  can  be  computed  in  terms  of  the  discrete  Fourier  transform, 
so  that  FFT  methods  apply  for  their  evaluation.  In  addition,  decomposition  at  the 
level  only  involves  2  terms  from  the  higher  level,  as  opposed  to  the  2n  terms  one  needed 
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with  the  earlier  periodic  wavelets.  A  similar  result  applies  for  reconstruction.  Finally, 
for  many  PDFs,  the  corresponding  scaling  functions  and  wavelets  have  good  localization 
properties.  One  of  the  main  results  in  this  paper  is  that  both  the  scaling  functions  and 
wavelets  derived  from  a  large  class  of  PDFs  have  “periodic”  standard  deviations  that  are 
0{l/n)  at  level  n,  where  the  notion  of  periodic  standard  deviation  is  that  described  by 
Breitenberger.  In  fact,  a  periodic  “uncertainty  principle”  applies  in  this  setting,  and  one 
can  show  that  the  angle/frequency  uncertainty  “window”  is  0{y/n)  at  level  n,  at  least  for 
scaling  functions  wavelets  associated  with  a  certain  class  of  PBFs. 


De-noising  Using  Wavelets  and  Cross  Validation 

N.  Weyrich  and  G.  T.  Warhola* 

Air  Force  Institute  of  Technology 
Wright-Patterson  AFB,  Ohio 

We  present  numerical  results  for  noise  reduction  performed  by  modified  wavelet  re¬ 
construction  and  address  the  automatic  choice  of  the  related  threshold  and  shrinkage  pa¬ 
rameter  without  any  prior  knowledge  about  the  noise  variance.  We  show  that  the  cross 
validation  method  can  be  a  helpful  tool  for  making  this  choice. 

As  an  application,  we  consider  the  problem  of  removing  noise  from  speech.  Our 
model  is  s{ti)  —  x{ti)  +  n{ti)  for  uniform  time  samples,  U,  where  x{ti)  represents  the  clean 
speech  to  be  estimated  and  the  n{ti)  are  independent,  identically  distributed  noise  samples 
having  zero  mean  and,  possibly  unknown,  variance,  (t^.  We  investigate  both  thresholding 
and  shrinkage  of  wavelet  coefficients  to  perform  this  task  with  both  prior  and  no  prior 
knowledge  of  cr^.  Wavelet  coefficients  of  s  are  nonlinearly  driven  toward  zero  based  on 
their  relation  to  a  threshold.  The  estimate,  x,  reconstructed  from  these  modified  wavelet 
coefficients  is  supposed  to  approximate  the  underlying  speech,  x.  Our  investigations  include 
the  use  of  wavelets  on  the  whole  real  line,  periodized  wavelets,  and  wavelets  constructed 
on  an  interval,  the  latter  two  types  being  applied  to  windowed  speech  segments.  We  also 
consider  orthogonal  (spline  and  Daubechies’  wavelets),  semi-orthogonal  (spline)  wavelets, 
orthogonal  trigonometric  wavelets  and  wavelet  packets. 

With  prior  knowledge  of  we  choose  a  threshold  parameter  using  the  method  of 
Donoho,  based  on  Stein’s  unbiased  estimate  of  risk. 

When  the  variance  is  unknown,  we  estimate  the  optimal  choice  of  threshold  using 
methods  of  cross  validation,  ordinary  and  generalized.  We  demonstrate  that  this  estimation 
is  close  to  the  ideal  threshold  parameter  for  test  cases  involving  knowledge  of  the  true 
speech  signal. 
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Smooth  Functional  Subdivision  Methods  for 
Irregular  Triangulations  of  the  Plane 

Joe  Warren 

nice  University 
Houston,  Texas 

For  a  wide  class  of  subdivision  methods  over  irregular  triangulations,  we  derive  nec¬ 
essary  conditions  for  these  schemes  to  produce  continuous  limit  surfaces.  We  next 
give  sufficient  conditions  for  these  schemes  to  produce  continuous  limit  surfaces.  We 
conclude  with  an  example  of  such  an  irregular  subdivision  method. 


An  Approximation  Problem  in  the  Analysis  of  Measurements 

G.  A.  Watson 

University  of  Dundee 
Dundee,  Scotland 

At  the  1992  Texas  meeting,  Zwick  introduced  a  planar  minimax  algorithm  for  the 
analysis  of  measurements  related  to  the  acceptance  or  rejection  of  a  manufactured  object. 
The  corresponding  least  squares  problem  can  be  solved  as  an  orthogonal  Procrustes  prob¬ 
lem,  and  the  minimax  problem  can  be  interpreted  as  a  weighted  least  squares  problem. 
Algorithms  based  on  this  are  considered,  as  well  as  related  methods  which  apply  to  various 
extensions  of  the  underlying  mathematical  problem. 


Non-Stationary  Multiresolution  Analysis  and 
Corresponding  Biorthogonal  Systems 

Ilona  Weinreich 

RWTH  Aachen 
Aachen,  Germany 

A  generalization  of  the  concept  of  Multiresolution  Analysis  is  obtained  by  dropping 
the  assumption  that  all  spaces  are  generated  by  the  dilates  of  a  single  function.  A  typical 
example  arises  in  connection  with  exponential  splines.  Our  objective  is  to  construct  a 
biorthogonal  wavelet  basis  corresponding  to  exponential  splines  and  to  indicate  possible 
applications.  A  central  issue  in  this  context  is  the  efficient  computation  of  the  dual  mask 
coefficients. 
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Fast  Algorithms  for  Trigonometric  Wavelet  Packets 
and  Applications  in  Speech  Processing 

Ewald  Quak  and  Norman  Weyrich* 

Air  Force  Institute  of  Technology 
Wright-Patterson  Air  Force  Base,  Ohio 

The  aim  of  this  talk  is  to  describe  wavelet  packet  functions  and  spaces  for  a  trigono¬ 
metric  multiresolution  analysis  based  on  fundamental  Lagrange  interpolants.  The  corre¬ 
sponding  algorithms  for  wavelet  packet  decomposition  and  reconstruction  are  investigated 
in  detail.  As  an  example,  an  application  of  trigonometric  wavelet  packets  to  speech  pro¬ 
cessing  is  considered. 


Characterization  Theorem  of  Best  Lp  (1  <  p  <  -f  oo) 
Approximation  with  Restricted  Ranges  of  Derivatives 

Shusheng  Xu 

Jiangnan  University 
Wuxi,  Jiangsu,  PRC 

We  assume  that  =  span{(^i, . . . , 0n}  is  a  n-dimensional  subspace  of  Lp[a, b], 
. . ,  (j>^^  exist  (s  =  0,1,...,  k)  and  have  a  maximal  linearly  independent  subset 
which  is  an  extended  Chebyshev  system.  Write 

=  {?  e  :  Isix)  <  q^"Hx)  <  u^ix),x  G  [a, 6]}  , 

where  Ig  and  Ug  are  extended  real  valued  functions  satisfying  -oo  <  lg{x)  <  Ug{x)  < 
-1-00.  In  this  paper  we  give  a  characterization  for  uniform  approximation  to  /  G  C[a,b] 
by  Ks  :  The  result  is  widely  applicable  and  containing  many  special  cases 

such  as  approximation  with  HB  interpolatory  restrictions,  comonotone  approximation, 
and  approximation  by  polynomials  with  bounded  coefRcients,  etc.  In  case  of  1  <  p  <  -l-oo, 
we  call  q  =  a.j<j>j  a  vector,  define  Euclidean  inner  product  in  $n,  and  write  K\  = 

{q  E  :  {q,hx)  <  dx,X  G  A},  where  dx  is  a.  real  number  and  hx  a  nonvanishing  vector 
in  for  each  A  G  A.  We  say  Ka  is  a  “local  convex  cone”  at  qo  G  Ka  if  there  exists 

a  ^  >  0  such  that  {g  G  :  11?  —  9o||  <  5}  C  {g  G  :  (9>^a)  ^  dx,X  e  A\A'}  with 

A'  :=  {A  G  A  :  (qo,  hx)  =  dA}-  Then  we  get  a  characterization  theorem  of  Lp  approximating 
from  the  product  of  Ks  and  a  “local  convex  cone”.  This  result  is  widely  applicable  as  well. 
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Computing  Continuous  Wavelet  Transforms  By  Subdivision 

Shuzhan  Xu 

University  of  Alberta 
Edmonton,  Canada 

In  real  applications,  efficient  algorithms  are  needed  to  compute  wavelet  transforms. 
For  the  continuous  wavelet  transform,  Holschneider,  Kronland-Martinet,  Morlet  and 
Tchamitchian  proposed  the  efficient  algorithme  d  trous  which  is  based  on  a  hole  filling 
process.  The  hole  filling  process  is  limited  in  the  sense  that  it  uses  spline  interpolation 
that  is  restricted  to  piecewise  constant  and  piecewise  linear  splines.  In  search  of  an  efficient 
algorithm  with  higher  order  accuracy,  we  extend  the  algorithme  d  trous  by  just  approxi¬ 
mating  the  analyzing  wavelet  by  a  linear  combination  of  refinable  functions;  approximating 
by  quasi-interpolation  in  the  spline  case.  The  refinability  permits  subdivision  schemes  to 
be  applied  in  the  computation  which  lead  to  a  very  efficient  new  algorithm.  The  algorithm 
has  an  almost  identical  implementation  as  the  algorithme  d  trous  but  the  starting  point 
is  different.  Our  algorithm  coincides  with  the  algorithme  d  trous  in  the  one  dimensional 
case  if  the  initial  approximation  is  given  by  piecewise  constant  or  piecewise  linear  spUne 
interpolation.  This  coincidence  reveals  the  connection  between  subdivision  and  hole  fill¬ 
ing,  which  are  two  basic  approaches  to  render  given  data  and  provide  approximation.  The 
algorithm  is  an  extension  of  algorithme  d  trous  for  it  provides  higher  order  approximation 
and  easily  generalizes  to  higher  dimensions. 


Fej6r  Means  for  Multiple  Fourier  Series 
Yuan  Xu 

University  of  Oregon 
Eugene,  Oregon 

This  is  joint  work  with  Hubert  Berens.  The  following  result  is  an  analogue  of  Fejer’s 
theorem  on  Fejer  kernel: 

Theorem.  Let  f  €  C{T*^)  and  let  denote  the  n-tb  partial  sum  of  its  Fourier 

Series,  where  the  sum  is  taken  w.  r.  t.  all  multi-indices  a  E  with  modulus  |ali  := 
|ai|  -f  . . .  -4-  lodi  <n,nE  Nq.  Then  the  Cesaro  (C,  2d  —  1)  means  of  3^]^  define  a  positive 
approximate  identity;  the  order  of  summability  is  best  possible. 

The  proof  of  the  theorem  depends  heavily  on  results  firom  the  theory  of  special  func¬ 
tions.  We  also  discuss  the  analogous  results  on  Fourier  integrals. 
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Construction  of  Wavelets  on  Finite  Domains 

Charles  A.  Micchelli  and  Yuesheng  Xu* 

IBM  T.  J.  Watson  Research  Center 
Yorktown  Heights,  New  York  10598 

and 

University  of  North  Dakota 
Fargo,  North  Dakota 

We  will  present  the  construction  of  wavelets  on  finite  domains  and  their  applications. 


On  Copositive  Polynomial  Approximation  in  Lp[-l,l] 

Y.  K.  Hu,  K.  A.  Kopotun,  and  X.  M.  Yu* 

Southwest  Missouri  State  University 
Springfield,  Missouri 

For  function  /  e  Lp[-1, 1],  0  <  p  <  oo,  with  finitely  many  sign  changes,  we  construct 
a  sequence  of  polynomials  P„  €  n„  which  are  copositive  with  f  and  such  that 

ll/-/Vllp<Ca,*(/,(ri+ir')^, 

where  denotes  the  Ditzian-Totik  modulus  of  continuity  in  Lp  metric. 


The  Integral  Representation  of  Bivariate  Splines 

Guanglu  Zhang 

University  of  Petroleum 
Dongying,  PRC 

This  is  joint  work  with  H.  W.  Liu.  In  this  paper  we  give  an  integral  representation 
of  bivariate  splines.  As  an  application,  we  use  such  representation  of  splines  to  study  the 
dimension  problem  of  the  quadratic  spline  space  <52  (A)  over  some  triangulation  A. 


Approximation  in  Norm  and  in  Sobolev  Seminorm 

Kang  Zhao 

University  of  Utah 
Salt  Lake  City,  Utah 

For  any  (f)  in  the  shift-invariant  space  S{(f>)  generated  by  ^  provides  simulta¬ 

neous  approximation  order  (m.  A:),  with  k  >  m,  iff  it  provides  approximation  order  k  —  fnin 
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the  Sobolev  seminorm  |  •  \m,2‘  If  each  S(D°‘<f>)  provides  approximation  order  k  for  |a|  =  m 
then  provides  approximation  order  k  +  m. 


On  Signal  Processing  by  Means  of  Periodic  Spline  Wavelets 

Valery  A.  Zheludev 

Military  Institute  for  Construction  Engineering 
St.  Petersbiug,  Russia 


The  objective  of  the  presentation  is  an  outline  of  techniques  for  adaptive  signal  pro¬ 
cessing  (SP)  based  on  the  spline  wavelet  analysis. 

Early  examples  of  wavelets  were  based  on  spline  functions.  Later,  spline  wavelets  were 
shadowed  by  the  compactly  supported  orthonormal  wavelets  by  I.  Daubechies.  However, 
in  a  good  deal  of  situations  spline  wavelets  offer  advantages  before  these  latter.  Especially 
these  manifest  itself  in  SP,  where  the  property  of  linear  phase  is  required,  which  fails  when 
employing  the  wavelets  by  I.  Daubechies.  In  recent  years  cardinal  spline-wavelets  were 
subjected  to  detailed  studying,  mainly  by  C.  K.  Chui  with  collaborators.  In  particular,  the 
authors  succeeded  in  constructing  compactly  supported  wavelets  in  the  spaces  of  cardinal 
splines.  The  remarkable  property  of  such  spline  wavelets,  in  addition  to  the  symmetry,  is 
that  their  dual  wavelets  belong  to  the  same  spaces  as  the  original  ones. 

We  discuss  here  periodic  spline  wavelets.  Our  approach  to  the  spline  wavelet  analysis 
is  based  on  an  original  computational  technique — ^the  so-called  Spline  Harmonic  Analysis 
(SHA)  which  is  a  version  of  the  Harmonic  Analysis  (HA)  in  spline  spaces.  SHA  in  some 
sense  bridges  the  gap  between  the  continuous  and  the  discrete  versions  of  HA.  It  is  a  rather 
universal  technique  applicable  to  a  great  variety  of  numerical  problems,  not  necessarily  to 
wavelet  analysis. 

The  employment  of  SHA  techniques  in  wavelet  analysis  is  found  to  be  remarkable 
fruitful.  This  approach  has  given  a  chance  to  construct  a  rich  diversity  of  spline  wavelets 
including  periodizations  of  the  above  mentioned  compactly  supported  spline  wavelets  and 
their  duals  as  well  as  the  selfdual  wavelets  by  Battle  and  Lemarie.  There  constructed 
the  so-called  high-  and  low-frequency  mother  wavelets  including  locally  supported  (up  to 
periodization)  ones. 

These  mother  wavelets  generate  a  library  of  semi-orthogonal  wavelet  packet  bases  of 
spline  spaces.  Moreover,  the  SHA  approach  provides  an  efficient  algorithm  for  changing 
from  the  H-spline  basis  to  a  wavelet  and  wavelet  packet  basis  of  the  spline  space  and  back. 

These  latter,  together  with  some  quadrature  formulas  established,  enable  us  to  suggest 
adaptive  fast  methods  for  digital  processing  of  one-  and  multi-dimensional  signals  by  means 
of  spline  wavelets  and  wavelet  packets.  We  emphasize  that  the  relations,  being  derived  by 
means  of  SHA  methods  for  periodic  signals,  provide  algorithms  for  local  processing  non¬ 
periodic  signals  in  real  time  conditions  as  well. 

The  algorithms  suggested  allow  simple  fast  numerical  implementation.  There  devel¬ 
oped  a  software  for  SP  on  the  basis  of  these  ones. 
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A  Good  Interpolation  Basis,  the  Marcinkiewicz-Zygmund  Inequality 

and  Ap  Weights 

Lefan  Zhong 

Peking  University 
Peking,  China 

This  is  joint  work  with  C.  K.  Chui.  The  equivalence  of  a  good  interpolation  basis 
and  the  Marcinkiewicz-Zygmund  inequality  is  established.  It  follows  from  the  theory  of 
singular  integral  operators  and  Carleson’s  interpolation  theorem  that  they  are  equivalent 
to  the  separation  of  sample  points  and  the  Ap  weight  condition  on  the  corresponding 
fundamental  polynomial. 


Asymptotic  Expansions  of  the  Representation  Formulae 
for  (C'o)m  Paremeter  Operator  Semigroups 

Mi  Zhou*  and  George  A.  Anastassiou 

The  University  of  Memphis 
Memphis,  Tennessee 

In  this  paper  we  expand  asymptotically  the  general  representation  formulae  for  {Co)m 
parameter  operator  semigroups.  When  we  consider  special  semigroups,  our  results  yield 
the  asymptotic  expansions  for  multivariate  Feller  operators.  In  particular,  the  asymptotic 
expansions  for  univariate  and  multivariate  Bernstein  operators  are  reobtained.  See  the 
related  examples  at  the  end. 


On  the  Approximation  of  Continuous  Functions  of 
Bounded  Variation  in  the  Space  C{R) 

V.  V.  Zhuk 

St. -Petersburg  State  University 
St.-Petersburg,  Russia 

Let  C  be  the  space  of  continuous  bounded  functions,  /  :  R  — R  with  the  norm 

OO 

ll/ll  =  sup  |/(a;)|,  V  the  set  of  functions  /  :  R  ^  R  with  the  variation  V  if)  <  oo,  and 
i€R 

a;(/,  h)  =  sup  ||/(-  +  t)  -  /II  the  modulus  of  continuity. 

|t|<a 

CX) 

Theorem  1.  Let  A,a,0>  0,  functions  f,geCnV,  f  dg  =  1.  Then  the  series 

—  OO 

OO 

:=  f(oo)g(oo)  -  f(-oo)g(-oo)  -  ^  {f{lk,p)  -  -  x)/q), 

k^—oo 
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where  'yk,0  sire  arbitrary  Sxed  points  of  the  segment  [kp,  (A:+ 1)^0]  converges  uniformly  over 
X  e  R  and  the  inequality 


11/ - «{/.?) II  <w(/,a^)  K(9)+u.(/,«  V  (9)  + 211/IK  V(9)  +  V^{9)} 

— i4  —00  —00  A 


holds. 

Corollary  1.  Let  A,q,P  >  0,f  €  C  DV.  Then 


RaMf)-=  ll/W-/(-°°)- 


(x-k/3)/a 

—00 


<  w(f,  oA)  +  a.(/, «  +  ^  y  dt. 

A 

If  K  is  the  class  of  functions  f  e  C  D  V  satisfying  the  conditions  sup  ||/||  < 

f€K 

00,  lim  sup  a>(/,  h)  =  0,  then  lim  sup  Ra  p{f)  =  0.  The  theorem  adduced  develops 
^-►o+  a,P^Q+  f^x 

in  the  several  directions  some  results  of  V.I.Zubov  (Dokl.Akad.Nauk  SSSR.  316.  (1991).  N 
6.  1298-1301)  concerning  to  the  approximation  of  probability  distributions.  The  analogs 
of  the  theorem  1  and  corollary  1  are  also  established  for  the  case  when  the  estimations  are 
carried  out  by  the  means  of  the  second  modulus  of  continuity. 
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